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FOREWORD 


This  report  was  prepared  by  the  Engineering  Technical  Materials  Group, 
Engineering  Department,  Convair,  A  Division  of  General  Dynamics  Corpo¬ 
ration,  San  Diego,  California  under  Contract  AF  33(6l6)-65o7,  "A  Study 
of  Structural  Design  Criteria,  Covering  Creep  Prediction  Methods  for 
Structural  Configuration.'1  This  contract  was  initiated  under  Project 
1367,  "Structural  Design  Criteria",  Task  No.  13534,  "Design  Criteria 
For  Heated  Airframes."  This  contract  was  administered  under  the  direc¬ 
tion  of  the  Flight  Dynamics  Laboratory,  Wright  Air  Development  Division, 
with  Lt.  J.  M.  O'Connor  and  later  Mr.  R.  M.  Bader  acting  as  project 
engineer . 


This  report  covers  the  period  of  work  from  April  15,  1959  to  April  15* 

I960. 

The  Convair  personnel  who  contributed  to  the  project  were:  Dr.  W.  H. 
Steurer,  Chief  of  Engineering  Materials,  Dr.  Carlo  Riparbelli,  Structural 
Research  Design  Specialist,  C.  W.  Alesch,  Engineering  Metallurgy  Design 
Specialist,  VI.  Robe,  Structural  Engineer,  G.  de  Vries,  Design  Engineer, 

}{.  H.  Stier,  J.  M.  Harvey  and  W.  M.  Parker,  Metallurgists,  and  A.  T. 

Green  and  L.  F.  Minkler,  Structural  Test  Engineers. 
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ABSTRACT 


A  step-by-step  numerical  process  for  introducing  thermal  and  creep  effects  into  a 
redundant  structural  analysis  was  developed.  This  process  involves  seven  consecu¬ 
tive  steps:  (1)  elastic  thermal  stress  determination;  (2)  unit  displacement  deter¬ 
mination  in  terms  of  stiffness  patterns;  (3)  relaxation  increment  determination; 

(,4J  real  thermal  stress  determination  resulting  from  summation  of  thermal  and 
relaxation  stresses;  (5)  incremental  stress  determination;  (6)  determination  of 
stress  redistribution  attendant  to  creep  and  (7)  summation  of  elastic,  thermal  and 
creep  stresses.  Calibrations  with  model  materials  and  structural  elements,  seeking 
verification  of  structural  criteria  introduced  into  numerical  processes  are  demon¬ 
strated.  An  example  problem  for  a  delta  wing  illustrates  method  of  calculation. 

A  critical  re-examination  of  interpretative  approaches  for  reducing  creep  data 
derived  from  materials  tests  into  simplified  form  suitable  for  numerical  analysis 
of  redundant  structures  is  introduced. 

The  creep  test  independent  variables,  time,  temperature  and  stress,  were  examined 
with  reference  to  their  influence  upon  the  dependent  variable,  strain.  These 
examinations  explored  three  fields  of  creep  activity  operative  under  conditions  of 
(1J  constant  stress  and  constant  temperature,  (2)  varying  stress  and  constant 
temperature,  and  (3)  constant  stress  and  varying  temperature.  Comparisons  of 
effects  arising  from  constant  stress  and  constant  temperature  creep  exposure  with 

those  arising  from  cyclic  stress  and  constant  temperature  exposures  also  are  pre- 
s  6nx e& • 


PUBLICATION  REVIEW 


This  report  has  been  reviewed  and  is  approved. 
FOR  THE  COMMANDER 
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PREDICTION  OF  CREEP  EFFECTS 
IN 

AIRCRAFT  STRUCTURES 
PART  I 

NUMERICAL  METHOD  FOR  INTRODUCTION 
OF  CREEP  EFFECTS  INTO  STRUCTURAL  ANALYSIS 


INTRODUCTION 

A  method  for  introducing  creep  factors  into  numerical  structural  analysis  is  presented. 
The  creep  effect  is  assumed  to  he  both  small  and  local  so  that  creep  deformation  may 
be  tolerated,  to  an  extent,  in  structures  exposed  to  the  combined  action  of  loads  and 
high  temperatures.  The  introduction  of  creep  factors  into  numerical  structural  anal¬ 
ysis  aims  primarily  to  extend  the  limits  of  velocity  and  acceleration  to  which  a  ma¬ 
terial  and  a  structure  built  from  it  can  be  submitted  for  a  limited  time. 

The  word  "creep"  as  used  herein  connotes  a  time  dependent  plastic  deformation  acti¬ 
vated  by  heat.  The  dependence  of  creep  upon  time,  a  characteristic  implied  in  all 
plastic  deformations,  is  essential  in  the  definition  of  creep.  Rigorously  speaking, 
time  is  always  a  factor  in  plastic  flow.  Practically  speaking,  creep  becomes  signi¬ 
ficant  when ^ the  rate  at  which  it  takes  place  is  appreciable  to  our  senses  and  measure¬ 
ments.  By  creep  rate  the  partial  derivative  of  the  plastic  deformation  with  res¬ 
pect  to  time  for  constant  stress"  is  meant,  but  this  does  not  mean  that  in  a  prac¬ 
tical  structure  the  stress  must  actually  remain  constant  in  conformance  with  mathe¬ 
matical  visualization.  Also,  since  creep  (time  rate  of  plastic  deformation  at  con¬ 
stant  stress)  and  relaxation  (time  rate  of  stress  decrease  at  constant  strain)  are 
practically  the  same  thing,  distinction  in  this  regard  is  not  made  herein.  In  fact 
plastic  deformation  is  introduced  into  the  numerical  process  in  the  form  of  a  relax¬ 
ation  term. 

The  nature  of  plastic  deformation,  including  creep  and  relaxation  is  presently  ex¬ 
plained  in  the  dislocation  theory  of  crystals  (l)(2)*,  and  its  discussion  is  beyond 
the  scope  of  this  report.  Its  essential  features  need,  however,  to  be  borne  in  mind 
for  a  rational,  even  if  crude  representation  needed  for  a  numerical  analysis. 

The  constitution  of  the  equilibrium  of  a  dissipative  system,  such  as  a  structure 
undergoing  deformation,  depends  at  any  time  upon  the  sequence  of  the  deformations 
involved  in  the  history  of  the  structure.  Accordingly  a  numerical  solution  for 
equilibrium  must  necessarily  be  found  for  successive  time  intervals.  The  effects  of 
various  loads  cannot  be  superimposed  irrespective  of  their  sequence  as,  by  way  of 
contrast,  can  be  done  in  a  linear  elastic  system. 


*  Numbers  in  parentheses  indicate  literature  references  listed  on  Page  10. 

Manuscript  released  by  the  authors  30  June  i960  for  publication  as  a  WADD 
Technical  Report. 
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If  the  plastic  strain  component  is  appreciable,  it  cannot  be  introduced  in  one  lump 
because  the  stresses  depend  on  the  plastic  deformation  and  further  plastic  deforma¬ 
tion  depends  on  the  stresses.  This  physical  fact  finds  its  counterpart  in  the  time 
dependence  of  the  plastic  flow.  Even  if  plastic  strain  of  appreciable  magnitude 
appears  within  a  negligible  time  compared  to  the  total  assumed  duration  time,  it  is 
still  necessary  to  divide  it  into  fractions  and  to  introduce  these  in  the  structural 
analysis,  taking  into  account  the  change  caused  in  the  stresses  at  each  step. 

BASIC  ASSUMPTIONS 


When  a  material  segment  of  finite  length  is  inscribed  on  a  structural  element,  this 
length  undergoes  a  change,  AL,  after  a  time  lapse.  At,  at  a  temperature,  T.  In  the 
elongation,  AL,  three  terms  can  be  distinguished:  a  thermal  expansion,  QT;  an  elas¬ 
tic  strain  component,  and  a  plastic  strain  component,  €p.  In  terms  of  unit 

elongation 


of  which: 


—  =  eP  +  e  +  OAT 
L  e  p 


(1) 


e 

e 


e 

P 

OAT 


represents  reversible,  time  independent, 
represents  irreversible,  time  dependent, 
represents  reversible,  time  independent. 


temperature  dependent 
temperature  dependent 
temperature  dependent 


elongation, 

elongation, 

elongation. 


and 


A  unit-elastic  elongation  is  related  to  stress  by  means  of  Hooke's  law  and  the  strain 
is  simultaneous  with  the  stress.  Plastic  strain  is  the  result  of  "slip"  at  certain 
surfaces  within  crystals  and  depends  on  the  history  of  shear  stress  at  these  slip 
surfaces.  Upon  application  of  load,  plastic  strain,  if  any,  follows  in  time: 


e 

P 


(2) 


and  it  is  a  deformation  at  constant  volume. 

In  a  region  of  a  structure  under  a  given  time  dependent  load  and,  in  which  the 
temperature  and  the  state  of  stress  are  also  known,  only  elastic  strain,  perfectly 
reversible,  will  be  observed  up  to  a  certain  limiting  stress.  The  limiting  stress 
is  dependent  on  the  state  of  stress,  on  the  temperature  and  on  the  past  history  of 
the  stresses.  When  the  stress  is  increased  above  the  limiting  stress,  or  what 
amounts  to  the  same  thing,  the  state  of  stress  is  maintained  and  the  temperature 
is  increased,  plastic  flow,  delayed  in  time  with  respect  to  the  stress  increment 
(3)  is  observed  to  take  place  at  a  rate  usually  decreasing,  and,  provided  the  stress 
is  constant,  this  decrease  proceeds  'until  it  becomes  negligible.  In  the  method  at 
hand  those  conditions  which  define  the  occurrence  of  negligible  creep  are  said  to 
define  the  points  of  creep  threshold  in  the  stress  plane:  0  (e^). 

If  stress  is  diminished,  under  conditions  of  creep,  the  elastic  strain  incorporated 
in  the  overall  deformation  diminishes  according  to  Hooke's  law.  The  plastic  defor¬ 
mation,  however,  does  not  reverse  appreciably  until  the  direction  of  the  shear  at 
the  slip  surfaces  is  deliberately  reversed.  This  view  omits  consideration  of  self 
induced  plastic  deformation  reversal  under  conditions  of  stress  diminution  (strain 
recovery) .  This  omission  is  made  in  view  of  a  need  for  a  simplified  representation. 
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THE  CREEP  THRESHOLD 


In  order  to  introduce  the  creep  effect  into  numerical  analysis,  the  usual  diagrams 
showing  -creep  deformation  with  respect  to  time  under  constant  stress  and  temperature 
conditions  are  used.  In  using  this  information,  data  in  the  region  of  primary  creep 
are  required,  but  under  current  conditions  are  seldom  available.  When  data  plots 
are  available  which  show  creep  rates  in  stress-strain  planes  for  given  temperatures, 
as  shown  in  Figure  1  which  is  extrapolated  from  data  given  in  (4),  they  also  may  be 
used.  In  this  approach  lines  are  traced  through  points  which  represent  stresses  at 
which  the  creep  rate  has  a  small  value  which  is  arbitrarily  called  negligible  and  is 
known  as  "threshold".  The  arbitrary,  negligible  creep  rate  utilized  in  Figure  1  is 
e  =  3  x  10"5  min"  .  The  "threshold"  is  defined  as  a  =  a  (e  ). 


OVERSTRESS 


Below  the  threshold,  creep  is  considered  negligible  and  is  set  equal  to  zero. 

Above  the  threshold  creep  occurs  at  a  rate  increasing  as  the  distance  of  points  of 
observation  recede  from  the  threshold.  In  plots  of  the  kind  shown  in  Figure  1  lines 
of  equal  creep  rate  can  be  interpolated  to  produce  a  diagram  of  the  kind  shown  in 
Figure  2  in  which  the  "overstress"  (a  -  a)  is  indicated  for  a  certain  stress  and 
strain  related  by  point  A  of  Figure  2,  for  example.  In  charts  of  the  kind  shown  by 
Figure  2  the  lines  of  equal  creep  rates  corresponding  to  equal  increments  in  the 
creep  rates  are  not  equally  spaced  because  the  creep  rate  is  not  a  linear  function 
of  the  overstress.  For  instance  the  representation 

ep  =  A  sinh  (cr  -  a)  (3) 

was  proposed  in  the  past  (5),  (6). 

VISCOSITY  REPRESENTATION  OF  CRE^P 

It  is  convenient  to  use  a  representation  of  creep  by  means  of  an  equation  of  the 
viscosity  type  (7)  (8)  such  as: 

£  =  0  for  |cr|<|c| 

ir 

(4) 

e  =  K  (0  -  a)  for  f  ct | >| ct | 

sr 

In  general,  K,  is  not  a  constant.  It  can,  however,  be  considered  a  constant  in  well 
behaved  materials  wherein  phase  changes  in  the  temperature  ranges  considered  are 
absent  and  the  overstress  and  strain  increments  are  small.  A  representation  accor¬ 
ding  to  a  viscosity  law  of  the  type  of  equation  (4),  obtained  by  linearizing  the 
true  creep  law  in  the  neighborhood  of  the  threshold  of  creep,  is  useful  for  compu¬ 
tation  because  it  offers  a  rational  basis  for  the  use  of  the  condition  of  plasticity 
in  dealing  with  creep. 

LIMITS  OF  USEFUL  RANGES  OF  PLASTIC  STRAINS 

The  plastic  strain  components  in  the  stress  and  deformation  analysis  described  here 
is  limited  to  small'  plastic  strains,  that  is  to  plastic  strain  components  of  the 
same  order  as  the  maximum  elastic  strain  that  can  be  imposed  at  the  same  region  of 
the  structural  element  considered.  For  instance  if  a  uniaxial  stress  of  180,000 
lbs/in  in  a  steel  member  is  allowed,  an  elastic  strain  of  6  x  10_3  follows  and  it 
may  be  decided  to  limit  analysis  to  plastic  deformations  in  the  order  of  1$  or  at 
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most  2 ’jo.  The  allowable  maximum  creep  rates  corresponding  to  the  Vjo  to  2 <jo  plastic 
strain  mentioned  will  be  determined  by  the  total  time  of  the  high  speed  and  high 
load  factor  phases  of  flight  during  which  creep  may  be  allowed.  For  instance,  if  a 
plastic  strain  of  1%  during  a  combat  phase  of  a  flight  of  30  minutes  duration  is 
allowed,  it  will  be  necessary  to  consider  those  combinations  of  stresses  and  temper¬ 
atures  which  will  produce  creep  rates  averaging  3  x  lO-^  min“l  throughout  aircraft 
life. 

The  reasons  for  such  limitations  are  the  following: 

1.  In  applications  involving  aircraft  structures,  plastic  strain  components  of  the 
same  order  as  the  elastic  ones  are  sufficient  to  cause  basic  changes  in  stress 
distributions  and  in  deflections.  This  range  of  strain  relations  is  therefore 
important  in  structural  analysis. 

2.  Aerodynamic  and  dynamic  considerations  limit  permissible  deformations  in  ordi¬ 
nary  flight  conditions. 

3.  Because  of  the  complexity  of  creep,  a  feasible  representation  of  its  action 
needs  to  be  limited  to  a  narrow  range,  essentially  to  a  fraction  of  the  primary 
creep  range,  with  exceptions. 

4.  No  plastic  component  of  strain  was  ever  allowed  previously  in  aircraft  struc¬ 
tures  within  limit  load,  except  a  residual  2  x  10-^  considered  negligible. 

This  convention  will  be  maintained  herein.  In  other  words  no  appreciable  plas¬ 
tic  deformation  is  assumed  under  conditions  wherein  plastic  deformation  occurs 
in  negligible  time,  especially  when  this  time  is  negligible  in  comparison  with 
that  assumed  for  creep  action.  On  the  same  ground  creep  rates  apt  to  produce 
plastic  deformations  of  2  x  10“3  strain  within  assumed  times  are  neglected. 
Consequently,  combinations  of  temperatures  and  states  of  stress  are  defined 
wherein  creep  rates  are  assumed  to  be  negligible.  Such  combinations  are  limi- 
ted  by  a  "threshold"  as  depicted  in  Figures  1,  2,  3  and  4. 

-'""''BIAXIAL  AND  TRIAXIAL  STATES  OF  STRESS 


While  most  available  test  data  pertaining  to  creep  are  obtained  under  uniaxial  and 
constant  stress,  biaxial  and  triaxial  strain  rates  need  be  introduced  into  numerical 
stress  analysis.  For  this  purpose  the  von  Mises  condition  of  plasticity  is  applied 
to  creep  according  to  Prager's  suggestion  to  furnish  a  working  assumption  (8),  (9) 
and  (10).  In  essence  the  von  Mises  condition  of  plasticity  defines  a  relation  among 
the  stress  components  for  which  plastic  deformation  will  appear.  All  that  is 
required  in  addition  to  this  is  the  amount  of  plastic  strain  which  can  be  considered 
non-negligible  and  how  long  stresses  may  act  before  creep  strains  evidence  them¬ 
selves  . 


Specifically,  when  the  principal  stresses  o^,  and  a ^  are  given,  the  value  of  the 
invarient  of  the  stress  tensor,  T,  is  computed  as  follows: 

T  =  (°i  “  °2^2  +  ^ °2  “  °3')2  +  (°3  “  ai^2  (5) 

2 

and  this  is  equated  to  its  value,  2a  ,  in  the  corresponding  uniaxial  stress  condi¬ 
tion.  Then  aQ  is  entered  in  a  diagram  as  shown  in  Figure  1  or  Figure  2,  and  the 
creep  rate,  £p0,  in  the  corresponding  uniaxial  condition  is  found.  In  addition  the 
hydrostatic  term  of  the  given  state  of  stress  is  computed  as  follows: 

P  =  1/3  {o±  +  a2  +  a  ) 
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along  with  the  components  of  the  stress  deviator 


VP 


S2  =  °2"P 


S3  =  °3~p 


(7) 


Then  since  the  power  dissipated  per  unit  volume  in  the  uniaxial  specimen  is 


S  =  c  •€ 
o  o 


(8) 


it  is  possible  to  proceed  under  the  following  working  assumptions t 

1.  The  power  dissipated  per  unit  volume  under  the  von  Mises-Prager  rule  is  the 
same  value  as  the  invarient  in  equation  ( 5 ) • 

2.  For  a  given  state  of  stress  the  strain  rate  components  are  such  that  the  power 
dissipated  is  a  maximum  (ll),  (12)  and  (13)*  If  the  deviator  components  are 
visualized  as  the  components  of  a  vector  and  the  strain  rate  components  as  the 
components  of  another  vector,  the  power  dissipated  is  a  maximum  if  the  scalar 
product  of  these  two  vectors  is  a  maximum.  Therefore  the  creep  rates  will  be 
proportional  to  the  components  of  the  deviator. 

These  two  conditions  are  written: 


Consequently 


S  =  s..  e  +  s^e^  +  s_e_  =  a  € 
11  22  33  oo 


s,  s„  s_ 

Jl _ 2  _  _3  _ 


=  7  =  const. 


2  ,  2  2 
sl  +  S2  +  s3 

?  =  — n — 1 

o  o 


(9) 

(10) 


(ii) 


from  which  the  creep  rates  are  determined. 
RELAXATION  OF  BIAX TALLY  STRESSED  PANELS 


In  the  case  of  biaxially  stressed  panels,  the  stress  normal  to  the  surface  is  zero. 
Assuming  the  principal  stresses,  and  a  ,  in  the  plane  of  the  panel  are  known  at 

a  certain  time,  while  =  o,  equation  (5)  is  written 


2  2 
°1  +  °2 


ala2  =  % 


(12) 


Then  for  a  given  value  of  one  of  the  two  principal  stresses,  for  instance  the  larger 
in  absolute  value,  any  plane  stress  distribution  may  be  represented  by 


a2  =  |a1  -1<  6  <  1 


(13) 
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Giving  recognition  to  particular  cases: 


|=1  a2  =  °1  Hydrostatic  plane 

stress 

|=o  0^=0  Uniaxial  stress 

|  =  -1  °2  =  -01  Pur6  shear 

equation  (12)  may  he  written: 

o±2  (1  +  I2  -  l)  =  oq2 

from  which 


(HO 

(15) 

(16) 

(17) 

(18) 


Equation  (l8)  relates  the  maximum  absolute  stress  in  the  biaxial  stress  distri¬ 
bution  to  the  uniaxial  stress,  o  ,  to  be  introduced  in  the  diagram  giving  the  creep 
rate. 


Equation  (l8)  is  plotted  in  Figure  5.  If  op  =  Op,  hydrostatic  plane  stress,  the 
condition  of  plasticity  is  the  same  as  if  op  =  o,  uniaxial  stress.  In  the  case  of 
pure  shear  ap  =  -op  =  C  the  condition  is  reached  for  a  value  of  the  shear  stress 
smaller  than  for  the  uniaxial  stress  by  a  factor  of  1.73* 

The  deviator  components  for  the  case  of  biaxial  stress  are  (equation  7) 


S1 

2 

"  3 

al  "  3  a2 

(19) 

s  _ 

2 

1 

(7  -  ~  0- 

(20) 

2 

“  3 

2  3  1 

S3 

l 

“  3 

(Op  +  ag) 

(21) 

and  in  general  all  three 

differ  from  zero. 

RELAXATION  OF  PANELS  RIGIDLY  HELD  AT  THEIR  EDGES 


In  this  case  the  assumption  is  made  that  the  plane  state  of  stress  is  uniform  in  a 
panel  and  that  its  components  are  known  and  expressed  in  terms  of  the  principal 
stresses  and  the  orientation  of  their  axes.  There  by  means  of  equation  (l8)  the 
equivalent  uniaxial  stress,  oQ,  is  computed.  This  stress  is  next  introduced  into 
the  experimental  diagram  shown  in  Figure  1  where  £pQ  is  obtained.  Following  this 

the  deviator  components  are  computed  from  equations  19,  20,  21  and  the  known  compo¬ 
nents  of  the  plane  state  of  stress.  Finally  from  equations  (10)  and  (ll)  the  creep 
rates  in  the  directions  of  the  principal  axes  in  the  plane  of  the  panel  are  com¬ 
puted.  These  creep  rates  are  assumed  to  remain  constant  during  a  time  interval.  At, 
which  is  chosen  to  be  sufficiently  short  for  the  assumed  condition  to  be  fulfilled. 
In  addition  the  temperature  is  also  assumed  to  be  constant  during  the  time  interval. 
At. 
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The  conditions  of  zero  elongation  in  both  directions,  denoted  subscript  1  and  2, 
give  rise  to 


(22) 


(23) 


which  when  solved  for  the  stress  increments  in  the  direction  of  the  principal  axes, 
result  in 


(24) 


(25) 


The  stress  increments  are  determined  by  the  increments  of  elastic  strain,  according 
to  Hooke's  law.  Poisson's  ratio  has  the  value  of  v  =  .3  as  in  most  materials  for 
the  application  of  Hooke's  law.  The  increments  of  plastic  strain  appear  in  equa¬ 
tions  (24)  and  (25)  inasmuch  as  they  are  equal  and  opposite  to  the  increments  of 
elastic  strain  in  the  relaxation  process  considered.  From  A and  An2  the  incre¬ 
ments  of  the  stress  components  in  the  direction  of  the  sides  of  the  panel  or  in  any 
other  chosen  direction  may  be  computed. 

The  relaxation  process  due  to  creep  is  an  essential  part  of  the  numerical  analysis. 
In  the  form  described  here  it  is  only  a  rational  guide  to  the  systematic  interpre¬ 
tation  of  test  data. 

SUBSTRUCTURES 

As  used  herein,  the  word  "substructure"  denotes  a  region  of  a  structure  limited  by 
an  arbitrarily  chosen  boundary  which  is  so  selected  that  the  stiffness  measured  on 
the  substructure,  taken  separately,  is  the  stiffness  of  the  real  structure  at  the 
region  considered.  A  substructure  is  composed  of  several  structural  elements  and 
its  stiffness  depends  on  the  stiffness  of  the  elements  which  compose  it  and  includes 
their  interaction,  local  instability,  connections,  etc.  For  example,  a  grid  of 
conveniently  spaced  planes  may  be  traced  by  means  of  lines  drawn  on  the  planform  of 
a  wing,  thus  defining  boxes  which  are  considered  substructures.  Each  box,  in  turn 
is  comprised  of  panels,  stringers  webs,,  etc.  In  the  assessment  of  this,  a  number 
of  boxes  can  be  built  and  tested  to  observe  the  displacement  versus  time  of  a  node 
or  of  a  section  with  respect  to  a  reference  plane  during  the  time  known  external 
loads  and  a  known  heat  input  are  imposed. 

NUMERICAL  ANALYSIS:  SYNTHESIS  APPROACH:  METHOD  OUTLINE 

The  case  of  a  structure  whose  mechanical  model  is  composed  of  single  elements  such 
as  stringers,  plates,  etc.  is  considered  first.  These  elements  are  joined  at  nodes 
at  which  internal  forces  equivalent  to  the  stresses  and  external  forces  equivalent 
to  the  given  loads  are  applied.  Consequently  the  structure  may  be  viewed  as  being 
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synthesized  from  elementary  components.  The  deformation  of  a  resultant  structure 
is  accordingly  represented  by  means  of  displacements  of  the  nodes.  The  states  of 
stress  and  the  temperature  in  each  element  are  represented  by  nominal  values  as  if 
they  vere  uniform  throughout  the  element . 

In  analysis,  the  stiffness  coefficients  for  each  element  are  computed  first  and  they 
depend  on  the  modulus  of  elasticity  of  the  material  used.  The  modulus  of  elasticity 
in  turn  depends  on  the  temperature.  The  stiffness  coefficients  relate  to  forces 
acting  on  elements  with  the  elastic  deformation  component  only.  Continuity  of 
structure  is  preserved  at  all  times  because  the  nodes  are  common  to  the  joined  ele¬ 
ments.  The  equations  used  in  analysis  fulfill  the  condition  of  equilibrium  and 
define  the  positions  of  the  nodes. 

The  stress  analysis  is  performed  by  a  step  by  step  process.  Each  step  of  the  analy¬ 
tical  process  covers  a  time  interval,  At,  and  consists  essentially  of  five  phases. 
Analysis  is  started  at  time,  t,  at  which  the  structure  is  in  equilibrium  and  its 
configuration  and  state  of  stress  are  known  at  all  elements.  The  analytical  steps 
are  as  follows: 


PHASE  I.  All  the  nodes  of  the  structure  are  locked  to  a  jig  of  infinite  rigi¬ 
dity.  The  temperature  increment,  AT,  corresponding  to  the  time 
interval,  At,  is  applied  at  each  structural  element.  The  corresponding 
"hydrostatic"  stress  terms  are  computed.  These  have  the  form 

Aa  =  -E  CAT  (26) 


for  bars  and  stringers  and  the  form 


Ac 

x 


Ea 

1-v 


AT 


(27) 


for  panels . 

Using  these  stress  terms,  the  stresses  acting  within  the  structural 
elements  are  computed  and  at  the  node  points  these  stresses  are  re¬ 
placed  with  equivalent  forces  acting  at  the  nodes.  Then  increments  of 
external  loads  are  applied  during  the  same  time  interval  used  in 
arriving  at  the  stress  terms  and  then  these  are  replaced  with  equiva¬ 
lent  concentrated  forces  at  the  nodes.  Finally  the  resultant  forces 
at  the  nodes  and  their  components  in  the  direction  of  chosen  coordinate 
axes  are  computed. 

PHASE  II.  Starting  with  the  conditions  established  in  Phase  I,  which  reflect  con¬ 
ditions  of  non-equilibrium,  the  node  displacements  necessary  to  reesta¬ 
blish  equilibrium  on  the  "unlocked"  structure  are  determined.  Following 
this  the  stresses  corresponding  with  the  displaced  node  conditions  are 
computed.  At  this  point  the  stress  analysis,  considering  thermal 
effects  but  neglecting  creep  in  the  time  interval  considered,  has  been 
performed . 

PHASE  III.  The  average  stresses  in  the  time  interval  considered  are  estimated 
according  to  creep  rates  which  are  recorded  on  calibration  diagrams 
derived  from  experiment .  The  creep  rates  are  interdependent  with 
stresses.  A  stress  increment  corresponding  to  a  convenient  fraction 
of  the  overstress  and  the  time  interval  consistent  with  this  overstress 
may  be  alternatively  chosen  for  stabilization. 
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In  this  step  it  is  necessary  to  take  into  account  the  "weight"  of  the 
section  areas  under  creep  versus  the  remaining  section  areas  in  order 
to  predict  the  change  of  stress  caused  by  creep. 

PHASE  IV.  The  structure  is  again  "locked"  at  the  nodes  and  then  the  force  incre¬ 
ments  corresponding  to  the  stress  relaxation  terms  over  the  whole  time 
interval  are  introduced. 

PHASE  V .  The  deformation  analysis  of  the  structure  is  repeated  as  in  Phase  II 

on  the  basis  of  the  force  increments  determined  in  Phase  IV.  The  cor¬ 
responding  stress  increments  due  to  relaxation  are  added  to  those 
found  in  Phase  II. 

The  process  deals  with  stress  relieving  increments  and  therefore  it  is 
inherently  convergent . 

NUMERICAL  ANALYSIS:  SUBSTRUCTURES  APPROACH:  DISCUSSION 

The  numerical  analysis  is  simpler  and  more  reliable  if  calibration  tests  directly 
obtained  in  substructures  are  available.  Then,  for  a  given  heat  input  per  time 
unit,  or,  for  a  given  boundary  layer  temperature,  direct  data  on  thermic  deforma¬ 
tion  and  on  creep  rates  which  depend  on  thermic  effect  and  on  the  applied  load  are 
provided.  In  this  situation  the  numerical  analysis  is  done  by  the  same  essential 
steps  already  outlined. 

The  numerical  process  described  above  takes  advantage  of  the  digital  computer  and 
is  preferred  because  of  the  labor  saving  advantage  gained  thereby.  Otherwise  a 
solution  could  be  obtained  by  more  laborious  relaxation  methods  which  introduce 
the  groups  of  internal  "forces"  with  zero  resultant  in  those  regions  where  creep 
occurs . 
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STRESS-LBS. /SQ.  IN. 


Figure  1.  Creep  Rates  at  Constant  Stress  Levels  (Schematic-See  Figures  23 
through  26  Numerical  Data) 
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SI  RAIN  RATE  PER  MINUTE 


Figure  2.  Creep  Rate  as  a  Function  of  Stress  and  Plastic  Strain  at  Constant 
Temperature  (Schematic) 
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STRESS 


Figure  3.  "Threshold"  and  Lines  of  Equal  Initial  Creep  in  the  Stress-Temp- 
perature  Plane  (Schematic) 
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STRESS 


Figure  4 


STRAIN 


"Threshold"  Surface  in  the  Stress-Strain  Temperature  Space 
(Schematic) 
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°0  UNIAXIAL  STRESS  IFnillVAIFNT) 
»  j  PRINCIPAL  STRESS  (MAXIMUM) 


Figure  5.  Condition  of  Plasticity  for  Biaxial  Stress 
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PREDICTION  OF  CREEP  EFFECTS 
IN 

AIRCRAFT  STRUCTURES 
PART  II 

CALIBRATION  METHODS  IN 
SUPPORT  OF  NUMERICAL  METHOD 


INTRODUCTION 

The  special  technique  specifically  aimed  to  determine  threshold  and  initial  creep 
rates  to  he  applied  to  panels  and  shell  elements  is  fully  described  in  reference 

(9). 

CALIBRATION  METHOD:  RESTRAINED  PANELS 

In  effecting  a  calibration,  a  panel  of  the  material  under  study  is  clamped  into  a 
frame  in  the  manner  shown  in  Figure  6.  The  material  used  in  the  clamping  frame  is 
different  from  that  under  study  and  is  so  selected  that  it  has  a  different  expan¬ 
sion  coefficient,  usually  smaller,  than  the  material  under  examination.  The  first 
calibration  step  involves  heating  the  clamped  assembly  in  an  oven  or  furnace  at  a 
predetermined  temperature  for  a  stated  time.  Since  the  panel  expands  at  a  greater 
rate  during  heating  than  the  clamping  frame,  it  is  put  under  compression  stress  when 
heated.  Providing  the  stresses  generated  are  sufficient,  creep  takes  place  in  the 
specimen  upon  heating.  Upon  completion  of  heating  the  assembly  is  allowed  to  slow 
cool  to  room  temperature  in  still  air. 

After  the  clamped  assembly  is  at  room  temperature,  strain  gages  are  applied,  as. 
illustrated  in  Figure  6,  and  the  bolts  which  attach  the  test  panel  to  the  clamping 
frame  are  released.  The  difference  in  strain  gage  readings  caused  by  release  of  the 
bolts  measures  the  stress  relaxation  which  occurs  from  removal  of  restraint  on  the 
panel.  The  stress  generated  in  a  panel  can  be  assigned,  within  limits,  by  choosing 
the  ratio  of  the  section  areas  of  the  panel  and  of  the  frame. 

The  residual  stress  for  a  given  constant  temperature,  and  therefore  constant  stress, 
has  the  form  of  a  coexponential  as  a  function  of  the  exposure  time.  This  is  shown 
in  terms  of  strain-time  variables  in  Figure  J.  The  coexponential  relationship  con¬ 
firms  the  assumption  of  viscosity  at  least  for  well  behaved  materials. 

In  these  tests  the  information  sought  consists  of  a  threshold  line  limiting  the  com¬ 
binations  of  force  or  moment  and  temperature  for  which  creep  is  negligible  as  illu¬ 
strated  in  Figure  8,  and  those  combinations  Of  forces  and  heat  input  for  which  creep 
rates  within  acceptable  ranges  are  observed.  Thus  calibration  diagrams  of  the  same 
type  as  those  shown  in  Figures  2  and  3  are  obtained. 

CALIBRATION  METHOD:  BEAMS  UNDER  PURE  BENDING 

In  these  tests  the  required  information  sought  consists  of  a  threshold  line  limiting 
totally  reversible  strain  in  a  beam,  the  onset  of  plastic  flow  and  lines  of  equal 
strain  rates  or  curvature  rates  for  combinations  of  temperature  and  bending  moments. 
The  growth  of  the  region  of  a  beam  which  undergoes  plastic  deformation,  the  amount 
of  plastic  deformation  in  those  regions  undergoing  plastic  flow  and  the  instability 
of  the  material,  if  any,  determine  the  curvature,  and  therefore,  the  strain  at  the 

lo 
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upper  and  lower  surfaces  of  the  beam  under  test.  In  the  results  from  team  tests  the 
difference  in  materials  properties  in  tension  and  compression,  and,  resultant  changes 
in  the  neutral  axis  are  implied.  The  neutral  axis  of  a  rectangular  section  should 
remain  at  the  centroid  of  the  section  providing  geometric  changes  do  not  occur.  In 
the  event  that  the  geometry  remains  constant  and  the  neutral  axis  shifts  during  test, 
the  indication  is  that  material  properties  changes  took  place.  To  determine  the 
rate  of  rotation  of  the  beam  with  respect  to  time,  incremental  strains  are  noted  as 
they  result  from  incremental  increases  in  moment.  Figure  9  schematically  diagrams 
the  bend  test  apparatus  used  for  tests  described.  Figure  10  displays  a  side  view 
of  a  test  beam  assembly  and  its  supports.  Figure  11  indicates  the  instrumentation 
and  support  provided  a  beam  test  specimen. 

CALIBRATION  METHOD:  SHEAR  PANELS 


Shear  panel  calibration  tests  are  made  in  conjunction  with  uniaxial  calibration  tests 
in  order  to  provide  correlation  data  representative  of  both  uniaxial  and  biaxial 
states  of  loading.  Uniaxial  loading  state  tests  are  made  with  flat  strips  with  the 
apparatus  pictured  in  Figure  12.  The  uniaxial- test  performance  consists  of  deter¬ 
mining  by  trial  the  loading  required  to  produce  0.0001  inch  per  inch  strain  in  ten 
minutes . 

The  biaxial  or  "pure  shear"  tests  are  made  with  square  specimens  mounted  in  a  steel 
shear  frame  as  shown  in  Figure  13.  Strains  are  measured  during  test  in  the  direc¬ 
tion  of  load  application  by  means  of  a  strain  gage  located  at  the  center  of  the 
square  plate.  In  order  to  avoid  buckling  of  the  test  plates,  restraining  plates  as 
shown  in  Figure  l4  are  lightly  pressed  upon  the  test  plates  and  serve  to  distribute 
loads  from  the  loading  connections  to  the  plate  in  addition  to  offering  the  desired 
buckling  restraint. 

In  the  analysis  of  data  from  the  biaxial  tests  described,  two  general  yield  criteria 
which  are  widely  accepted  are  those  of  Tresca  and  von  Mises.  The  mathematical 
statement  of  the  two  theories  are: 

Tresca 


a 

o 


(28) 


where  is  numerically  between  and  and 
von  Mises 

°0  =  (o1-o2)2  +  (o2-a3)2  +  (o3-o1)2^p 


In  the  case  of  pure  shear  where  a,  =  -aQ  and  a_  =  0,  the  two  conditions  become 


Tresca 


°o  =  2al 


von  Mises  a_  =  =  1.73a1 


(29) 

(30) 

(31) 


It  is  not  possible,  however,  to  evaluate  frictional  effects  of  the  plates  which  re¬ 
strain  the  shear  test  plates  from  buckling. 
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An  effective  strain  in  the  test  material  is  expressed  by 


de 


P 


+  (ae2 


+  (dU 


(32) 


by  using  the  fact  that  plastic  (creep)  action  takes  place  at  constant  volume  where 
de1  =-de2  and  de  =  0.  Accordingly,  in  the  shear  test,  the  effective  strain  in 
terms  of  strain  measured  in  the  direction  of  the  load  is 


de  =  -2—  de  =  1.15e  (33) 

0  TV  1  1 

In  order  to  check  the  concept  of  overstress,  where 

e  =  K  (o  -  oQ)  (3*0 

the  tensile  and  shear  specimens  are  subjected  to  a  load  ten  percent  greater  than  that 
required  to  cause  a  plastic  strain  of  0.0001  inch  per  inch  for  five  minutes.  Under 
this  condition  strain  increments  are  read  at  five  minute  intervals  and  average  creep 
rates  are  computed  from  these  as  functions  of  time  and  the  shapes  of  the  strain  curves 
for  the  uniaxial  and  biaxial  conditions  are  compared  to  note  the  differences  in  out¬ 
come  from  both  conditions  of  test. 

CALIBRATION  RESULTS:  RESTRAINED  PANELS 


Calibrations  employing  circular  restrained  plates,  although  planned,  were  not  accom¬ 
plished  because  of  material  procurement  diff iculties  which  arose  out  of  the  1959 
steel  strike. 

CALIBRATION  RESULTS:  BEAMS  UNDER  PURE  BENDING 


The  results  deriving  from  incremental  applications  of  moment  to  a  cantilever  beam  of 
the  kind  shown  in  Figure  9  are  given  in  Table  1.  The  beam  material  was  Armco  PH  15-7 
Mo  Condition  RH  950  stainless  steel.  The  beam  dimensions  were  0.125  inch  wide  by 
1.0  inch  deep  by  10  inches  long.  The  incremental  moment  loads  and  the  strains  at 
the  top  and  bottom  of  the  beams  were  measured  during  the  test  and  are  shown  in  Table 
1  along  with  the  calculated  neutral  axis  locations  and  radii  of  curvature.  These 
tabulations  show  the  tendency  for  the  neutral  axis  to  shift  as  loading  increases  and 
offers  some  evidence  of  non-uniformity  in  the  material's  reaction  to  tensile  and  com¬ 
pressive,  forces  .  Figures  15  and  l6  plot  the  results  from  a  second  series  of  bend 
tests  run  to  expand  upon  the  findings  of  the  initial  series.  These  curves  plot 
bending  moment  against  incremental  strain  and  show  a  distinct  cyclic  pattern  in  the 
stress-strain  relationships  obtaining  in  a  beam.  The  out  of  phase  relationship  of 
the  tension  and  compression  portions  of  the  beams  also  is  noteworthy.  These  circum¬ 
stances  suggest  action  taking  place  within  the  metal  itself  which  is  intimately 
associated  with  slip,  rotational  and  gliding  mechanisms  operative  within  the  metal¬ 
lic  lattice  structures.  Beam  tests  thus  present  themselves  as  a  useful  tool  for 
correlating  simultaneous  interactions  of  materials  and  possibly  structural  elements 
to  tensile  stresses,  which  ordinarily  are  measured  in  simplified  testing  practice, 
and  compressive  stresses,  which  usually  are  measured  under  conditions  of  experimental 
difficulty  and  complexity. 
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CALIBRATION  RESULTS:  SHEAR  PANELS 


,n  1^1  callb-yatlons  were  “a^e  with  copper,  lead  and  aluminum  shear  panels  and  in¬ 
tended  to  verify  the  usefulness  of  the  von  Mises  and  Tresca  criteria  for  study  of 
axiality.  Table  2  shows  data  obtained  from  tests  made  with  restrained  shear 
panels  in  the  manner  described  above.  In  performing  these  tests,  the  strain  and 
time  were  held  constant,  so  that  a  0.001  strain  was  attained  in  10  minutes,  while 
the  stresses  were  allowed  to  increase.  The  Table  2  data  sums  the  stresses  parallel 
and  transverse  to  the  loading  direction  for  differing  thicknesses  of  each  of  the 
three  materials  tested  and  also  shows  the  ratio  of  the  sums  of  the  directional 
stresses.  In  considering  the  outcome  from  the  individual  test  cases,  a  wide  varia¬ 
tion  in  the  stress  ratio  was  apparent.  This  immediately  suggested  the  existence  of 
varying  states  of  strain  arising  from  the  reaction  of  the  metal  itself  to  applied 
force  and  from  those  effects  arising  from  thickness  variations.  On  this  basis  veri¬ 
fications  of  the  two  theoretical  criteria  were  not  apparent.  However,  when  all 
stresses  for  all  materials  were  summed  and  ratioed,  the  "average"  result  fell  with¬ 
in  the  range  of  value  encompassed  by  the  two  theoretical  criteria. 


The  foregoing  experiment  emphasizes  a  need  for  judicious  selection  of  "model  mate¬ 
rials  for  demonstration  of  the  validity  of  theoretical  criteria  for  use  in  prac¬ 
tical  studies.  This  apparently  arises  from  strain  reactions  peculiar  to  a  base 
metal  subjected  to  stress.  However,  in  the  long  run,  the  need  for  a  suitable  model 
is  expressed  in  the  desirability  for  test  conditions  under  which  thermal  expansion 
and  buckling  effects  can  be  minimized  and  thus  bring  focus  of  attention  upon  the 
object  of  verification  rather  than  the  "side  effect". 


Tables  3  and  k  show  results  obtained  from  test  runs  of  30  minutes  duration  made  to 
examine  stress  and  strain  relations  in  tension  and  shear  panels.  The  test  data  them¬ 
selves  are  given  in  Figures  17  through  22.  In  both  Tables  3  and  4  a  wide  range  of 
I^rioni0n  in  effec^ive  Plastic  strain  and  in  strain  rate  is  reflected  in  a  summary  of 
the  30  minute  condition.  An  outstanding  characteristic  shown  in  Figures  17  through 
22  is  the  reaction  of  the  materials  to  load  application  as  is  indicated  by  the  varia¬ 
tions  in  ordinate  values  shown  in  the  several  charts. 
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TABLE  1  POSITION  OF  NEUTRAL  AXIS  AND  RADIUS  OF  CURVATURE 
OF  A  BEAM  UNDER  SIMPLE  BENDING 


Temperature 

Moment 
lbs . 

Strains 

Neutral  Axis 
Position 

Y 

t! 

Radius  of 
Curvature 

R 

Tt 

fcT 

n-j^  n 

B 

It-^  M 

CT+CB 

”1  !T 

75 

3000 

4o42x10" 

'b  343 5x10 “b 

74 77x10 "b 

.460 

133.74 

75 

3150 

4i4o 

3545 

7685 

.461 

130.12 

75 

3300 

4250 

3645 

7895 

.462 

126 . 66 

75 

3450 

4380 

3770 

8150 

.463 

122.70 

75 

3600 

4600 

3980 

8580 

.464 

116.55 

75 

3750 

5070 

4365 

9435 

.463 

105.99 

Material 

-  Armco 

PH  15-7  Mo,  Cond.  RH  950 

Symbols 

-  eT' 

strain  at  top  of  beam 

V 

strain  at  bottom  of  beam 

Yo> 

neutral  axis  distance 

2 7 

(see  sketch) 

R, 

radius  of  curvature  of  beam 

£T 

\ 

M 

N.  A  | 

/ 

Y2 

CB 

1 

B 

L-  1 

-** - '  JLU  '  r 
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TABLE  2 


RELATIONSHIP  OF  UNIAXIAL  TO  BIAXIAL  LOADING 


-  CRITERIA  VERIFICATION  DATA 


Material 

Thickness 

in. 

a 

0 

Stress  (l) 
lh/sq.  in. 

°1 

Q 

O 

Q 

H 

Copper 

.0216 

5,850 

M 

IT 

5,850 

-- 

-- 

Tt 

11 

— 

4,550 

-- 

tt 

— 

4,550 

— 

TT 

11,700 

8,830 

1.325 

Total 

tr 

.037 

10, 500 

_  _ 

_  _ 

IT 

TT 

11,1+50 

-- 

-- 

TT 

— 

5,000 

-- 

It 

— 

5,230 

— 

It 

21,950 

10,230 

2.140 

Total 

II 

.06J+7 

5,150 

— 

_  _ 

11 

IT 

5,200 

— 

Tt 

It 

2,650 

— 

TT 

TT 

_____ 

_  _ 

__ 

TT 

TT 

10,350 

— 

1.950 

Total 

Sum  of  Totals 

l+l+,000 

24,360 

1 . 806  ( Ave . ) 

Lead 

.0312 

1,125 

_  _ 

_  _ 

Tf 

Tt 

1,070 

-- 

-- 

TT 

11 

— 

835 

— 

TT 

TT 

— 

835 

— 

tt 

2,195 

i,64o 

1.338 

Total 

TT 

•  039 

1,000 

_  _ 

__ 

11 

tt 

1,000 

-- 

-- 

n 

IT 

— 

1,250 

-- 

IT 

— 

1,000 

-- 

IT 

TT 

2,000 

2,250 

0.888 

Total 

Tt 

.0652 

880 

_  __ 

TT 

IT 

835 

-- 

-- 

IT 

11 

— 

730 

— 

IT 

IT 

— 

675 

-- 

TT 

1,715 

1,405 

1.221 

Total 

Sum  of  Totals 

5,910 

5,295 

1.118  (Ave.) 
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TABLE  2  RELATIONSHIP  OF  UNIAXIAL  TO  BIAXIAL  LOADING 
(Cont'd)  -  CRITERIA  VERIFICATION  DATA 


Material 

Thickness 

in. 

a 

0 

Stress  (l) 
lb/sq.  in. 

°1 

Q 

O 

Q 

H 

Aluminum  (3) 

.020 

22,200 

— 

— 

tt 

tt 

19,500 

— 

— 

Tf 

tt 

— 

9,650 

— 

tt 

tt 

10,150 

— 

Tt 

tt 

31,700 

19,800 

1.601  Total 

tt 

.o4o 

18,000 

_  _  _ 

_ 

t! 

tt 

16,800 

— 

— 

tt 

tt 

— 

7,700 

— 

tt 

tt 

tt 

tt 

3M00 

15,400  (2) 

2 . 260  Total 

Sum  of  Totals 

66, 500 

35,200 

I.889  (Ave.) 

Sum  of  Totals  of  Stress 

for  All 

Materials 

O  = 

116, 410 

0 

°1  » 

64,855 

Vi 

1-795 

Value 

of  a  /a 

0  1 

from  von  Mises  criterion 

1-795 

Value 

of  a  /a 

0  1 

from  Tresca's  criterion 

2.000 

Notes : 


(1)  Stresses  given  are  those  in  the  direction  of  the  load  required 
to  produce  0.001  strain  in  10  minutes  where 

o  =  load/  2  x  thickness 
o 

=  load/  V"2"  x  length  of  side  x  thickness 

(2)  Twice  the  value  of  is  taken  as  the  total  since  only  one  shear 
test  was  made  for  the  material  and  thickness. 

(3)  7178-0  aluminum  alloy 
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TABLE  3  EFFECTIVE  PLASTIC  STRAINS 
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TABLE  4 


STRAIN  RATE  COMPARISON 


Material  Thickness 

Total  Strain 

Average  Strain 

Deviation 

Deviation  of 

in. 

in  30  min. 

in  30  min. 
n 

1  v 

Of  Strain 

Strain  + 

in .  x  10 

From 

Average 

(1) 

n  A  £i 
(2) 

Average 

Strain 

Copper 

.0216 

22 5T 

198 

27 

.136 

Copper 

.0216 

180T 

198 

-  18 

-.091 

Copper 

.0216 

190S 

198 

-  8 

-•043 

Copper 

.037 

2540T 

1468 

1072 

•  731 

Copper 

•  037 

2395T 

1468 

927 

.632 

Copper 

•  037 

585S 

1468 

-883 

-.602 

Copper 

•037 

425S 

1468 

-1043 

-.711 

Copper 

.0647 

43  5T 

348 

87 

.250 

Copper 

.0647 

250S 

3^8 

-  88 

-.250 

Lead 

.0312 

90  5T 

812 

-  93 

.114 

Lead 

.0312 

1380T 

812 

568 

.699 

Lead 

.0312 

150s 

812 

-662 

-.816 

Lead 

.039 

675T 

728 

-  53 

-.073 

Lead 

.039 

790T 

728 

62 

.085 

Lead 

•  039 

720s 

728 

-  8 

-.Oil 

Lead 

.0652 

44ot 

1014 

-574 

- .  566 

Lead 

.0652 

1588s 

1014 

574 

.  566 

Aluminum 

.020 

65T 

72.5 

-7.5 

-.135 

Aluminum 

.020 

80s 

72.5 

7-5 

•  135 

Aluminum 

.o4o 

21 5T 

222.5 

-7-5 

-.034 

Aluminum 

.o4o 

230s 

222.5 

7-5 

.034 

Notes : 

(l)  Suffix 

"T"  denotes 

tension  test: 

suffix  "S" 

indicates  shear 

test . 

(2)  Where  n  represents  material  thickness,  and  €.  represents 
incremental  strain.  1 
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POSITION  OF  STRAIN  GAGE 


l  1 
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L_ 

STEEL  RINGS  1 

k_ iri 
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1_ 1  Yl  ! 
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45" 

1 

i  h=.16"J 

^ _ _ _  ' 

i 

i 

- — i 

U - x- - ib 

5.2" 

SPECIMEN 


Figure  6.  Test  Specimen  Assembly 
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PLASTIC  STRAIN 


Figure  7.  Residual  Strain  in  Aluminum  Plates 
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BENDING  MOMENT 


Figure  8.  Creep  Bending  of  a  Beam  -  Lines  of  Equal  Time  Rate  of  Curva¬ 
ture  (Schematic) 
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Figure  9.  Beam  Test  Apparatus  (Schematic) 
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DYNAMOMETER 
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Figure  10.  Test  Beam  Assembly  and  Supports 
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Figure  12.  Uniaxial  Loading  State  Test  Apparatus 
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Figure  13.  Biaxial  Loading  State  Test  Apparatus 
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Figure  14.  Test  Panel  Instrumentation  and  Restraining  Plates  Used  in 
Biaxial  Tests 
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Figure  15.  Incremental  Strain-Bending  Moment  Relationships  in  Test  Beams 
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120 


□ - □ 


LOWER  GAGE -TENSION 
UPPER  GAGE-  COMPRESSION 


Figure  16.  Incremental  Strain-Bending  Moment  Relationships  in  Test  Beams 
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Figure  17.  Stress-Strain-Time  Relationships  for  Various  Thicknesses  of 
Copper  in  Tension 
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Figure  18.  Stress-Strain-Time  Relationships  for  Various  Thicknesses  of 
Copper  in  Shear 
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Figure  19.  Stress-Strain-Time  Relationships  for  Various  Thicknesses  of 
Lead  in  Tension 
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Figure  20.  Stress-Strain-Time  Relationships  for  Various  Thicknesses  of 
Lead  in  Shear 
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Figure  21.  Stress-Strain-Time  Relationships  for  Various  Thicknesses  of 
Aluminum  in  Tension 
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PREDICTION  OF  CREEP  EFFECTS 
IN 

AIRCRAFT  STRUCTURES 
PART  III 

DEMONSTRATION  OF  NUMERICAL 
METHOD 


PLAN  FOR  PROBLEM  SOLUTION 

The  plan  of  attack  for  problem  solution  must  necessarily  follow  a  time  sequence 
because  of .the  dissipative  nature  of  plastic  flow.  A  relaxation  pattern  was  used 
for  convenience  and  possibly  some  of  the  intermediate  outlined  phases  can  be  lumped 
together  in  a  new  numerical  scheme. 

In  a  visualization  of  the  sequence  of  physical  processes  identified  at  each  phase 
of  problem. attack,  the  following  events  are  pictured.  First,  a  structure  is  locked 
into  a  rigid  restraining  frame  at  some  uniform  reference  temperature,  say  room 
temperature.  The  structure  is  locked  at  node  points  and  is  in  a  state  of  equili¬ 
brium.  Then  the  structure  is  considered  heated  in  such  a  manner  that  temperature 
distributions  within  it  replicate  a  given  condition  of  structural  operation.  Dur¬ 
ing  heating  the  following  events  are  visualized  as  taking  place  in  two  successive 
steps:  an  accumulation  of  thermal  expansion  stress  followed  by  an  accumulation  of 
dilation  effects.  In  the  analysis,  stress  build-up  is  visualized  as  taking  place  in 
a  locked  structure,  and,  the  expansion  or  dilation  as  taking  place  upon  "unlocking" 
of  the  structure.  In  the  analysis  of  stresses  and  deformations  taking  place  under 
these  hypothetical  conditions,  the  structure  is  considered  an  elastic  unit  moving 
from  a  condition  of  equilibrium  at  low  temperature  to  a  new  condition  of  equilibrium 
dependent  upon  conditions  of  a  specified  elevated  temperature  distribution.  Finally 
in  completion  of  problem  solution,  readjustments  in  stresses  and  deformations  taking 
place  at  elevated  temperatures  as  a  consequence  of  creep  are  superimposed  upon  those 
stresses  and  deformations  which  result  from  elastic  action  taking  place  because  of 
thermal  expansion. 

In  detail,  the  attack  for  problem  solution  follows  the  sequence  of  steps  given 
below:  & 


1.  The  stresses  in  the  heated  and  locked  structure  are  computed. 

The  resultant  forces  and  moments  at  nodes  are  computed  from  stress  distribu¬ 
tions  determined  in  the  preceeding  step. 

3.  The  forces  and  moments,  so  far  as  they  are  carried  by  the  locks,  are  applied 
on  the  structure  as  external  loads.  At  this  stage  of  procedure  forces  and 
moments  are  treated  separately.  Moments  are  treated  according  to  existing 
routine  for  the  stress  analysis  of  wings.  Accordingly,  the  treatment  at  hand 
is  only  concerned  with  forces  acting  in  the  plane  of  a  wing. 

4.  The  stiffness  coefficients,  applicable  to  various  structural  elements,  such  as 
panels,  spar  caps  and  spar  webs,  are  computed  to  represent  temperature  condi¬ 
tions  under  which  individual  elements  operate. 

5.  The  stiffness  coefficients  for  each  node  of  the  structure  are  computed. 
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Using  the  loads  acting  upon  the  nodes,  which  were  determined  in  step  3,  and, 
the  stiffness  coefficients,  determined  in  step  5,  the  displacements  taking 
place  at  the  nodes  are  computed. 

The  stress -increments  arising  from  load  displacements  determined  in  step  6  are 
computed. 

8.  Stress  increments  found  in  step  7  are  superposed  upon  stresses  determined  for 
the  locked  structure  in  step  1.  These  are  thermal  stresses  present  in  a  per¬ 
fectly  elastic  structure  before  creep  is  introduced  to  it. 

9-  Using  thermal  stresses  determined  in  step  8,  equivalent  uniaxial  stresses  are 
determined  in  accordance  with  equation  (12)  given  under  Numerical  Method. 

These  equivalent  uniaxial  stresses  are  introduced  into  creep  charts  of  the  kind 
shown  in  Figures  35  to  38  inclusive  to  determine  the  possibility  of  creep  action 
being  present  in  various  structural  elements.  Where  creep  is  activated,  compo¬ 
nents  of  the  deviator  are  computed  and  corresponding  creep  rates  are  ascertained 
from  data  of  the  kind  shown  in  Figures  35  to  38  inclusive. 

10.  By  applying  simple  uniaxial  relaxation  to  one  section  and  by  proportioning 
relaxation  in  this  section  and  proportioning  then  relaxation  the  stress 
increments  in  terms  of  overstress,  a  suitable  time  interval  is  determined. 

Using  this  time  value,  relaxation  stress  increments  are  computed  in  terms  of 
principal  and  resolved  stress  increments. 

11.  Using  resolved  stress  increments,  corresponding  forces  at  the  locked  nodes  are 
determined  and  these  are  used  to  compute  corresponding  node  displacements. 
Finally,  stress  increments  due  to  these  displacements  are  computed. 

12.  The  incremental  stresses  computed  in  step  11  are  superimposed  upon  the  stresses 
determined  in  step  8.  These  totalized  stresses  are  the  thermal  stresses  in  the 
structure  and  reflect  elastic  and  creep  actions. 

EXAMPLE  PROBLEM  STATEMENT 

A  conventional  delta  wing  structure  is  assumed  for  this  example  and  is  shown  in 

Figure  23.  The  surfaces  are  equal  and  have  the  same  temperature  distributions.  The 

temperatures  assumed  are: 


Panels  450°F 
Spar  Caps  300^ 
Spar  Webs  150°F 
Front  Spar  (caps  and  webs)  450°F 


For  purposes  of  this  example,  the  structure  shown  in  Figure  24  has  its  node  [l] 
hinged  to  a  rigid  reference  frame.  The  nodes  [5],  [8]  and  [10]  are  hinged  to 
carriages,  so  that  displacements  in  the  "x"  direction  are  prevented  and  those  in 
the  "y"  direction  are  unrestricted.  In  proceeding,  the  method  of  stiffness  coeffi¬ 
cients  is  used  and  forces  acting  on  nodes  in  the  locked  structure  are  assumed  as 
input  for  solution  as  shown  in  Figure  25*  The  final  force  and  stress  distribution 
arrived  at  in  problem  solution  results  from  superposition  of  force  and  stress  incre¬ 
ments,  developed  in  the  course  of  redundant  analysis,  upon  those  initial  forces  and 
stresses  introduced  at  the  onset  as  shown  in  Figure  25- 


WADD  TR  60-411  PART  I 


44 


ASSUMPTIONS  INCORPORATED  IN  EXAMPLE  PROBLEM 


Assumptions  immediately  applicable  in  analysis  are: 

1*  Linear  relations  between  forces  and  displacements  exist.  Buckling  and  plastic 
deformations  are  excluded  from  the  general  assumption. 

Actions  of  forces  in  the  plane  of  the  wing  and  of  moments  bending  the  spars  are 
separately  computed.  Their  effects  are  additive. 

3.  The  deformation  of  the  entire  structure  in  its  middle  plane  is  described  in 
terms  of  displacement  of  its  nodes.  The  spar  elements  are  considered  as  bars 
hinged  at  nodes.  These  bars  remain  straight  between  nodes. 

4-.  Buckling  of  panels,  bending  of  spar  elements  between  nodes,  shear  lag,  and, 

effects  of  local  thermal  gradients  are  assumed  negligible  in  initial  analysis. 
These  assumptions  aim  to  obtain  average  stresses  in  the  structural  elements. 

In  the  event  consideration  need  be  given  these  conditions,  their  effects  may 
be  introduced  into  the  solution,  at  an  appropriate  point,  in  the  form  of  correc¬ 
tions  . 

EXAMPLE  PROBLEM  MATERIALS 


The  material  assumed  for  the  problem  is  7075-T6  aluminum  alloy.  The  value  of  Young's 
modulus,  E,  as  a  function  of  temperature  is  taken  from  NACA-TN-3584,  Figure  15,  page 
23.  The  coefficient  of  expansion,  also  a  function  of  temperature,  is  assumed  con- 
stant^but  is  used  only  in  the  presence  of  vertical  stress  and  strain.  The  value, 


a 


'  -  Zo  -j  *7  w-L  vcii/itai  ctliU.  Burain.  TJ 

12.8  x  10  F  ,  is  used  in  the  absence  of  more  precise  information. 


Poisson 's 


modulus  is  taken  at  v  =  .3  throughout  the  example.  Detailed  modulus  of  elasticity 
and  thermal  expansion  coefficient  data  are  shown  in  Appendix  II. 


EXAMPLE  PROBLEM  SOLUTION 


Initially  all  nodes  of  the  structure  are  "locked"  against  all  displacement  and  rota¬ 
tion  as  shown  in  Figures  24  and  25-  In  this  initial  condition  all  elements  of  the 
structure  are  at  a  uniform  temperature,  T0,  assumed  as  reference  temperature.  The 
temperature,  T±,  for  each  element,  i,  is  given  in  Table  5.  Temperature  increments 
are  computed  according  to  the  equation,  XT  =  Ti  -  TQ .  The  Young's  modulus,  E,  and 
the  expansion  coefficient,  a,  are  determined  as  functions  of  the  temperature,  T- 
according  to  the  respective  equations,  a  -  a  (Tp)  and  E  =  E  (T-^).  X> 

STRESSES  IN  THE  LOCKED  STRUCTURE 

The  axial  stress,  cr,  and  the  Force  N,  due  to  the  temperature  increment,  ATi,  in  each 
bar  (spar  cap)  having  cross  section  area.  A,  is  computed  while  the  bar  ends  are  res¬ 
trained  against  displacement  according  to  the  formulas,  a  -  -EoXT  and  N0  =  AEoAT. 

The  hydrostatic  stress,  a,  and  the  uniform  edge  thrust,  p,  in  each  plate  of  thick¬ 
ness,  t,  are  computed  according  to  the  formulas. 


0  =  0-,  =  o_  =  -  - -  oAT  and 

1  2  1  -  v 

pi  -  p2 =  -  rh  “r- 

These  thrust  components  are  the  force  inputs  for  the  stress  analysis.  They  are  re¬ 
duced  to  nodes  (panel  points)  in  the  usual  way  by  integrating  and  distributing  them 
to  adjacent  nodes  in  proportion  to  distances  to  the  nodes. 
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The  loads  to  the  locked  nodes  from  adjacent  panels  are  at  least  partially  balanced. 
The  resultant  forces  at  the  locked  nodes  are  the  force  inputs  for  the  stress  analy¬ 
sis  which  will  give  the  second  term  for  the  stresses. 

The  computed  stresses  and  forces  in  the  spar  elements  of  the  locked  structure  are 
given  in  Table  6  and  those  for  the  plate  elements  or  panels  are  given  in  Table  J. 

The  resultant  forces  acting  on  the  nodes  of  the  locked  structure  are  shown  in  Figure 
25  which  presents  the  force  input  for  solution  of  the  redundant  structure.  Figure  2 
augments  Figure  25  and  Table  7  in  that  it  diagrams  the  node  forces  balancing  hydro¬ 
static  pressure  in  triangular  panels  incorporated  in  the  delta  wing  plan  shown  in 
Figure  25* 

STIFFNESS  COEFFICIENTS  FOR  SINGLE  ELEMENTS 

In  a  spar  element  of  length,  L,  and  cross  section  area,  A,  which  is  fixed  at  its  left 
end  and  free  to  move  at  its  right  end,  as  shown  in  Figure  27,  the  stiffness  coeffi¬ 
cients,  measured  by  forces  acting  on  the  bar,  are  given  in  terms  of  horizontal  dis¬ 
placement,  U,  and  vertical  displacement,  V.  In  most  cases  vertical  displacements 
are  small  and  are  neglected  in  the  example  problem.  In  the  example  problem  stiffness 
coefficients  for  spar  and  rib  elements  in  the  "x"  and  "y"  directions  are  evaluated  in 
Figure  28 . 

In  consideration  of  oblique  frame  elements,  Figure  29  illustrates  a _ spar  element  of 
length,  L,  whose  projection  on  the  x-axis  is,  a,  and  on  the  y-axis  is,  b,  so  that 

L  =^\j( a2  +  b2).  When  a  unit  displacement  is  imposed  upon  node  B  in  the  x-direction, 
the  relation, 


DE  _  6L  a 

DB  "  1  L  > 

arises  from  similarity  of  the  triangles  ABC  and  EDB.  Therefore  the  strain  follows 
the  relation, 


5L  a 


and  the  axial  force,  N,  follows 

N  =  EAe  =  EA  . 


The  force  component  in  the  x-direction,  Nx> 

2 

N  =  N  f  =  EA 
x  L  l2 

and  the  force  component  in  the  y-direction, 


follows  the  relation 


N  follows  the  relation 

y 


N 

y 


k6 
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In  the  same  way,  when  unit  displacement  is  imposed  on  the  node  B  in  the  y-direction, 
the  axial  force  and  its  components  follow  the  relations, 


N  =  -EA 


75 


N  = 
x 


-EA 


ah 


and 


\  ^ 


The  stiffness  coefficients  for  the  oblique  bar,  shown  in  Figure  30,  which  arise  from 
the  relations  given  above  are  tabulated  in  Table  8  where  the  subscripts,  "B"  and  "C", 
denote  applicable  node  points.  In  addition,  the  stiffness  coefficients  for  the 

oblique  bar  can  be  written  analytically  as  it  follows  with  L 


8l 

3a 


2a 

2L 


a 

L 


-V 


(a  + 


b2) 


that 


and  for  a  displacement,  u,  of  the  node  B,  that  is  an  increment,  u,  of  length  a,  the 


relation 


6L  = 


3l 


u  = 


au 

L 


obtains,  and  therefore, 

6L 

e  =  _  = 


au 
T  2 


and,  for  an  incremental  displacement  of  u  =  1,  the  results  shown  in  Table  8  can  be 
derived.  The  numerical  results  obtained  for  the  example  problem  conditions  through 
use  of  the  stiffness  coefficients  given  in  Table  8  are  shown  in  Figure  31. 

The  stiffness  coefficients  for  the  rectangular  plates  are  computed  from  the  equations 
given  in  Table  9  which  are  taken  from  J.  H.  Argyris,  "Energy  Theorems  and  Structural 
Analysis",  Aircraft  Engineering,  April  1955>  page  125- 

The  stiffness  coefficients  for  the  triangular  plates  are  computed  from  the  equations 
given  in  Table  10  and  which  are  taken  from  Turner,  M.  J.,  Clough,  R.  W.,  Martin, 

H.  C.  and  Topp,  L.  S.,  "Stiffness  and  Deflection  Analysis  of  Complex  Structures", 
Journal  of  Aeronautical  Sciences,  September  1956.  The  numerical  stiffness  coef¬ 
ficients  applicable  to  the  triangular  plates  of  the  example  problem  are  given  in 
Table  11. 


The  computation  of  stiffness  coefficients  for  each  node  of  the  structure,  and  each 
displacement,  y,  taking  place  at  respective  nodes  is  summarized  in  Table  12.  The 
assembled  matrix  is  given  in  Table  13 .  The  degrees  of  freedom  up,  vp,  u^,  u q  and  upo 
at  the  node  points  [1]  [5]  [8]  and  [10]  are  set  equal  to  zero  according  to  the  boun¬ 
dary  conditions  established  and  shown  in  Figure  24.  The  loads  applied  are  those 
acting  on  the  nodes  of  the  locked  structure  as  diagrammed  in  Figure  25*  The 
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displacements  of  the  nodes  resulting  from  digitial  computer  calculations  for  the  load 
input  given  in  Figure  25  are  summarized  in  Table  14  and  are  indicated  in  Figure  31 
in  an  exaggerated  scale  vith  respect  to  that  used  to  diagram  the  structure. 

STIFFNESS  COEFFICIENTS  FOR  WING  STRUCTURE 


The  stiffness  coefficients  defined  for  unit  displacements  of  each  node  of  the  wing 
structure  in  the  directions,  u,  and  v,  are  computed  as  sums  of  the  stiffness  coeffi¬ 
cients  common  to  a  given  node.  The  assembled  matrix  of  stiffness  coefficients  is 
given  in  Table  13.  The  force  inputs  are  those  shown  in  Figure  25-  The  result  of 
combining  these  coefficients  with  the  force  Inputs  appears  in  terms  of  displacement 
of  nodes  as  shown  in  Table  l4  and  Figure  31« 

In  proceeding  to  wing  structure,  stress  increments  due  to  the  load  displacements  are 
computed.  The  internal  forces  are  next  determined  as  products  of  the  displacements 
and  the  element  stiffness  coefficient  matrix.  This  matrix  multiplication  is  per¬ 
formed  on  a  digital  computer  and  the  resulting  product  is  shown  in  Table  15  which 
lists  the  forces  on  rectangular  panels  or  plates  and  in  Table  1 6  which  lists  the 
forces  on  triangular  panels . 

Following  this,  stress  increments  for  panels  are  computed  from  the  forces  given  in 
Tables  15  and  1 6.  Final  stresses  are  found  by  superposition  of  stress  increments 
upon  stresses  in  the  locked  structure  shown  in  Table  5*  The  computations  for 
arriving  at  incremental  stresses  in  panels  of  the  wing  are  given  in  Table  17*  Com¬ 
putations  of  stress  increments  in  spar  caps  and  webs  are  shown  in  Table  18.  Final 
stresses  resulting  from  thermal  expansion  effects  and  effects  arising  from  displace¬ 
ments  attendant  to  thermal  expansion  are  given,  for  the  spar  caps  in  Table  19,  and, 
for  the  spar  webs  in  Table  20.  The  final  thermal  stresses  in  the  spar  caps,  spar 
webs  and  wing  panels  are  shown  in  Figures  32,  33  and  34  respectively.  Those  stresses 
shown  In  Figures  32,  33  and  3^  represent  thermal  stresses  in  perfectly  elastic 
structure  prior  to  introducing  effects  arising  from  creep  action. 

CREEP  EFFECT 


The  computation  of  the  creep  effect  proceeds  in  accordance  with  the  method  developed 
in  the  Numerical  Method  section  of  this  report.  The  creep  calculations  are  entered 
by  transforming  the  component  thermal  stresses  shown  in  Figures  32,  33  and  34  into 
equivalent  uniaxial  stresses  by  means  of  the  von  Mises-Prager  method  discussed  in 
the  Numerical  Method  section.  The  results  of  such  calculations  with  respect  to  the 
panels  in  the  example  problem  wing  are  given  in  Tables  21,  22  and  23  and  apply  at 
the  beginning  of  a  creep  period  where  time  is  considered  at  zero  value.  In  parti¬ 
cular,  the  uniaxial  stresses  sought  are  tabulated  in  column  4  of  Table  23.  When 
these  stresses  are  entered  into  creep  charts  of  the  type  shown  in  Figures  35  through 
39  inclusive,  the  creep  rates  shown  in  column  6  of  Table  23  are  developed  for  the 
temperatures  shown  in  column  5  of  the  table.  Then  returning  to  the  stress  components 
shown  in  Table  21  the  components  of  the  deviator  are  computed  in  accordance  with 
equations  (19),  (20)  and  (21)  of  the  Numerical  Method  section.  The  outcome  of  these 
computations  is  given  in  Table  24.  Employing  the  stress  deviator  values  given  in 
Table  24,  the  creep  rates  in  the  principal  directions  are  calculated  in  accordance 
with  equations  (10)  and  (ll)  of  the  Numerical  Method  section  to  develop  creep  rates 
corresponding  to  components  of  the  stress  deviator  as  shown  in  Table  25* 

A  numerical  value  for  the  time  interval,  At,  was  determined  as  shown  in  Table  26  by 
applying  a  simple  uniaxial  relaxation  to  one  section  and  by  imposing  the  condition 
that  the  relaxation  stress  increments  would  be  of  the  order  of  10%  of  the  overstress. 
The  time  interval  resulting  is  At  =  4  minutes.  The  relaxation  stress  increments 
in  the  locked  panels  of  the  wing  are  computed  as  shown  in  Table  27  in  terms  of  the 
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principal  stresses  A0j_  and  Ao'2*  The  corresponding  stress  increments  Acx,  Adj  and 
ATXy  are  calculated  as  shown  in  Tables  28  and  29* 

Continuing  with  the  calculations,  the  force  inputs  at  the  locked  nodes  are  computed 
for  the  rectangular  and  triangular  plates  as  shown  in  Tables  30  and  31  according  to 
formulas  given  in  respective  tables.  These  force  inputs  are  subsequently  resolved 
into  resultant  forces  (See  Tables  32  and  33)  acting  at  the  locked  nodes.  These 
resultant  forces  are  introduced  as  known  terms  in  the  system  of  equations  represented 
in  Table  13.  The  corresponding  displacements  obtained  by  combining  Table  31  and  32 
data  with  that  given  in  Table  13  through  matrix  multiplication  by  means  of  a  digital 
computer  are  found  in  Table  3^-*  The  increments  of  stresses  due  to  these  displace¬ 
ments  which  represent  relaxation  of  the  nodes  are  computed  as  shown  in  Tables  35 

to  38. 

The  total  stresses,  the  sum  of  the  initial  thermal  expansion  stress  given  in  Table 
5>  the  relaxation  increment  applied  to  the  thermal  expansion  stress  condition  given 
in  Table  16,  and  the  stress  increments  arising  subsequent  to  creep  as  given  in 
Table  38*  are  given  in  Table  39- 
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WADD  TR  60-411 


PART  I 


TABLE  5 


THERMAL  STRESSES  IN  EXAMPLE  PROBLEM  WING 


Element 

Temperature 

°F 

Modulus 

Of  Elasticity 
PSI  x  106 
(1) 

Thermal 

— 

Elements 

PSI 

(2) 

Stresses 

"Plate" 

Elements 

PSI 

(3) 

Panel 

450 

8.4 

-  68,400 

Cap  -  Spar 

300 

9.6 

-  36,900 

-- 

Web  -  Spar 

150 

9.8 

-  18,750 

-  26,700 

Caps  &  Webs 
(Front  Spar) 

450 

8.8 

-  50,700 

-- 

Notes: 

(1)  The  moduli  were  found  slightly  in  error.  Since  a  method  is  demonstrated 
corrections  are  not  made  to  correct  the  moduli  given. 

(2)  Where  the  thermal  stresses  in  "bar"  elements  are  computed  from  the  formula 

-  EOT. 

(3)  Where  the  thermal  stresses  in  "plate"  elements  is  computed  from  the  formula 

EOT  , 

-  - —  ,  and 
1-v  7 

E  represents  modulus  of  elasticity 
a  represents  coefficient  of  thermal  expansion 
T  represents  temperature  differential 
v  represents  Poisson's  ratio 
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TABLE  6  COMPUTED  SPAR  ELEMENT  STRESSES  AND  FORCES 
IN  EXAMPLE  PROBLEM  WING 


(1) 


Spar  Caps 


Area,  A,  =  1  sq.  in  where  each  cap  area  is  0.5  sq.  in. 
Temperature  Differential,  T,  =  300°F 

Axial  Stress,  a  =  -  EOT  =  -36,900  lhs/sq  in  (See  Table  5) 
Axial  Force,  Nq  =  -  36,900  lbs.  for  both  caps. 

Spar  Web 


Area,  A,  =  1  sq.  in.  where  depth  =  5  in.  and  thickness  =  .2  in 
Temperature  differential,  T,  =  150°F 

Axial  stress,  a  =  -EOT  =  -18,750  lbs/sq  in.  (See  Table  5) 
Axial  Force,  Nq  =  -18,750  lbs. 

Axial  Force  in  Locked  Spar 


Axial  Force,  N 
’  o 


-36,900  lbs.  (spar  caps)  -18,750  lbs.  (spar  web) 
-55,650  lbs 


Front  Spar 

Area,  A  =  2  sq.  in. 

Temperature  differential,  T,  =  4-50°F 

Axial  stress,  a  =  -EOT  =  -50,700  lbs/sq.  in.  (See  Table  5) 

Axial  Force,  N  =  2  cr  =  101, too  lbs. 

’  o  o  ’ 


Note: 


(l)  Horizontal  and  vertical  directions. 
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TABLE  7  STIFFNESS  COEFFICIENT  MATRIX  MULTIPLICATION  FOR  FRONT  SPAR 


IN  EXAMPLE  PROBLEM  WING. 


UB 

VB 

uc 

vc 

UB 

a2 

ab 

a2 

ab 

+  3 
L3 

+  3 

L3 

L3 

L3 

ab 

L3 

+  3 

L3 

bf 

L3 

uc 

2 

a 

ab 

2 

a 

ab 

L3 

L3 

+  3 

L3 

L3 

vc 

+  3 

L3 

_bf 

L3 

ab 

L3 

b2 
+  L3 
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TABLE  8  COMPUTED  RECTANGULAR  PLATE  STRESSES  AND  FORCES 


IN  EXAMPLE  PROBLEM  WING 


DATA 


Temperature  Differential,  T,  =  450  °F 
Modulus  of  Elasticity,  E,  =  8.4  x  10°  lbs/sq.  in 


Coefficient  of  Thermal  Expansion,  a, 
Thickness,  t,  =  0.1  in. 


12.8  x  10"6  in/in/°F 


SOLUTION 

Then  where  the  horizontal  stress,  o-^,  equals  the  vertical  stress,  Og, 

E  T 

u,  -  an  -  -  ^ —  =  -68, 400  lbs/sq.  in. 

1  2  1— v  7  ' 

and  the  edge  thrust,  p,  equals  -68,400  lbs/sq.  in. 

Hence  the  thrust  on  the  upper  edge,  Pu,  and  on  the  lower  edge,  P  ,  sums 
as 

Pp.  +  VI  -  136,800  lbs/sq.  in. 

Considering  a  10  in  (x-direction)  by  20  in.  (y-direction)  panel 

the  reaction  in  the  x-direction,  R^,  =  237>600  lbs,  and 

the  reaction  in  the  y-direction,  Fh,  =  138,600  pounds. 
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TABLE  9  STIFFNESS  COEFFICIENTS  FOR  RECTANGULAR  PLATES 
IN  EXAMPLE  PROBLEM  WING 


List  of  Coefficients 


K. 


13 


1  -  V 


E  d  •  t  G  •  t  •  i 

2  3f  6d 


K, 


23  1  -  V2 


E  a  •  t  G  •  t 

5d" 


K 


33  1  _  v2  31 


E  d  •  t  G  •  t 


3d 


K- 


E  d  •  t  G  •  t  •  i 


'43  ~  2  6i 

1  -  v 


3d 


K63 

K73 

k83 


v  E  t 


4(1  - 

v2) 

v  E 

t 

4(1  - 

v2) 

v  E 

t 

4  (1  • 

2 

-  V 

v  E 

t 

4(1  - 

v2) 

G  t 

TT 


G  t 

“T 


G  t 

"  T 


G  t 
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TABLE  9  STIFFNESS  COEFFICIENTS  FOR  RECTANGULAR  PLATES 
(Cont'd)  IN  EXAMPLE  PROBLEM  WING 


List  of  Symbols 

K  represents  the  stiffness  coefficient  where  the  first  digit  of  the 
subscript  denotes  the  node  point  at  which  load  application  is  con¬ 
sidered  and  the  second  digit  of  the  subscript  denotes  the  node  at 
which  displacements  are  under  consideration. 

E  represents  Young's  modulus. 

v  represents  Poisson's  ratio 

d  represents  dimension  (see  sketch) 

1  represents  dimension  (see  sketch) 

t  represents  thickness 

G  represents  shear  modulus 

Reference 

Argyres,  J.  H.,  "Energy  Theorems  and  Structural  Analysis",  Aircraft  Engineering, 
Vol  2,  No.  1)-,  April  1955*  page  125- 

Sketch 


(D 


(J) 
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TABLE  10  STIFFNESS  COEFFICIENT  MATRIX  FOR  TRIANGULAR  PLATES  IN  EXAMPLE  PROBLEM  WING 


CM 


-P 

W 


i — I 

OJ 


<0 

> 

d  Id 

CM 

> 

rd  I'd 

1 - 1 

0 

1 - 1 

> 

d  I'd 

+ 

d  |d 

d  I'd 

I 

-did 

CO 

:=> 

d  Id 

dH 

1 — 1 

1 

1 — 1 

0 

OJ 

!=> 

d 

0 

> 

1 

0 

> 

1 — 1 

t=> 

d 

1 — 1 

+ 

did 

^  I'd 

1 

d  Id 

1 — 1 

1 

CM 

d 

I — I 

1 

> 

1 

CM 

d 

CO 

!=> 

1 - 1 

> 

CM 

> 

CO 

> 

•H 

-p 
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0 

o 

H 

I: 

w 

o 

o 


03 

03 


ft 

Eh 

03 


cci 

-P 

P 

O 

N 

•H 

u 
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i — i 

cd 

o 

•H 
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> 
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1 
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H  <p 
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•H 
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ft 
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-p 
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p 
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.P 

to 

p 
<D 
i — i 

• — i 
p 

to 

CD 

e 

oj 
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p 
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-P 

-P 

4-> 

-P 

d 

d 

d 

d 

0 

0 

0 

0 

0 
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d 
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reference:  Turner,  M.  J.,  Clough,  R.  W.,  Martin,  H.  C.  and  Topp,  L.  S.  "Stiffness  and  Deflection 
Analysis  of  Complex  Structures",  Journal  of  Aeronautical  Sciences,  Vol.  2^  No  0 
September  1956,  page  805.  '  y’ 


TABLE  11  STIFFNESS  COEFFICIENTS  FOR  TRIANGULAR  PLATES  IN  EXAMPLE  PROBLEM  WING 


VO 

o 


x 


on 

> 

b- 

b- 

OJ 

0 

1 

b- 

b- 

OJ 

O 

0 

OJ 

VO 

-4" 
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on 
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-4" 

> 

on 

O 

on 

VO 

VO 

O 

• 

• 

• 

• 

0 

1 

0 

0 

1 

0 

b- 

on 

00 

VO 

OJ 

b- 
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OJ 

-4" 

-4* 

c— 

CO 

CO 

0 

eg 

0 

OJ 

« 

• 

1 — 1 

1 

1 — 1 

O 

1 

0 

in 

in 

b- 

VO 

1 — l 

t*- 

JZ) 

0 

-4" 

vo 

on 

0 

CO 

1 — 1 

VO 

OJ 

OJ 

• 

• 

• 

• 

• 

OJ 

1 — 1 

1 

O 

l 

0 

on 

1 

0 

1 

OJ 

on 

1 — 1 

eg 

on 

to 

to 

to 

t> 

> 

> 

WADD  TR  60-411  PART  I 


58 


TABLE  11  STIFFNESS  COEFFICIENTS  FOR  TRIANGULAR  PLATES  IN  EXAMPLE  PROBLEM  WING 
(Cant'd) 
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represents  thickness 


TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  1 


STIFFNESS  COEFFICIENTS  x 

:  106 

NODE 

ul 

DISPLACEMENT 

V1 

DISPLACEMENT 

STRINGER 

PANEI. 

1256 

TOTAL 

STRINGER 

PANEL 

1256 

TOTAL 

U1 

1.94 

1.337^ 

3 . 2744 

— 

0.3 

0.3 

V1 

— 

0.3 

0.3 

0.97 

0.7382 

1 . 7082 

U2 

-1.94 

-1.1758 

-3.1158 

— 

0.0231 

0.0231 

V2 

— 

0.0231 

-0.0231 

— 

-0.277 

-0.277 

U5 

— 

0.5072 

0 . 5072 

— 

-0.0231 

-0.0231 

V5 

— 

0.0231 

0.0231 

-0.97 

-0.0921 

-1.0621 

u6 

— 

-0 . 6688 

-0.6688 

— 

-0.3 

-0.3 

v6 

— 

-0.3 

-0.3 

— 

-O.369I 

-O.369I 
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SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


TABLE  12 
(Cont 'd) 


NODE  2 


NODE 

STIFFNESS  COEFFICIENTS  x 

106 

U2  DISPLACEMENT 

V2  DISPLACEMENT 

STRINGER 

PANEL  PANEL 
1256  2367 

TOTAL 

STRINGER 

PANEL 

1256 

PANEL 

2367 

TOTAL 

U1 

-1.94 

-1.1758  ~~ 

-3.1158 

— 

-0.0231 

-0.0231 

V1 

— 

0.0231  - 

0.0231 

— 

-0.277 

— 

-0.0277 

U2 

3-88 

1.3374  1.3374 

6.5548 

— 

-0.3 

0.3 

0 

V2 

— 

-0.3  0.3 

0 

0.97 

O.7382 

0.7382 

2.4464 

U3 

-1.94 

-  -1.1758 

-3.1158 

— 

— 

0.0231 

0.0231 

V3 

— 

-  -0.0231 

-0.0231 

— 

— 

-0.277 

-0.277 

U5 

— 

-0 . 6688  - 

-0.6688 

— 

0-3 

— 

0-3 

V5 

— 

-0.3  - 

0.3 

— 

-O.369I 

— 

-O.369I 

u6 

— 

O.5072  O.5072 

1.0144 

— 

0.0231 

-0.0231 

0 

V6 

— 

-0.0231  0.0231 

0 

-0.97 

-0.0921 

-0.0921 

-1.1542 

U7 

— 

-  _0 . 6688 

-0 . 6688 

— 

- - 

-0.3 

-0.3 

V7 

.... 

-  -0.3 

-0.3 

— 

— 

-O.369I 

-0.3691 

WADD  TR  bO-4ll  PART  I 


6l 


TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  3 


NODE 

STIFFNESS  COEFFICIENTS  x  10 

U  DISPLACEMENT 

3 

DISPLACEMENT 

STRINGER 

PANEL 

2367 

PANEL 

347 

TOTAL 

STRINGER 

PANEL 

2367 

PANEL 

347 

TOTAL 

U? 

-1.94 

-1.1758 

-3-1158 

— 

-0.0231 

— 

-0.0231 

V2 

— 

0.0231 

— 

0.0231 

— 

-0.277 

— 

-0.277 

U0 

3-88 

1.3374 

2.0075 

7.2249 

— 

-0.3 

0.6 

0.3 

3 

V3 

— 

-0.3 

0.6 

0.3 

•  97 

0.7382 

1.108 

2.8162 

U) 

-1.94 

_ 

-1.846 

-3.786 

— 

— 

-0.277 

-0.277 

4 

V4 

— 

— 

-0.323 

-0.323 

— 

— 

-0.646 

-0.646 

-0 . 6688 

______ 

-0.6688 

— 

0.3 

— 

0.3 

6 

V6 

— 

0.3 

— 

-0.3 

— 

-O.369I 

— 

-0.3691 

Ury 

0 . 5072 

-0.1615 

-0.3457 

— 

0.0231 

-0.323 

-0.2999 

7 

v7 

— 

-0.0231 

-0.277 

-0.3001 

-0.97 

-  .0921 

-0.462 

-l . 5241 
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TABLE  12  SUMMARY  OF  STRIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  4 


NODE 

STIFFNESS  COEFFICIENTS  x  10 

DISPLACEMENT 

DISPLACEMENT 

STRINGER 

PANEL 

347 

TOTAL 

STRINGER 

PANEL 

347 

TOTAL 

U3 

-1.94 

-1.846  ■ 

-3.786 

— 

-0.323 

-0.323 

V3 

— 

-0.277 

-0.277 

— 

-0 . 646 

-0.646 

U4 

2.114 

1.846 

3.960 

-0.348 

0 

-0.348 

V4 

-0.348 

0 

-0.348 

0.696 

0.646 

1.342 

°7 

-0.174 

0 

-0.174 

0.348 

0.323 

0.671 

V7 

0.348 

O.277 

0.625 

-O.696 

0 

-0.696 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  5 


NODE 

STIFFNESS  COEFFICIENT  x  10 

U_  DISPLACEMENT 

5 

V5 

DISPLACEMENT 

STRINGER 

PANEL  PANEL 

1256  5689 

TOTAL 

STRINGER 

PANEL 

1256 

PANEL  TOTAL 

5689 

U1 

_ 

0.5072 

0 . 5072 

— 

0.0231 

0.0231 

V1 

— 

-0.0231 

-0.0231 

-0.97 

-0.0921 

-1.0621 

u? 

_ 

-0 . 6688  — 

-0.6688 

— 

0.3 

0.3 

c. 

V2 

— 

0.3 

0.3 

— 

-O.369I 

—  -O.369I 

U5 

1.94 

1.337^  1-3374 

4.6148 

— 

-0.3 

0.3  0 

\ 

— 

-0.3  0.3 

0 

1.94 

0 . 7382 

0.7382  3.4i64 

UA 

-1.94 

-1.1758  -1.1758 

-4.2916 

_ _  _ 

-0.0231 

0.0231  0 

v6 

— 

0.0231  -0.0231 

0 

— 

-0.277 

-0.277  -0.554 

UA 

V8 

0.5072 

0 . 5072 

_ _ 

— - 

-0.0231-0.0231 

— 

0.0231 

0.0231 

-  -97 

— 

-0.0921-1.0621 

n 

-0.6688 

-0.6688 

_____ 

_ 

-0.3  -0.3 

9 

V9 

— 

-0.3 

-0.3 

—  —  —  — 

“  “  “ 

-0.3691-0.3691 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  6 


NODE 

STIFFNESS  COEFFICIENT  x  10^ 
HORIZONTAL  (U^)  DISPLACEMENT 

STRINGER 

PANEL 

PANEL 

PANEL 

PANEL 

TOTAL 

1256 

2367 

5689 

679 

U1 

— 

-0 .6688 

— 

— 

- - 

-0 . 6688 

V1 

— 

-0.3 

— 

— 

— 

-0.3 

U2 

— 

0.5072 

0.5072 

— 

— 

1.0144 

V2 

— 

0.0231 

-0.0231 

— 

— 

0 

U3 

— 

— 

-0 . 6688 

— 

— 

-0.6688 

V3 

— 

— 

0.3 

— 

— 

-0.3 

U5 

-1.94 

-1.1758 

— 

-1.1758 

— 

-4 . 2916 

V5 

— 

-0.0231 

— 

0.0231 

— 

0 

u6 

3.88 

1-3374 

1-3374 

1-3374 

2.0075 

9.8997 

V6 

— 

0.3 

-0.3 

-0.3 

0.6 

0.3 

°7 

-1.94 

— 

-1.1758 

— 

-1.846 

-4 . 9618 

V7 

— 

— 

0.0231 

— 

-0.323 

-0.2999 

U8 

— 

— 

— 

-0.6688 

— 

-0 . 6688 

V8 

— 

— 

— 

0.3 

— 

0.3 

°9 

— 

— 

— 

0 . 5072 

-0.1615 

0.3457 

V9 

-0.0231 

-0.277 

-0.3001 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  6 


NODE 

STRINGER 

STIFFNESS  COEFFICIENT  x  10 
VERTICAL  (vg)  DISPLACEMENT 

PANEL  PANEL  PANEL  PANEL 

1256  2367  5689  679 

TOTAL 

U1 

— 

-0.3 

— 

— 

— 

-0.3 

V1 

— 

-O.369I 

— 

— 

— 

-O.369I 

U2 

— 

-0.0231 

0.0231 

— 

— 

0 

V2 

-0.97 

-0.0921 

-0.0921 

— 

— 

-1.1542 

u3 

— 

— 

0.3 

— 

— 

-0.3 

V3 

— 

— 

-0.3691 

— 

— 

-O.369I 

U5 

— 

0.0231 

— 

0.0231 

— 

0 

V5 

— 

-0.277 

— 

-0.277 

— 

-0.554 

U6 

— 

0.3 

-0.3 

-0.3 

0.6 

0.3 

V6 

1.94 

0 . 7382 

0.7382 

0 . 7382 

1.108 

5 . 2626 

U7 

— 

— 

-0.0231 

— 

-0.277 

-0.3001 

V7 

— 

— 

-0.277 

— 

-0.6A6 

-O.923 

U8 

— 

— 

— 

0.3 

— 

0.3 

V8 

— 

— 

— 

-0.3691 

— 

-0.3691 

U9 

— 

— 

— 

0.0231 

-0.323 

-0 . 2999 

V9 

-0.97 

-0.0921 

-0.462 

-l . 5241 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  7 


NODE 

c 

STIFFNESS  COEFFICIENT  x  10 
HORIZONTAL  (U7)  DISPLACEMENT 

STRINGER 

PANEL 

2367 

PANEL 

347 

PANEL 

679 

TOTAL 

U2 

— 

-0.6688 

— 

— 

-0.6688 

V2 

— 

-0.3 

— 

— 

-0.3 

U3 

— 

0.05072 

-O.I615 

— 

0.3457 

V3 

— 

0.0231 

-0.323 

— 

-0.2999 

U4 

-0.174 

— 

0 

— 

-0.174 

V4 

0.348 

— 

0.323 

— 

0.671 

U6 

-1.94 

-1.1758 

— 

-1.846 

-4.9618 

V6 

— 

-0.0231 

— 

-O.277 

-0.3001 

U7 

2.288 

1.3374 

0.1615 

1.846 

5.6329 

V7 

-0 . 696 

0.3 

0 

0 

-0.396 

U9 

-0.174 

— 

— 

0 

-0.174 

V9 

0.348 

0.277 

0.625 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  7 


NODE 

£ 

STIFFNESS  COEFFICIENT  x  10 
VERTICAL  (V^)  DISPLACEMENT 

STRINGER 

PANEL 

2367 

PANEL 

347 

PANEL 

679 

TOTAL 

u2 

_____ 

-0.3 

_ 

— 

-0.3 

V2 

— 

-O.369I 

— 

— 

-O.369I 

u. 

— 

-0.0231 

-0.277 

— 

-0.3001 

j 

V3 

-0.97 

-0.0921 

-0.462 

— 

-1 . 5241 

u4 

0.348 

— 

0.277 

— 

0.625 

\ 

-0.696 

— 

0 

— 

-0.696 

u6 

— 

0.0231 

— 

-0.323 

-0.2999 

V6 

— 

-0 . 277 

— 

-0.646 

-0.923 

UY 

V7 

-0 . 696 

0.3 

0 

0 

-0.396 

2.362 

0.7382 

0.462 

0.646 

4.2082 

U9 

0.348 

— 

— 

0.323 

0.671 

y 

V9 

-O.696 

0 

-0.696 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  8 


NODE 

STIFFNESS  COEFFICIENT  x  106 

UQ  DISPLACEMENT 

V 

Q  DISPLACEMENT 

STRINGER 

PANEL 

5689 

PANEL 

8-9-10 

TOTAL 

STRINGER 

PANEL 

5689 

PANEL 

8-9-10 

TOTAL 

U5 

— 

0 . 5072 

— 

0 . 5072 

— 

0 . 0231 

— 

0.0231 

V5 

— 

-0.0231 

— 

-0.0231 

-0.97 

-0.0921 

— 

-1.0621 

u6 

— 

-0 . 6688 

— 

-0 . 6688 

— 

0-3 

— 

0.3 

V6 

— 

0.3 

— 

0.3 

— 

-O.369I 

— 

-0.3691 

U8 

1.94 

1.3374 

2.0075 

5 . 2849 

— 

-0.3 

0.6 

0.3 

V8 

— 

-0.3 

0.6 

0.3 

1.94 

0.7382 

1.108 

3 • 7862 

U9 

-1.94 

-1.1758 

-1.846 

-4 . 9618 

— 

-0.0231 

-0.277 

-3.001 

V9 

— 

0.0231 

-0.323 

-0.2999 

— 

-0.277 

-0.646 

-0.923 

O 

t — 1 

1=) 

— 

— 

-0.1615 

-0.1615 

— 

— 

-0.323 

-0.323 

c 

H 

O 

-0.277 

-0.277 

-0.97 

-0.462 

-1.432 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  9 


NODE 

STIFFNESS 

HORIZONTAL 

COEFFICIENT  x  106 
(Ug)  DISPLACEMENT 

STRINGER 

PANEL 

5689 

PANEL 

679 

PANEL 

8-9-10 

TOTAL 

U5 

— 

-0 . 6688 

— 

— 

-0 . 6688 

V5 

- - 

-0.3 

—  —  —  — 

-0.3 

U6 

— 

0.5072 

-0.l6l5 

— 

0.3457 

V6 

— 

0.0231 

-0.323 

— 

-0.2999 

U7 

V1 

-0.174 

— 

0 

— 

-0.174 

0.348 

— 

0.323 

— 

0.671 

U8 

-0.194 

-1.1758 

— 

-1.846 

-4.9618 

0.0231 

— 

— 

-0.277 

-0.3001 

uo 

2.288 

1.3374 

0.1615 

1.846 

5.6329 

y 

V9 

-0.696 

0.3 

0 

0 

-0.396 

U10 

-0.174 

— 

— 

0 

-0.174 

V10 

0.348 

0.277 

0.625 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  9 


NODE 

xr 

STIFFNESS  COEFFICIENTS  x  10 
VERTICAL  (Vg)  DISPLACEMENT 

STRINGER 

PANEL 

5689 

PANEL 

679 

PANEL 

8-9-10 

TOTAL 

U5 

— 

-0.3 

— 

— 

-0.3 

V5 

— 

-O.369I 

— 

— 

-O.369I 

u6 

— 

-0.0231 

-0.277 

— 

-0.3001 

V6 

-0.97 

-0.0921 

-0.462 

— 

-1.5241 

U7 

0.348 

— 

0.277 

— 

0.625 

V7 

-O.696 

— 

0 

— 

-O.696 

U8 

— 

0.0231 

— 

-O.323 

-0.2999 

V8 

— 

-0.277 

— 

-0.646 

-0.923 

U9 

-0.696 

0.3 

— 

0 

-0.396 

V9 

2.362 

O.7382 

0.462 

0.646 

4 . 2082 

U10 

0.348 

— 

— 

— 

0.671 

V10 

-O.696 

-O.696 
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TABLE  12  SUMMARY  OF  STIFFNESS  COEFFICIENTS  AND  DISPLACEMENTS 
(Cont'd)  AT  EACH  NODE  OF  EXAMPLE  PROBLEM  WING 


NODE  10 


NODE 

STIFFNESS  COEFFICIENTS  x 

106 

U10 

DISPLACEMENT 

V  DISPLACEMENT 

STRINGER 

PANEL 

8-9-10 

TOTAL 

STRINGER 

PANEL 

8-9-10 

TOTAL' 

U8 

V8 

-0.l6l5 

-0.l6l5 

— 

-0.277 

-0 . 277 

— 

-0.323 

-0.323 

0.97 

-0.462 

-1.432 

uQ 

-0.174 

0 

-0.174 

0.348 

0.277 

0.625 

9 

V9 

0.348 

0.323 

0.671 

-0.696 

0 

-0.696 

U10 

V10 

0.174 

0.1615 

0.3355 

-0.348 

0 

-0.348 

-0.348 

-0.348 

1.666 

0.462 

2.128 

WADD  TR  60-4ll 


PART  I 


72 


TABLE  13  STIFFNESS  COEFFICIENT  MATRIX  FOR  EXAMPLE  PROBLEM 

SOLUTION  INPUT 

NODES  1,  2  &  3 


STIFFNESS  ( 

COEFFICIENTS 

x  106 

DISPLACE¬ 

MENT 

U1 

vi 

U2 

V2 

U3 

T3 

U1 

3-2744 

0.3 

-3.H58 

-0.0231 

— 

— 

V1 

0-3 

1.7082 

0.0231 

0.277 

— 

— 

U2 

-3.1158 

0.0231 

6.5548 

0 

-3.U58 

-0.0231 

V2 

0.0231 

-0.277 

0 

2 . 4464 

0.0231 

-0.277 

U3 

— 

— 

-3.H58 

0.0231 

7.2249 

0.3 

v3 

— 

— 

-0.231 

-0.277 

0-3 

2.8162 

u4 

---- 

— 

— 

— 

-3.786 

-0.277 

V4 

— 

— 

— 

— 

-0.323 

-0.646 

U5 

0.5072 

-0.231 

-0 . 6688 

0.3 

— 

— 

V5 

0.0231 

-1.0621 

-0.3 

-0.3691 

— 

— 

u6 

0.6688 

-0.3 

1.0144 

0 

-06688 

0.3 

v6 

-0.3 

-0 . 3691 

0 

-1.1542 

0.3 

-0.3691 

UT 

— 

— 

-0 . 6688 

-0.3 

0.3457 

-0.2999 

Y7 

— 

— 

-0.3 

-0.3691 

-0.3001 

-1.5241 

U8 

----- 

— 

— 

— 

— 

— 

V8 

— — 

— 

— 

— 

— 

— 

U9 

---- 

— 

— 

— 

— 

— 

V9 

— 

— 

— 

---- 

— 

— 

u10 

— 

— 

— 

— 

— 

V10 

.... 

“  “  “  “ 

“  “  “  “ 
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TABLE  13  STIFFNESS  COEFFICIENT  MATRIX  FOR  EXAMPLE  PROBLEM 
(Cont'cL)  SOLUTION  INPUT 

NODES  4,  5  &  6 


STIFFNESS 

COEFFICIENTS 

x  106 

DISPLACE¬ 

MENT 

U4 

v4 

U5 

A 

u6 

v6 

U1 

— 

— 

0.5072 

0.0231 

0.6688 

0.3 

V, 

- - 

— 

0.0231 

1.0621 

0.3 

0 

U2 

---- 

— 

0.6688 

0.3 

l.oi44 

-1.1542 

VP 

— 

---- 

0.3 

-0 . 3691 

0 

0.3 

4 

-3.786 

-0.323 

— 

— 

-0.6688 

-0.3691 

V3 

-0.277 

-0. 646 

---- 

---- 

0.3 

— 

3.960 

-0 . 348 

— 

— 

— 

— 

\ 

u 

-0-348 

1.342 

— 

— 

---- 

— 

_ 

_  _  _  _ 

4.610 

0 

4.2916 

0 

5 

v 

0 

3.41 64 

0 

-0.544 

5 

u^ 

_____ 

_  _  _  _ 

-4.2916 

0 

9.8997 

0.3 

O 

- - 

— 

0 

-0.554 

0.3 

5 . 2626 

O 

U7 

V7 

U8 

V8 

-0.174 

0.671 

— 

— 

-4.9618 

-0.3001 

-0.625 

-0.696 

— 

— 

-0.2999 

-9.23 

— 

-0.5072 

-0.0231 

-0.6688 

-0.3 

— 

— 

0.0231 

-1.0621 

0.3 

-0.3691 

*9 

_  _  _  — 

—  -  —  - 

0. 6688 

-0.3 

0.3457 

-0 . 2999 

v_ 

_  _  _  _ 

-0.3 

-0.3691 

-0.3001 

-1.5241 

9 

TT 

---- 

----- 

— 

— 

—  —  — 

10 

V 

..MM  — 

_ 

— 

— 

— 

—  —  ~  — 

'10 
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TABLE  13  STIFFNESS  COEFFICIENT  MATRIX  FOR  EXAMPLE  PROBLEM 
(Cont'd)  SOLUTION  INPUT 


NODES  7,  8  &  9 


STIFFNESS 

COEFFICIENTS 

x  106 

DISPLACE¬ 

MENT 

D7 

V7 

u8 

V8 

U9 

VQ 

.7 

U1 

— 

— 

— 

— 

— 

- -  - 

V1 

— 

— 

— 

— 

-  - - 

_ 

U2 

-0.6688 

-0.3 

— 

— 

— 

_ _  _ 

V2 

-0.3 

-0.3691 

— 

— 

— 

_  _  _  _ 

U3 

0 . 3^57 

-0.3001 

— 

— 

- - 

_  _ _ 

V3 

-0.2999 

-1.5241 

— 

-  - _ 

_____ 

\ 

-0.174 

0.625 

— 

— 

-  -  -  - 

_  _  _  _ 

0.671 

-O.696 

— 

— 

----- 

- - 

U5 

— 

---- 

0.5072 

0.231 

-0. 6688 

-0.3 

V5 

— 

— 

0.0231 

-1.0621 

-0.3 

-0.3691 

u6 

-4.9618 

-0.2999 

-0. 6688 

0.3 

0.3457 

-0.3001 

V6 

-0.3001 

-0.923 

-0.3 

-0.3691 

-0.2999 

-1.5241 

U7 

5.6329 

-0.396 

— 

— 

-0.174 

0.625 

V7 

-0.396 

4.2082 

— 

— 

0.671 

-0. 696 

U8 

— 

— 

5-2849 

0-3 

4.9618 

-0.2999 

V8 

— 

— 

0.3 

3.7862 

-0.3001 

-0.923 

U9 

-0.174 

O.67I 

-4.96l8 

-0.3001 

5.6329 

-0.396 

V9 

0.625 

-O.696 

-O.2999 

-O.923 

-0.396 

4.2082 

U10 

— 

— 

-0.1615 

-O.323 

-0.174 

0.671 

V10 

-O.277 

-1.432 

0.625 

-0.696 
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TABLE  13  STIFFNESS  COEFFICIENT  MATRIX  FOR  EXAMPLE  PROBLEM 

SOLUTION  INPUT 

NODE  10 


DISPLACE¬ 

MENT 


U„ 


U2 

V. 

u, 

V, 

U|. 


U, 

v' 


u. 


u7 

v7 

ue 

Uc 

1' 


u. 


10 


10 


STIFFNESS  COEFFICIENTS  x  10 6 


U10 

V10 

LOADS  x  10 

- _ 

— 

-0.19245 

---- 

— 

-0.12405 

0 

-0.19245 

0 

-0.19245 

— 

— 

_  _  _ 

— 

— 

0.2379 

— 

— 

-0.0909 

— 

— 

-0.32925 

0 

0 

0 

0.32925 

— 

— 

— 

— 

0.19245 

-0.1615 

0.277 

-0.32925 

-0.323 

-1.432 

0 

-0.174 

0.625 

0.32925 

0.671 

-0.696 

0.19245 

0.3555 

-0. 348 

0.04545 

-0.348 

2.128 

0.21495 

WADD  TR-411  PART  I 


76 


TABLE  14  NODE  DISPLACEMENTS  IN  EXAMPLE  PROBLEM  WING 


Assumption 


=  U8  = 


u. 


10 


=  0 


Displacement 


Displacement 

Magnitude 

Mode 

in.  x  10 

U2 

0. 042428 

v2 

-0.011423 

U3 

O.O85230 

-0.022482 

0.130206 

-0.040740 

V5 

0.079823 

u6 

0.051957 

v6 

0.075393 

U7 

0.103587 

V 

0.068127 

*8 

0.162927 

u9 

0.051068 

v9 

0.161250 

vio 

0.24839 

Symbols 


U  represents  horizontal  displacement 
V  represents  vertical  displacement 
Subscript  denotes  node  point  location 
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TABLE  15  STIFFNESS  COEFFICIENTS,  FORCES  AND  DISPLACEMENTS  FOR  RECTANGULAR 

PANELS  IN  EXAMPLE  PROBLEM  WING 


Q 


ch 

Ch 

<D 

O 

O 

W 

m 

<u 

Pl 

'H 

Ch 

-H 

-P 

CO 


o 


m 


< 


00 

I 


rH 

PO 

CM 

O 


CO 

Lf\ 

t~- 

I — 1 


1 — 1 

1 — 1 

1 — i 

rH 

o\ 

on 

CM 

on 

vo 

CM 

ON 

OJ 

on 

0 

O 

0 

1 — 1 

rH 

rH 

on 

CM 

,qn 

CM 

ON 

VO 

0 

O 

O 

on 

PO 

• 

* 

• 

• 

CM 


OO 

CM 

O 


1 — I 

I 


PO 


on 


-4" 

m 

PO 


PO 

1 


CO 

co 


t- 

c 

CM 


CO 

CM 

1 - 1 

co 

rH 

$ 

f- 

O 

on 

OJ 

LON 

b~ 

on 

OJ 

VO 

on 

ITN 

0 

1 — 1 

0 

CM 

CO 

PO 

t- 


PO 


OJ 

rH 

1 — 1 

■ — 1 

CO 

on 

b- 

on 

CM 

on 

OJ 

OJ 

ON 

< 

on 

b- 

0 

CVl 

0 

O 

• 

• 

• 

• 

• 

• 

PO 


-4" 

PO 


PO 

CM 

O 


OJ 

rH 

CO 

co 

b- 

0 

on 

OJ 

$ 

b— 

on 

LA 

b- 

LTN 

0 

VO 

on 

on 

on 

r — 1 

• 

1 — 1 

H 

OJ 

rH 

H 

CO 

on 

PO 

on 

VO 

CM 

b- 

on 

PO 

O 

CM 

t- 

O 

ICO 


PO 

I 


H* 

00 

H 

CM 

rH 

co 

b- 

LA 

on 

t— 

on 

oQ 

on 

b- 

OJ 

O 

OJ 

NO 

on 

on 

H 

0 

LTN 

0 

VO 

on 

• 

• 

• 

• 

• 

• 

• 

• 

H 

rH 

1 

1 

1 

1 

pq 

pq 

O 

O 

_P 

.  0 

> 

> 

to 

> 

to 

> 
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.3691  -  .0231  -  .0921  .0231  -  .277  .3  *7382 


TABLE  15  STIFFNESS  COEFFICIENTS,  FORCES  AND  DISPLACEMENTS 
(Cont'cL)  FOR  RECTANGULAR  PANELS  IN  EXAMPLE  PROBLEM  WING 


Forces 


PANEL 

UA 

VA  UB 

VB 

1256 

-0.10515 

■0.04662  +0.10873 

-0.05837 

2367 

-0.10800 

■0.05187  0.10941 

-0.05610 

5689 

-O.II765 

0.05060  0.11977 

-O.O5585 

PANEL 

UC 

VC  UD 

VD 

1256 

-0.11510 

0.05379  0.11152 

0.04920 

2367 

-0.11532 

0.05^:83  0.11706 

0.05162 

Displacements 

PANEL 

PANEL 

PANEL 

1256 

2367 

5689 

UA 

0 

0.042428 

0 

0 

-0.011423 

0.079823 

UB 

0.042428 

0.08523 

0.51957 

VB 

-0.011423 

-0.022482 

0.075393 

^C_ 

0 

0.051957 

0 

vc 

0.079823 

0.075393 

0.162927 

0.051957 

0.103587 

0.051068 

VD 

0.075393 

0.068127 

0.16125 
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TABLE  1 6  STIFFNESS  COEFFICIENTS,  FORCES  AND  DISPLACEMENTS 

FOR  TRIANGULAR  PANET.fi  IN  EXAMPLE  PROBLEM  WING 

Stiffness  Coefficient  Matrix 


UA 

VA 

UB 

VB 

la. 

la. 

UA 

2.0075 

.6 

-1.846 

-  .323 

-  .1615 

-  .277 

VA 

.  6 

1.108 

-  .277 

-  .646 

-  .323 

-  .462 

UB 

-1.846 

t— 

OJ 

1 

1.846 

0 

0 

.277 

VB 

-  -323 

-  .646 

0 

.646 

.323 

0 

uc 

-  .1615 

-  .323 

0 

•  323 

.1615 

0 

Zc_ 

-  .277 

1 

. 

ON 

ro 

.277 

0 

0 

.462 

Forces 


PANEL 

UA 

VA 

UB 

VB 

347 

-0.10520 

-0.04846 

0.10813 

-0.00587 

-0.00293 

0.05432 

679 

-0 . 11660 

-0.04899 

0 . 11909 

-0.00498 

-0.00249 

0.05397 

8-9-10 

-0.11740 

-0.05255 

0.11794 

-0.00108 

-0.00054 

0.05363 
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TABLE  1 6  STIFFNESS  COEFFICIENTS,  FORCES  AND  DISPLACEMENTS  FOR  TRIANGULAR 

(Cont'd)  PANELS  IN  EXAMPLE  WING 


Displacements 


PANEL 

PANEL 

PANEL 

3^7 

679 

8-9-10 

UA 

0.085230 

0.051927 

0 

VA 

0.022482 

0-75393 

0.162927 

UB 

0.130206 

0.103587 

0.51068 

VB 

0.040740 

0.069127 

0.161250 

uc 

0.103587 

0.051068 

0 

0.068127 

0 . 161250 

0.2^839 
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TABLE  17  STRESSES  IN  RECTANGULAR  AND  TRIANGULAR  PANELS  IN  EXAMPLE  PROBLEM  WING 


Panel  1256 


Aax  =  - 


Aoy  = 


U  U 

l  5  =  55061 


V  V 
1+2 


51495 


xy 


o  = 

X 


Vn  Vr  V  V. 

- l_t - - La - =  _  35g2 

-68400  +  Aax  =  -13339  lbs. /in. 2 


a  = 

y 


=  -68400  +  Aay  =  -I69O5  lbs. /in. 


Aax  = 


Aay  = 


Panel  2367 

U.  ,  U  ? 

=  55828  lbs. /in. 

V .  V 

6+7  


=  53988  lbs. /in/ 


rxy  = 


U  U  U 
6  +  7  -  2  -  3  


-14145  lbs. /in. ‘ 


ax 


68400  +  Aax  -  -  12572  lbs. /in.' 


°y 


68400  +  Aay  =  -  14412  lbs. /in. ‘ 


Aax  = 


Panel  5689 

U,  U0  Q 

■4  /  -2 —  =  59206  lbs. /in. 


Aay  =  ■+  2.  =  53226  lbs. /in. 2 

U  \J,  UQ  U  0 

Txy  - - ^  -  6g~ — — — ^  =  -211 6  lbs. /in. 2 


ax  =  -  68400  +  /Lax  =  -9194  lbs. /in/ 


oy 


68400  +  Aay  =  -15174  lbs. /in/ 


"82" 
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TABLE  17  STRESSES  IN  RECTANGULAR  AND  TRIANGULA  PANELS  IN  EXAMPLE  PROBLEM  WING 
(Cont 'd) 

Panel  347 

Uk  P 

Aax  =  —  =  5^063  lbs. /in. 

V7  P 

Any  =  =  54320  lbs. /in. 

\  ? 

Txy  =  —  =  -  2933  lbs. /in. 

ox  -  -68400  +  A ax  =  -1^337  lbs. /in. 2 
oy  =  -68400  +  Aay  =  -14080  lbs. /in. 2 

Panel  679 

U7  ,  0 

Aax  =  2  =  59546  lbs. /in. 

Aay  =  ^  =  53970  lbs. /in. 

V7  ;  ? 

Txy  =  2  =  “2490  lbs. /in. 

ax  =  -68400  +  Aax  =  -8854  lbs. /in. 2 

ay  =  -68400  +  Aay  =  -14433  lbs. /in. 2 

Panel  8-9-10 

UQ  p 

Aax  =  =  58971  lbs. /in. 

V 

Aay  =  =  53630  lbs. /in. 2 

v9  ,  2 
Txy  =  =  -541  lbs. /in. 

ax  =  -68400  +  Aax  =  -9429  lbs. /in. 2 

ay  =  -68400  +  Aay  =  -14771  lbs. /in. 2 
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TABLE  18  STRAIN  DERIVATION  AND  STRESS  INCREMENTS  IN  EXAMPLE  PROBLEM  WING 
SPARS  EXCLUSIVE  OF  FRONT  SPAR 


STRAIN  DERIVATION  STRESS-CAPS  STRESS -WEBS 

E  =  .96  x  10 '  E  =  .98  x  10 ' 
A 0  (lbs./in.^)  A 0  (lbs./in.^) 


"X"-direction  Components 


ex. 


1-2 


ex 


2-3 


ex 


3-4 


ex 


5-6 


ex 


6-7 


ex, 


8-9 


U2  - 

UI 

10 

U3  - 

U2 

10 

u4  - 

U3 

10 

U6  - 

U5 

10 

U7  - 

U6 

10 

1 

ON 

1 

U8 

10 

0.0042427 
0. 0042801 
0.0044977 
0.0051957 
0.0051630 
0.0051068 


40730 

41089 

43178 

49879 

49565 

49025 


41579 

41945 

44078 

50918 

50597 

50046 


"Y"-direction  Components 


£yl-5  = 

1  0 

OJ 

LfN 

1 

V. 

V 

6  - 

? 

6y2-6  = 

20 

V^ 

V 

ey3-7  = 

-7  -. 

20 

vQ 

V 

8  - 

-1 

€y5-8  = 

20 

V 

€y6-9  = 

JLz. 

20 

v 

6 

10 

_ 

ey8-l0  = 

20 

0.0039912 
0.0043407 
0.0045305 
0. 0041553 
0.0042929 
=  0.00427 


38316 

39114 

41821 

42539 

43492 

44398 

39891 

40722 

41211 

42070 

41023 

41877 

WADD  TR  60-4ll  PART  I 


84 


TABLE  18  STRAIN  DERIVATION  AND  STRESS  INCREMENTS  IN  EXAMPLE  PROBLEM  WING 
(Cont'd)  SPARS  EXCLUSIVE  OF  FRONT  SPAR 


STRAIN  DERIVATION 


STRESS 

E  =  .88  x  10< 
Aa  (IBs. /in.  ) 


_^L_ 

el+-7  500 

_  _l_ 

e7-9  ~  500 

=  _1_ 
6 9-10  “  500 


10  (u4  -  u?) 
10  (u7  -  u9) 
1°  (u9  -  u10) 


20  (V^  -  V?) 

20  (V?  -  V9) 

“  20  (y9  -  V 


0.0048870 
0.0047753 
=  0.0045069 


43006 

42023 

39661 
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TABLE  19  STRESSES  IN  EXAMPLE  PROBLEM  WING  SPARS 


MEMBER 

THERMAL  STRESS 
c  =EoT 
(Locked) 

STRESS 

INCREMENT 

Ac 

TOTAL  STRESS 

a=a  +Ao 

0 

WING  SPARS  EXCLUSIVE  OF  FRONT  SPAR  (l) 

1-2 

-36900 

40731 

3831 

2-3 

-36900 

41089 

4189 

3-4 

-36900 

43179 

6279 

5-6 

-36900 

49878 

12978 

6-7 

-369OO 

49565 

12665 

8-9 

-369OO 

49025 

12125 

1-5 

-369OO 

38316 

l4i6 

5-8 

-369OO 

41670 

4770 

8-10 

-369OO 

43492 

6592 

2-6 

-369OO 

39891 

2991 

6-9 

-369OO 

41211 

4311 

3-7 

-369OO 

41022 

4122 

FRONT 

SPAR  (2) 

4-7 

-50700 

43006 

-7694 

7-9 

-50700 

42022 

-8677 

9-10 

-50700 

39662 

-IIO38 

NOTES : 

(1)  T  =  200°F,  E  = 

:  .96  x  10^  lbs/in2,  a 

=  12.8  x  lo"6 

(2)  T  =  450°F,  E  = 

=  .88  x  10^  lbs/in2,  a 

=  12.8  x  10“6 
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TABLE  20 


STRESSES  IN  EXAMPLE  PROBLEM  WING  SPAR  WEBS 


MEMBER 

THERMAL  STRESS 

a  =EoT  (1) 

0 

(Locked) 

STRESS 

INCREMENT 

A0 

TOTAL  STRESS 

a=o  +Aa 

0 

1-2 

-18750 

41578 

22828 

2-3 

-18750 

4194  5 

23194 

3-4 

-18750 

44077 

25327 

5-6 

-18750 

50917 

32167 

6-7 

-I875O 

50598 

31848 

8-9 

-18750 

50046 

31296 

1-5 

-I875O 

39114 

20364 

5-8 

-18750 

42538 

23788 

8-10 

-18750 

44398 

25648 

2-6 

-I875O 

40722 

21972 

6-9 

-18750 

42070 

23320 

3-7 

-I8750 

41877 

23127 

NOTE: 

(1)  T  =  100 °F., 

E  =  .98  x  107, 

a  =  12.8  x  lO-6 
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TABLE  21 


PRINCIPAL  STRESS  COMPUTATIONS  FOR  EXAMPLE  PROBLEM 
WING  PANELS  AT  TIME  ZERO 


PANEL 

STRESS  COMPONENTS 

PRINCIPAL 

STRESSES 

0 

X 

a 

X 

T 

xy 

cri(l) 

°2(2) 

1-2- 5-6 

-  13339 

-  16905 

-  3582 

-  11124 

-  19122 

2-3-6-Y 

-  12571 

-  14412 

-  i4i4 

-  11805 

-  15180 

5-6-8-9 

-  9193 

-  15174 

-  2116 

-  8518 

-  15847 

3-4-7 

-  14337 

-  14085 

-  2932 

-  11274 

-  17145 

6-7-9 

-  7639 

-  14433 

-  2490 

-  7962 

-  15327 

8-9-10 

-  9^29 

-  14770 

-  541 

9375 

-  14823 

NOTES : 


(1) 


fa  +o 

[  _ y 

v  2 


\ 

; 


+ 


0 


1 


(2) 


2 

T 

xy 


a 


2 


WADD  TR  60-4ll  PART  I 


88 


TABLE  22 


ANGLE  OF  PRINCIPAL  AXES  COMPUTATIONS  FOR  EXAMPLE 
PROBLEM  WING  AT  TIME  ZERO 


PANEL 

STRESS  COMPONENTS 

a 

X 

a 

y 

T 

xy 

1-2-5-6 

-  13339 

-  16905 

-  3582 

2-3-6-7 

-  12571 

-  14412 

-  l4l4 

5-6-8-9 

-  9193 

-  15174 

-  2116 

3-4-7 

-  14337 

-  14085 

-  2932 

6-7-9 

-  7639 

-  14433 

-  2490 

8-9-10 

-  9429 

-  14770 

54i 

PANEL 

ANGLE 

OF  PRINCIPAL  AXES 

DATA 

Tan  2a 
(1) 

2a 

a 

Sin2a 

Cos2a 

Sin2a  Cos2a 

1-2- 5-6 

-  2.008 

63  °31 ' 

31°45' 

0.2769 

0.7231 

0.6677 

2-3-6-7 

1-537 

56°57' 

0 

CO 

VO 

0.2271 

0.7728 

0.4l89 

5-6-8-9 

-  0.7078 

35°17’ 

17°38’ 

0.0917 

0.9082 

0.2887 

3-4-7 

22.866 

87°29' 

43°45> 

0.4782 

O.5217 

0.4995 

6-7-9 

-  0.8928 

4i°45  ' 

20°53 ' 

0.1272 

0.8729 

0-3331 

8-9-10 

-  0.2028 

11°28' 

5°44' 

0.00998 

O.990 

0.0994 

NOTE: 


(1) 


Tan  2a  = 


2t 

SZ 

a  -  a 

x  y 
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TABLE  23 


EQUIVALENT  UNIAXIAL  STRESSES  IN  EXAMPLE  PROBLEM  WING  PANELS 
ACCORDING  TO  von  MISES-PRAGER  CRITERIA 


PANEL 

PRINCIPAL 

°1 

STRESSES 

°2 

EQUIVALENT  STRESS 

0 

0 

(1) 

TEMP. 

°F 

STRAIN 
RATE,  £ 

(2)  ° 

POWER 

DISSIPATION 

a  .  e  (3) 

0  0 

1-2-5-6 

-11124 

-19122 

16572 

450 

-4 

1.28x10 

2.12 

2-3-6- 7 

-11805 

-15180 

13806 

450 

-4 

0 . 71x10 

O.98O 

5-6-8-9 

-  8518 

-15847 

13746 

450 

0 . 70xl0_1+ 

O.962 

3-4-7 

-11274 

-17145 

15103 

450 

-4 

0.95x10 

1.43 

6-7-9 

-  7962 

-15327 

13276 

450 

-4 

0.13x10 

O.836 

8-9-10 

-  9375 

-14823 

12985 

450 

0.58x10"^ 

0.753 

NOTES: 


(2)  Taken  from  calibration  chart. 

(3)  Expressed  in  pound  minutes  per  square  inch. 
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TABLE  24  COMPUTATIONS  INVOLVING  DEVIATOR  COMPONENTS  IN  EXAMPLE  PROBLEM 


WING  PANELS 


PRINCIPAL 

STRESSES 


HYDROSTATIC 


DEVIATOR 


PANEL 

°l 

a2 

(1) 

si 

S2 

1-2 -5-6 

-  11124 

-  19120 

-  10081 

-  1042 

-  9040 

10082 

2-3-6- 7 

-  11805 

-  15180 

-  8995 

-  2809 

-  6i84 

8993 

5-6-8-9 

-  8518 

-  15847 

-  8122 

-  396 

-  7725 

8121 

3-4-7 

-  11274 

-  17145 

-  9472 

-  1801 

-  7672 

9473 

6-7-9 

-  7962 

-  15327 

-  7762 

-  200 

-  7564 

7764 

8-9-10 

-  9375 

-  14823 

-  8060 

-  1310 

-  6757 

8067 

PANEL 

ENERGY 

DISSIPATION 

COMPONENT 

POWER 

DISSIPATION 

(3) 

PROPORTION¬ 

ALITY 

FACTOR 

(4) 

1-2- 5-6 

184.46  x  108 

2.12 

O 

1.15  x  10 

2-3-6-7 

£ 

127.06  X  10 

O.98O 

0.771  x  10“8 

5_6-8-9 

125.83  x  106 

0.9262 

O.765  x  10"8 

3-4-7 

151.84  x  106 

1.43 

O 

0.942  x  10 

6-7-9 

r 

117.52  x  10 

O.836 

O 

0.712  X  10 

8-9-10 

r 

111.47  X  10 

O.753 

O.676  x  10'8 

NOTES : 


(l)  Hydrostatic  term  =  (a1  +  a2)/3 


(2)  S1  =  ax 


(3)  See  Table  XXIII 


(^)  Energy  Dissipation  Component 
Power  Dissipation 


(^) 
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TABLE  25  CREEP  RATES  IN  PRINCIPAL  DIRECTIONS  AND  IN  THE  ELASTIC  STRAIN 
INCREMENT  IN  "LOCKED"  EXAMPLE  PROBLEM  WING  PANELS 


PROPORTIONALITY  FACTOR 

PANEL 

DEVIATOR 

S1 

COMPONENTS 

S2  S3 

PROPORTIONALITY  FACTOR 

1/r 

1-2-5-6 

-  10*4-2 

-  9040 

10082 

1.15  x  10"^ 

2-3-6-j 

-  2809 

-  6184 

8993 

-8 

0.771  x  10 

5-6-8-9 

-  396 

-  7725 

8121 

,  -8 

0.765  x  10 

3-4-7 

-  1801 

-  7672 

9^73 

,  -8 

0.942  x  10 

6-7-9 

-  200 

-  7564 

7764 

-8 

0.712  x  10 

8-9-10 

-  1310 

-  6757 

8067 

.  „  -8 

O.676  x  10 

PANEL  CREEP  RATES  IN  PRINCIPAL 

DIRECTIONS  (I) 


€1 

e2 

_ 

1-2-5-6 

-  1.198x10" 5 

-  10. 3 96x10” 5 

11.594x10 

2-3-6- 7 

-  2.166x10“ 5 

-  4 . 768x10" 5 

6. 934x10" 5 

5-6-8-9 

-  0.303x10" 5 

-  5. 909x10" 5 

6. 212x10” 5 

3-4-7 

-  1.697x10” 5 

-  7. 227x10” 5 

8. 924x10" 5 

6-7-9 

-  0.l42xl0"5 

5. 385x10" 5 

5. 527x10" 5 

8-9-10 

-  0.885x10” 5 

-  4 . 568x10' 5 

5. 4 53x10" 5 
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TABLE  25  CREEP  RATES  IN  PRINCIPAL  DIRECTIONS  AND  IN  THE  ELASTIC  STRAIN 
(Coat'd)  INCREMENT  IN  "LOCKED"  EXAMPLE  PROBLEM  WING  PANELS 


PANEL 

ELASTIC  STRAIN 
INCREMENTS  (2) 

*1 

^2 

^3 

1-2-5-6 

-  4. 792x10“ 5 

41.584x10“ 5 

46. 376x10“ 5 

2-3-6- 7 

-  8.664x10  ^ 

-19. 072x10“ 5 

27. 73 6x10“ 5 

5-6-8-9 

-  1.212x10“ 5 

-23. 63 6x10“ 5 

24. 848x10“ 5 

3-4-7 

-  6. 788x10“ 5 

-28. 908x10“ 5 

35. 696x10“ 5 

6-7-9 

-  0.568x10“ 5 

-21.540x10“ 5 

22.108x10“ 5 

8-9-10 

-  3.540x10“ 5 

-18. 272x10“ 5 

21.812x10“ 5 

NOTES: 


(1) 

Vr 4  > 

S2/r  ■  e2  ; 

S3/r  ‘  4 

(2) 

^At  - 

e2At  =  Ae2j 

e^At  =  Ae 

where  At  = 

t ime  increment 

-  4  min. 
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TABLE  26  RATE  OF  CHANGE  CF  STRESS  WITH  RESPECT  TO  TIME  IN  PRINCIPAL  DIRECTIONS 

OF  EXAMPLE  PROBLEM  WING 


PANEL 

CROSS 

SECTIONAL 

AREA 

SQ..  IN. 

STRAIN  RATES 

(1) 

INCREMENTAL  STRESS- 
INCREMENTAL  TIME 
RATIO  (2) 

1 — 1  1 

3  k 

^2 

At 

Ao3 

At 

41 

e2 

'3 

1-2-5-6 

2 

1.198 

10.39 

n.594- 

44.44 

607.3 

651.0 

2-3-6-j 

2 

2.166 

4.76 

6.934 

122.7 

152.6 

274.5 

3-4-7 

1 

1.697 

7.22 

8.924 

84.76 

351.2 

435.3 

1-5 

2 

0 

0 

0 

2-6 

2 

0 

0 

0 

3-7 

2 

0 

0 

0 

4-7 

2 

0 

0 

0 

NOTE:  (1)  Multiply  values  shown  hy  10  J  to  arrive  at  true  value. 


(2) 


Ao^ 

At 


E 


zAi 

ZA. 
1 


lbs  I sq  in | min 


where  i  =  1,  2,  3  • 

c 

E  =  8.08  x  10 
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TABLE  27  PRINCIPAL  STRESS  INCREMENTS  ARISING  FROM  RELAXATION  IN  LOCKED  EXAMPLE 
PROBLEM  WING  PANELS  AFTER  A  CREEP  INTERVAL  OF  4  MINUTES 


PANEL 

ELASTIC  STRAIN 

INCREMENTS 

“l  ^2 

PRINCIPAL  STRESS 
INCREMENT  DUE 

TO  RELAXATION  (l) 

Aa1  Aa2 

I-2-5-6 

4. 792x10 

41.584x10”'’ 

+  153^ 

+  3820 

2-3-6-7 

-8. 662x10“ 5 

-19. 072x10“ 5 

+  1277 

+  1924 

5-6-8-9 

-1.212x10“ 5 

-23. 636x10” 5 

+  737.2 

+  2131 

3-^-7 

-6. 788x10“ 5 

-28.905x10“5 

+  1373 

+  2748 

6-7-9 

-0.568x10“ 5 

-21.540xlo“5 

+  642.2 

+  1928 

8-9-10 

-3.540x10" 5 

-18. 272x10“ 5 

+  901.1 

+  1717 

NOTE: 

(1)  A dx  =  -  — — g  (Ae1  -  VAe  ),  and 

1-P 

E 

A a 2  = - ^  -  VAe^)  where 

1-M- 

- 2  =  8.879  x  10  lb/sq  in.,  and 

1-H 

E  =  8.08  x  106  and 


P  =  0.3 
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TABLE  28  STRESS  INCREMENTS  DUE  TO  RELAXATION  IN  LOCKED 
PANELS  OF  EXAMPLE  PROBLEM  WING 


PANEL 

PRINCIPAL 

ANGLE 

STRESS 

DATA 

PRINCIPAL  STRESS 
COMPONENTS 
(INCREMENTAL) 

Ao1  A ag 

PRINCIPAL  STRESSES 
(INCREMENTAL) 

Sin  a 

Cos  a 

Aa 

X. 

(1)' 

Ac 

y 

(2) 

1-2- 5-6 

O.2769 

0.7231 

1534. 

3820 

2167 

3186.8 

2-3-6-7 

0.2271 

0 . 7728 

1277. 

1924 

1423.8 

1777- 

5-6-8-9 

0.0917 

0.9082 

737-2 

2131 

864.9 

2002.6 

3-4-7 

0.4782 

0.5217 

1373 

2748 

2030.3 

2090.6 

6-7-9 

0.1271 

O.8729 

624.2 

1928 

789.9 

1762.3 

8-9-10 

0 . 00998 

O.99OO 

801.1 

1717 

810.2 

1708 

NOTE : 

(1)  Ao  =  An,  •  Cos  a  •  Ac^  Sin  a. 

x  1  2 

(2)  Aa  =  A a_  Cos  a.  •  Ac-,  Sin  a 

y  2  1 
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TABLE  29 


SHEAR  STRESS  INCREMENTS  DUE  TO  RELAXATION  IN  LOCKED  PANELS  OF 
EXAMPLE  PROBLEM  WING 


PANEL 

PRINCIPAL 

STRESS 

(INCREMENTAL) 

A"  U- T~ 

txy 

PRINCIPAL  STRESS 
ANGLE  DATA 

Sin  a  .  Cos  a 

PRINCIPAL  STRESS 
COMPONENTS 
(INCREMENTAL) 

Aax  Act2 

1-2-5-6 

456.7 

0.4477 

2167 

3187 

2-3-6-J 

147-9 

0.4189 

1424 

1777 

5-6-8-9 

328.5 

0.2887 

865. 

2003 

3-4-7 

30.47 

0.4995 

2030 

2091 

6-7-9 

323.1 

0.3331 

790 

1762 

8-9-10 

89.26 

0.0994 

810 

1708 

NOTE: 


(l)  A  =  (Aa2  -  Act1)  ‘  Sin  a  Cos  a 
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TABLE  30  FORCE  INPUT  AT  LOCKED  NODES  OF  RECTANGULAR  PLATES  IN  EXAMPLE 


PROBLEM  WING 


PANEL 

FORCE 

COMPONENT 

AF 

-^(1) 

AS 

-^(2) 

AF 

2y(3) 

AS 

/(M 

1-2-5-6 

2167 

457 

6347 

913 

2-3-6-T 

1424 

148 

355^ 

296 

5-6-8-9 

865 

329 

4005 

657 

NOTES: 


(1)  AFx  Aa 

— —  =  —7=—  *Ax  since  when  A 

2  2  x 


=  10  x  0.2  =  2  sq.  in. 


AF  =  ^  .  2  =  Aa 
x  2  > 


(2)  AS 


2Aa 


(3)  AFir  Aa 


„  Ay  since  when  A 

2  y 


=  20  X  0.2  =  4  sq.  in. 


AF 

~2 


*  =  Aa 


(4)  s 
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TABLE  30  FORCE  INPUT  AT  LOCKED  NODES  OF  RECTANGULAR  PLATES  IN  EXAMPLE 
(Cont.)  PROBLEM  WING 


PANEL 

DEGREE  OF  FREEDOM^) 

Ul 

U2 

U3 

U4 

Vi 

V2 

V3 

V4 

1-2-5-6  - 

2624 

1710 

1710 

2624 

7287 

5466 

54 66 

7287 

2-3-6- 7  - 

1572 

1276 

1276 

1572 

3350 

3258 

3258 

3850 

5-6-8-9  - 

1194 

536 

536 

1194 

4662 

3676 

3676 

4662 

NOTE 


(5) 


The  degrees  of  freedom  are  computed  from  the  following  relations 


X 

%  w 

11 

AS 

X 

2 

V1 

AF 

=  _j:  _ 

2 

AS 

u 

2 

AF 

U  =  — - 

2  2 

AS 

X 

2 

V2 

AF 

=  _ 

2 

AS 

—I 

2 

AF 

V-v5 

1 

10  fe 

X 

V3 

AF 

_ y  _ 

2 

AS 

2 

X 

%  W 

II 

hr1" 

AS 

_  _£ 

2 

AF 

y  + 

2 

AS 

-JL 

2 

applicable  at  each  node  point  represented  at  the  corner  of  a  given 
panel. 
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TABLE  31 


FORCE  INPUT  AT  LOCKED  NODES  OF  TRIANGULAR  PLATES  IN  EXAMPLE 
PROBLEM  WING 


PANEL 

AF 

X 

2 

FORCE 

AS 

~ 2~(l) 

COMPONENT 

AF 

_ 1 

2 

/IS 

-^(2) 

3-4-7 

2030 

30 

4i8i 

6i 

6-7-9 

790 

323 

3525 

646 

8-9-10 

8io 

8  9 

34l6 

179 

NOTES J 


(1) 


AS 

_ x 

2 


(2) 


AS 

2 


AS 

2 


A  since  when  A  =  10  x  0.2  =  .2  sq.  in. 
x  x 


2  =  T 


xy 


At 


2 

At 


xy  .  A  since  when  A  =  20  x  0.2  =  .4  sq. 

—  y  x 


xy  .  4  =  2t 


xy 
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in. 


TABLE  31 
(Cont. ) 


FORCE  INPUT  AT  LOCKED  NODES  OF  TRIANGULAR  PLATES  IN  EXAMPLE 
PROBLEM  WING 


NOTE} 

(3) 


PANEL 

DEGREE  OF  FREEDOM  (3) 

U3 

X 

v3 

V7 

3-4-7 

-2060 

2060 

-L242 

4242 

6-7-9 

-1113 

1113 

-4171 

4171 

8-9-10 

-  899 

899 

-3595 

3595 

The  degrees  of  freedom  are  computed  from  the  following  relations 


u 


7 


=  o 


u 


3 


AF 

x 

2 


AS 

_ x 

2 


V 


3 


AS 

_£ 

2 


u4 


AF  AS 

_ x  +  _ x 

2  2 


AF  AS 

V  =  — ^  +  — Z 

7 


applicable  at  each  node  point  represented  at  the  corner  of  a 
given  panel. 
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TABLE  32  RESULTANT  FORCES  IN  "X"  DIRECTION  AT  LOCKED  NODES  (l)  OF  EXAMPLE  PROBLEM  WING 


i 


r _ , 

0 

0 

H 

1 — 1 

£3 

on 

o> 

£3 

•  CO 

ao 

£3 

•  £- 

£3 

CM 

, — , 

•  VO 

VO 

§ 

£3 

0 

£h 

P>4 

1 — 1 

O 

•  LTV 

LT\ 

O 

w 

0 

, - , 

■  -d" 

£3 

•  00 

or 

£3 

'  CM 

CM 

£3 

>•  H 

1 — i 

£3 

I 

I 


I  I  I 
I  I  I 
I  I  I 


o 


-d-  1 

1  ON 

OO 

ON  1 

|  • 

OV 

r— 1  1 

1  ON 

O 

1 — 1  1 

1  CO 

CM 

LTV 

1 

1  CSV 

OO 

VO  1 

1  .  ON 

-d" 

00  1 

1  00 

H 

Lf\  I 

1 

II  I  I  I  I 


CM  1 

1  OO 

1  LCN 

h-  l 

1  1 — 1 

1  CO 

LT\  | 

1  1 — 1 

1  VO 

1 — 1  1 

1  1 — 1 

1  CM 

-=t- 

VO 

VO 

CO 

1 — l 

1 

1  1 — 1 

CM 

£— 

00 

1 — 1 

1  D— 

VO 

CM 

LTV 

1 — 1 

1  I>- 

CM 

O 

1 — 1 

1 — 1 

1 

1 

1 

1 

o\ 

1 

1 

1 

1 

i  0 

£- 

1 

1 — 1 

1 

1  ON 

1 — 1 

1 

1 — 1 

1 

1  CM 

I  I  I 


O  1 

1  O 

VO  1 

1  VO 

O  1 

1  0 

CM  1 

1  CM 

0 

vo 

L 

1  vo  1 

1  O  1 

!  cS 

CM 

1  CM  1 

1  t>- 

1 — 1 

1  1 

1 

I 


O  CM  1  1 

1— 1  1  1 

LT\  I  I 

1 — I  1 — I  I  I 


QO 

CO 


-d*  I  l  I  I 

CM  I  1  1  1 

VO  l  1  1  1 

CM  l  1  l  I 

1  1  1  1 


-d- 

CM 

VO 

CM 


vo 

Lf\ 

c— 

1 

vo 

a\ 

CO 

h- 

on 

0 

1 — 1 

1 — 1 

1 

1 

1 

1 

1 

1 

a3 

CM 

00 

vo 

-d" 

£— 

o\ 

-P 

1 

1 

1 

1 

1 

1 

O 

r— i 

CM 

LT\ 

CO 

VO 

CO 

EH 

O 

-P 

O 

03 


aJ 

o 

iH 

Jh 

<1) 

& 

B 

u 

o 

ch 


-P 

I 

M 


s 

O 

-P 

O 

a) 

a 

tJ 

h 


£ 

o 

<u 

Q) 

P 

Eh 

o 


<u 

u 

bD 

& 


CM 


£3 
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TABLE  33  RESULTANT  FORCES  IN  "Y"  DIRECTION  AT  LOCKED  NODES  (l)  OF  EXAMPLE  PROBLEM  WING 


I 


o 

I — 1 


I 

I 

I 

I 


I  I  I  Lf\  LfN 

II  I  ON  ON 

I  I  I  LfN  LT\ 

i  i  i  oh  on 


ON 


I  I 

I  I 

I  I 

I  I 


CM 

1  H 

1  OH 

VO 

1  N- 

1  OH 

VO 

1  H 

1  CO 

-4 

1 

1  00 

CO 


I  I  IN-  I 

I  I  NO  I 

I  I  CO  I 


I  I 

I  I 

I  I 

I  I 


i — I 

CO 


O 

1  CM 

1  LTN 

CM 

LfN 

l  -4 

1  ON 

On 

CO 

l  CM 

1  LTN 

o 

OH 

l  -4 

1  OH 

CO 

I 


IN- 

CO 

CM 

IN¬ 


CH  NO 
LTN  N- 
CM  NO 
OH  OH 


:  & 

H  i  VO 

-4  I  CM 


VO 

MO 

-4 

LCN 

I 


I  CM  I  I 

I  VO  I  I 

I  MO  I  I 

i  I  r 


i 

i 

i 

i 


-4 

O 

CO 


l  i  l  l  l 

I  <1  I  l  o 

I  I  I  I  I 


& 

1  CM  1 

1  -4  1 

I  o 

l  o 

OJ 

1  CM  1 

1  LfN 

co 

'  -4  l 

1  IN- 

I  I  I 


MO 

VO 

-4" 

UA 

I 


O 

ir\ 

00 

co 


i  i  i 
i  t  i 
i  i  i 
i  i  I 


VO 

rH 

CO 

ON 


CO 

CM 

IH 

i 


l  C'- 

I  CO 

I  CM 

I  t>- 

i 


a 

o 

•H 

-P 

O 


CO 

w 

g 

§ 


CM 


I 

s 


NO 

C- 

ON 

LfN 

1 

vo 

• 

CO 

(N- 

ON 

o 

rH 

CM 

i 

co 

VO 

-d* 

C~- 

ck 

rH 

CM 

iA 

OH 

1 

VO 

1 

co 
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Input  for  numerical  solution. 


TABLE  34  SPAB  CAP  STRESS  ARISING  FROM  NODE  UNLOADING  IN  EXAMPLE  PROBLEM  WING  (l) 


NODE 

"X" 

NODE  DISPLACEMENTS 

DIRECTION 

"Y"  DIRECTION 

IDENTITY  MAGNITUDE 

IDENTITY 

MAGNITUDE 

[1] 

£lx 

0 

eiy 

0 

[2] 

e2x 

3.723 

CVJ 

o> 

9.024 

[3] 

€3x 

5.519 

€3y 

29.075 

w 

€4x 

6.63O 

[5] 

€5x 

0 

€5y 

46.591 

[6] 

e6x 

-0.1+48 

£6y 

57.182 

[7] 

e7x 

-1.865 

e7y 

66.375 

[8] 

€8x 

0 

e8y 

78.609 

[91 

69x 

0.863 

€9y 

90.742 

[10] 

£10x 

0 

€10y 

99.217 

104 
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TABLE  3I+  SPAR  CAP  STRESS  ARISING  FROM  NODE  UNLOADING 
(Coat'd)  IN  EXAMPLE  PROBLEM  WING  (l) 


SPAR 

CAP 

STRAIN 

"X"  DIRECTION-Ae^  (3) 

"Y"  DIRECTION-Ae  (3) 

....  .. .  Y 

STRESS 

A a 

FORMULA  (2) 

MAGNITUDE 

FORMULA  (2) 

MAGNITUDE 

y 

2-1 

€2x  “  €lx 

10 

_4 

0.3723  x  10 

e  -  e 

5y  ly 

20 

-1+ 

2.3295  x  10 

2050 

3-2 

€3x  £2x 

0.1796  x  lO"^ 

e8y  e5y 

1.6004  x  10"^ 

ll+08 

10 

20 

**~3 

£^x  "  €3x 

10 

0.1111  x  10-14- 

£10y  e8.y 

20 

1.0104  x  10"14' 

899.2 

6-5 

66x  65x 

10 

-0.01+48  x  lO-^ 

6y  2y 

2.1+029  x  lO-14" 

2115 

7-6 

£7x  e6x 

10 

-0.11+17  x  10 

9y  6y 

20 

I.678O  x  lO'^ 

ll+77 

9-8 

e9x  e8x 

10 

O.O863  x  lO-14" 

e7y  "  €3y 

1-8335  X  10-14- 

1613 

NOTE: 


(1) 

(2) 

(3) 


T  =  300 °F. 

ZT 

E  =  8.8  x  10  lbs/sq  in. 

First  digit  of  subscript  denotes  node  point:  the  second  the  direction. 

The  term  e. .  is  a  deformation  term. 
iJ 

3  s^ra^n  ^erm  where  j  denotes  direction  of  strain. 
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TAELE  35  STRESS  INCREMENTS  ARISING  IN  SPAR  CAPS  FROM  NODE  UNLOCKING 

IN  EXAMPLE  PROBLEM  WING 

FOR  DEFORMATION  PATTERN  -  Refer  to  Table  36 


SPAR 

STRAIN 

INCREMENT 

CAP 

DIFFERENTIAL 

INCREMENT 

EFFECTIVE 

INCREMENT 

Ae 

X 

Ae 

y 

Ae 

(1) 

4-7 

e4x  " 

e7x  = 

8.495 

e7y  ' 

'  €4y  = 

13 . 621 

O.7167 

x  10-4 

7-9 

€7x  ' 

G9x  = 

-2.728 

e9y  ' 

'  e7y  = 

4.367 

0.9201 

-4 

x  10 

9-10 

e9x  “ 

eiOx 

=  0.863 

€10y 

1 

=  8.475 

0.3563 

-4 

x  10 

SPAR  CAP  STRESS 

INCREMENT  (2) -An 

4-7 

579-1 

7-9 

743.4 

9-10 

287.9 

NOTES : 

(1)  Ae  =  (aAe  +  bAe  )/e2 

x  y 

where 

a  =  10  in. 

b  =  20  in. 

2  2 

e  =  500  in.  =  a  +  b 

(2)  Ac  =  EAe 

where 

£ 

E  =  8.08  x  10  lbs/sq  in. 

WADD  TR  60-411  PART  I  106 


TABLE  3 6  STRAIN  INCREMENTS  ARISING  FROM  NODE  UNLOCKING 

IN  EXAMPLE  PROBLEM  WING  PANELS 


NOTE: 


NODE  DEFORMATION  (l) 

"X"  DIRECTION  "Y"  DIRECTION 


[1] 

en 

=  0 

e ,  = 

0 

lx 

iy 

[2] 

€2x  = 

=  3.723 

e2y  = 

9.024 

[3] 

£3x  = 

=  5.519 

€3y  - 

29.705 

[4] 

e4x  = 

=  6.630 

52.704 

[5] 

e5x  = 

=  0 

e5y  = 

46.591 

[6] 

e6x  = 

:  -0.448 

€6y  " 

57.182 

[7] 

n 

£ 

u; 

:  -1.865 

11 

£ 

66.375 

[8] 

e8x  “ 

0 

e8y  ~ 

78.069 

[9] 

e9x  = 

0.863 

e9y 

90.742 

[10] 

€10x 

=  0 

eiOy  = 

=  99.217 

l)  Subscript  notation:  number  indicates  effective  node  point; 
letter  indicates  direction. 
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02 


S3 

H 

£ 


VO 

on  rd 


-4-  -4-  -4 

i  i  ' 


o 

o 

O 

I — 1 

r — l 

rH 

X 

X 

X 

CXI 

CVI 

OJ 

vo 

O 

o\ 

VO 

CM 

oo 

oo 

1 - 1 

vo 

OJ 

OJ 

1 — 1 

II 

■ 

II 

rO 

OJ 

OO 

vo 

UJ 

UJ 

UJ 

1 

i 

1 

►>5 

■ 

VO 

tH" 

ON 

uj 

uj 

uj 

O 

o 

+ 

■ft  + 

ft  H- 

>> 

t>5 

r~ 1 

OJ 

LTN 

UJ 

UJ 

UJ 

1 

1 

>> 

>> 

Lf\ 

VO 

00 

uj 

uj 

UJ 

-ft 

1 

-4 

1 

i 

o 

O 

O 

rH 

rH 

i — l 

X 

X 

X 

-ft 

LTN 

o 

o 

on 

CO 

on 

on 

t- 

o 

co 

VO 

• 

• 

rH 

i — ! 

i — 1 

■ 

■■ 

>> 

00 

vo 

uj 

CD 

O 

O  I 

4  1 

ft 

to 

o 

OX 

rH 

id 

UJ 

H 


2 


EH 

CQ 


X 

3 


o 

H 

Eh 

O 


P 


X 


P 


-ft  —t  --I’ 

i  i  i 


o 

o 

o 

rH 

[ — 1 

I — 1 

X 

X 

X 

CO 

On 

c— 

on 

CO 

o 

vo 

i—l 

OJ 

rH 

o 

o 

• 

• 

o 

o 

o 

n 

II 

II 

X 

X 

X 

LTN 

vo 

03 

UJ 

UJ 

U) 

1 

X 

1 

X 

X 

vo 

to¬ 

OX 

UJ 

ll) 

o 

U) 

o 

+ 

CM  + 

CM  + 

X 

X 

X 

r~ 1 

OJ 

LTX 

UJ 

U) 

U) 

| 

1 

X 

1 

X 

X 

0J 

on 

vo 

UJ 

UJ 

UJ 

vo 

t- 

ox 

1 

1 

LTN 

VO 

• 

CO 

1 

1 

OJ 

1 

on 

vo 

1 

rH 

OJ 

LTN 

-ft 

-ft 

1 

o 

-4 

O 

rH 

O 

i — 1 

X 

i — 1 

X 

tH- 

X 

1 — 1 

on 

on 

! - 1 

Cvl 

vo 

1 - 1 

i — l 

CO 

1 - 1 

• 

• 

o 

o 

o 

1 

ii 


X 

oo 

uj 


X 

-H- 

uj 


X 

X 

vo 

CO 

U) 

CD 

o 

o 

CM  1 

CM  1 

X 

X 

t 

OX 

U) 

CD 

O 

OJ 


t>- 

-ft 

I 

OO 


on 

I 

t- 

I 

VO 


o 

I — 1 

I 

ON 

I 

CO 
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NOTE:  (2)  Subscript  notation:  number  indicates  effective  node  point;  letter  indicates  direction 


TABLE  37  DISTORTION  INCREMENTS  ARISING  FROM  NODE  UNLOCKING 
IN  EXAMPLE  PROBLEM  WING  PANELS 


DEFORMATION  INCREMENT  PATTERN 


PANEL  DEFORMATION  INCREMENT -At 
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TABLE  38  STRESS  INCREMENTS  ARISING  FROM  NODE  UNLOCKING  IN  EXAMPLE  PROBLEM  WING  PANELS 
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TABLE  39  RESULTANT  STRESSES  IN  EXAMPLE  PROBLEM  WING 
AFTER  4  MINUTE  CREEP  INTERVAL 
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Figure  23.  A  Possible  Subdivision  of  a  Delta  Wing  into  Sub- structures 
(Schematic) 
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Figure  24.  Node  Point  Location  and  Attachment  Scheme  for  Example 
Problem  Wing 
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Figure  25.  Force  Input  at  Nodes  of  Locked  Example  Problem  Wing 
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REACTIONS 


FORCES  IN  PANEL 


NOTE:  ACTION  ON  LOCKS  DUE  TO  p  -136800  LBS./IN. 


Figure  26.  Node  Forces  on  Triangular  Panels  in  Example  Problem  Wing 
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HORIZONTAL  DISPLACEMENT 
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VERTICAL  DISPLACEMENT 


V2i=  -12 -p 


Figure  27.  Stiffness  Coefficients  for  Beam  Elements  in  Example  Problem 
Wing 
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.97  X  106 


A=  2  SQ.  IN. 

E  =  0.97  X  107 

e=AL± 

L  10 

0  =  ie  =0.97  X  106 
n=  EAe=1.94  X  106 


STIFFNESS  COMPUTATIONS 
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1.94  X  106  b - L  =  10 - 1.94  X  106 

- - I - l-H - - 

STIFFNESS  IN  "X"  DIRECTION  _L„ 

.97  X  106 


STIFFNESS  IN  "Y" 
DIRECTION 


Figure  28.  Stiffness  Coefficients  for  Beam  Elements  in  Example  Problem 
Wing 


WADD  TR  60-411  Part  I 


117 


Figure  29.  Oblique  Frame  Elements  in  Example  Problem  Wing 
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Figure  30.  Stiffness  Coefficients  for  Front  Spar  (Oblique  Bar)  in  Example 
Problem  Wing 
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Figure  31.  Mode  of  Displacement  of  Nodes  of  Example  Problem  Wing  after 
Unlocking 
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Figure  32.  Thermal  Stresses  In  Spar  Caps  of  Example  Problem  Wing 
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Figure  33.  Thermal  Stresses  In  Spar  Webs  of  Example  Problem  Wring 
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Figure  34.  Thermal  Stresses  In  Panels  of  Example  Problem  Wing 
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Creep  Rates  for  75  S-T  Aluminum  Alloy  at  300°F 
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Figure  36.  Creep  Rates  for  75  S-T  Aluminum  Alloy  at  400°F 
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Figure  37.  Creep  Rates  for  75  S-T  Aluminum  Alloy  at  600°F 
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Figure  38.  Creep  Rates  for  75  S-T  Aluminum  Alloy  at  450°F 


WADD  TR  60-411  Part  I 


127 


7 


Figure  39.  Creep  Rate  for  75  S-T  Aluminum  Alloy  at  450°F  at  One  Minute 
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PREDICTION  OF  CREEP  EFFECTS 
IN 

AIRCRAFT  STRUCTURES 
PART  IV 

MATERIALS  CHARACTERISTICS 
IN  RELATION  TO  CREEP 
AND  ITS  PREDICTION 


INTRODUCTION 

In  an  examination  of  characteristics  of  aeronautical  structural  materials  for  use  in 
elevated  temperature  structures,  a  re -consideration  of  basic  governing  conditions 
ruling  uses  of  these  materials  in  past,  "conventional"  temperature  regimes  is  helpful. 

In  former  efforts,  great  emphasis  has  been  put  upon  strength  capabilities  of  airframe 
constructional  materials  that  their  strength  to  weight  ratios  be  kept  low  and  thus 
promote  efficient,  light  weight  vehicles.  This  approach  prevailed  in  a  comparatively 
simple  environment  which  involved  a  narrow  range  of  temperatures  centering  around 
"room"  conditions.  In  this  environment  uses  of  materials  required  high  inherent 
stiffness,  or  moduli  of  elasticity  and  rigidity,  and,  the  preservation  of  stable 
internal  structures  in  the  materials  themselves  to  permit  wide  ranges  of  strength 
utilization  without  impairing  their  stiffness.  In  short  low  temperature  operations 
permitted  conditioning  materials  so  that  high  proportional  limits  could  be  realized 
in  them  and,  what  is  also  important,  could  be  depended  upon  to  serve  without  detri¬ 
mental  change  when  subjected  to  service  environments.  Briefly  the  circumstances 
outlined  allowed  material  use  under  conditions  of  elastic  action  and  this  coincided 
with  structural  design  assuming  small  deflections,  especially  in  structural  elements 
susceptible  to  buckling  effects  where  plastic  action  is  deliberately  restricted  to 
small  amount  s . 

In  the  presence  of  elevated  temperatures,  several  occurrences  take  place  which  run 
counter  to  the  fundamental  ideas  incorporated  in  previous  effort,  namely,  restriction 
of  advantages  arising  out  of  elastic  action.  Among  these  adverse  happenings  are: 

1.  Decay  of  moduli  of  elasticity  and  rigidity  resulting  simply  from  temperature 
increases . 

2.  Decay  in  proportional  limits  resulting  simply  from  temperature  increases  which 
weaken  base  metal  lattice  structures. 

3.  Decay  in  proportional  limits  resulting  from  the  deterioration  of  metal  improve¬ 
ment  steps,  such  as  cold  working,  heat  treat  strengthening,  etc.  ,  resulting 
from  thermal  activation  of  adverse  metallurgical  mechanisms. 

Viewing  these  things  from  the  reverse  side,  it  can  be  seen  that  the  core  idea  of 
alloy  design  for  elevated  temperature  materials  involves  steps  seeking  to  retain 
moduli  and  proportional  limit  values  inherent  in  base  metals  together  with  develop¬ 
ment  of  metal  improvement  schemes  which  require  high  temperatures  to  initiate  dete¬ 
rioration. 

SIGNIFICANT  VARIABLES  INFLUENCING  MATERIALS  PERFORMANCE  IN  HEATED  ENVIRONMENTS 

Four  variables  assert  themselves  in  relation  to  structural  design  when  materials 
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are  gaged  for  performance  in  heated  environments.  They  are  load,  temperature,  time 
and  deformation.  Assuming  that  load  and  temperature  are  items  which  can  he  chosen, 
they  may  he  classed  as  independent  variables:  the  remaining  two  variables,  time  and 
deformation,  then  depend  on  the  other  two,  load  and  temperature.  It  should  be 
realized  that  the  foregoing  classification  of  variables  is  arbitrary  and  is  chosen 
from  convenience.  It  can  readily  be  assumed  also  in  a  structural  design  problem 
that  temperatures  may  be  specified  by  vehicle  velocity  and  time  by  service  require¬ 
ments  and  hence  loads  and  deformations  become  dependent  variables. 

This  example  illustrates  the  point  that  although  operative  variables  in  a  given 
situation  may  be  readily  recognized,  this  recognition  quite  possibly  is  of  relative 
unimportance  uecause  of  the  many  approaches  that  may  conceivably  be  taken  in  putting 
them  into  use.  Accordingly  it  becomes  evident  that  interrelations  of  variables  are 
over-riding  in  importance  to  variables  themselves  when  the  information  utility  is 
approached  objectively. 

THE  DEFORMATION  VARIABLE 


In  approaching  interrelations  of  variables,  knowing  the  relative  importance  of  vari¬ 
ous  factors  entering  into  a  complex  is  valuable  in  essaying  effective  interrelations. 
Inasmuch  as  loads,  temperatures  and  time  are  measurable  quantities  in  terms  of 
pounds,  degrees  and  seconds,  direct  implications  of  magnitudes  are  readily  perceived. 
The  deformation  quantity,  however,  is  directly  measurable  as  a  change  in  length.  But 
because  the  stiffness  of  compression  members,  for  example,  is  not  directly  propor¬ 
tional  to  length  or  length  change,  a  simple  and  direct  relationship  cannot  be  assigned 
to  deformation  factors.  Possibly  the  nearest  approach  to  specification  of  deformation 
variables  lies  in  statement  of  limiting  amounts  by  which  effects  may  be  defined. 

In  a  foregoing  section  it  is  stated  that  "No  plastic  component  of  strain  was  ever 
allowed  previously  in  aircraft  structures  within  limit  load,  except  a  residual 
2  x  10“3  considered  negligible".  In  the  sense  used  above,  the  term  limit  load  ex¬ 
presses  working  load  imposed  on  a  member  and  when  limit  load  Is  multiplied  by  1.5> 
the  product  expresses  ultimate  or  failure  load  for  the  member. 

Considering  the  foregoing,  the  deformation  which  takes  place  in  a  heated  structure 
can  be  visualized  as  the  combination  of  reversible  and  irreversible  deformation 
actions  consisting  of  the  following: 

1.  Reversible  actions. 

a.  Elastic  action  taking  place  in  accordance  with  Hooke's  law  which  is  regu¬ 
lated  in  magnitude  by  the  temperature  present,  and,  in  extent  by  the  con¬ 
dition  of  the  material  at  hand  combined  with  the  extent  to  which  combinations 
of  time  and  temperature  effect  changes  in  the  material  condition. 

b.  Thermal  expansion  taking  place  under  conditions  of  negligible  structural 
restraint.  A  condition  of  negligible  structural  restraint  characterizes 
most  tests  aimed  to  assess  materials  properties.  In  short  the  thermal 
expansion  term  is  cancelled  out  in  tests  but  it  must  be  reckoned  into 
structural  problems. 

2.  Irreversible  actions 

a.  Time  independent  plastic  action  which  appears  immediately  upon  imposition 

of  a  load  on  a  structural  element  or  structure.  In  current  practice  invol¬ 
ving  room  temperature  operations,  "sets"  arising  from  irreversible  plastic 
action  are  commonly  recognized  in  the  term  "permanent  set  .  In  the  light 
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of  current  knowledge,  irreversible  plastic  action  is  governed  in  magnitude 
by  loads  imposed  upon  members  and  in  extent  by  the  rate  at  which  loads 
are  applied. 

b.  Time  dependent  plastic  action  which  usually  is  observed  immediately 
upon  completion  of  time  independent  plastic  action  taking  place  upon 
loading.  This  action  falls  into  three  categories  in  order  of  their 
appearance  with  respect  to  time.  They  are  termed  "primary  creep" 
wherein  decreasing  amounts  of  deformation  occur  as  time  progresses, 
"secondary  creep"  wherein  deformation  accumulates  in  direct  proportion 
to  time,  and  "tertiary  creep"  wherein  increasing  amounts  of  deformation 
appear  as  time  elapses. 

In  view  of  both  the  structural  significance  of  deformation  effects  and  the  involved 
deformation  relationships  occurring  in  actual  operations  of  structures,  the  material 
tests  to  be  described  were  studied  with  particular  attention  to  singular  effects  of 
time,  temperature  and  load  upon  deformation  characteristics. 

"STANDARD"  CREEP  TESTS 


The  "standard"  creep  tests  presently  described  were  performed  with  0.025  in  by  0.125 
in  by  15  in  long  strips  of  Armco  PH  15-7  Mo  Condition  RH-950  stainless  steel  in  a 
special  creep  testing  machine.  In  the  performance  of  these  tests,  the  practice  used 
involved  conventional  procedure  of  attaching  specimens  suitably  equipped  with  tem¬ 
perature  measuring  thermocouples  into  the  machine's  grips,  then  rapidly  applying  heat 
to  the  specimen  assembly  until  the  desired  temperature  was  attained,  and,  then  main¬ 
taining  this  temperature  throughout  the  duration  of  the  test.  Following  application 
of  heat  and  attainment  of  temperature  equilibrium,  loads  were  slowly  applied  so  as  to 
minimize  adverse  effects  arising  from  momentum  and  impact.  As  soon  as  full  load 
application  was  attained,  the  extension  of  specimens  was  observed  as  were  extension 
increments  subsequently  appearing  at  twenty  minute  intervals  throughout  duration  of 
the  test. 

A  noteworthy  item  with  this  test  procedure  is  the  manner  in  which  effects  of  thermal 
expansion  are  cancelled  through  application  of  load  subsequent  to  heating.  The  test 
results,  accordingly,  reflect  the  effects  of  time,  temperature  and  load  in  terms  of 
stress  and  their  effects  upon  only  elastic  and  plastic  deformations. 

In  the  performance  of  these  tests,  temperatures  for  test  were  selected  to  represent 
(l)  a  condition  (700°F)  where  small  metallurgical  change  in  the  test  material  was 
expected;  (2)  a  condition  where  experience  indicated  that  creep  effects  would  be 
marginal  (800°F);  (3)  a  condition  where  moderate  metallurgical  change  in  the  material 
might  be  expected  in  a  long  run  (900°F);  and  (^)  a  condition  where  severe  metallur¬ 
gical  change  in  the  material  could  be  expected  in  a  comparatively  short  run  (1000°F). 
The  loads  (stresses)  applied  were  estimated  from  Larson-Miller  stress  -  rupture  plots 
arbitrarily  proportioned  for  estimating  a  time  of  about  fifty  hours  required  to  pro¬ 
duce  about  one  percent  total  creep  deformation.  The  specific  conditions  for  individual 
tests  are  given  in  Table  40.  Figures  40  to  55  show  results  from  these  tests.  Figure 
56  shows  details  of  deformation  patterns  taken  during  these  tests.  Figure  57  shows 
the  elastic  recovery  and  permanent  set  patterns  taken  by  the  test  specimens  upon  un¬ 
loading.  Figures  58  and  59  compare  the  total  time  independent  deformation  comprised 
of  combined  elastic  and  time  independent  inelastic  deformation  taking  place  upon 
loading  with  the  elastic  recovery  appearing  upon  unloading. 

LIMITING  DEFORMATION:  STRESS  STRAIN  DIAGRAM  APPROACH 


In  the  study  of  data  showing  deformation  information  confusion  may  arise  in  distin¬ 
guishing  between  the  importance  of  those  deformations  arising  from  the  simple 
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application  of  load  and  those  permanent  deformations  remaining  after  loads  are  re¬ 
leased.  In  practice  illustrated  by  ordinary  stress  strain  curves,  the  actual  deflec¬ 
tions  taking  place  upon  loading,  especially  those  occurring  below  the  yield  points, 
are  of  minor  concern  in  that  they  do  not  serve  to  establish  limits  for  material 
usages.  The  limiting  factor  is  the  "offset"  or  permanent  set  arising  from  plastic 
flow.  In  a  sense  loading  deflections  are  indicative  of  causes  and  permanent  sets 
are  indicative  of  results. 

In  the  study  of  creep  curves  the  subtlety  of  the  interrelations  of  cause  and  effect 
assert  themselves,  especially  when  such  study  is  approached  with  a  specific  limit¬ 
ing  permanent  set,  or  what  is  analagous,  plastic  strain,  in  mind.  The  first  thing 
that  appears  in  creep  testing  is  the  deformation  that  occurs  upon  loading  and  is 
independent  of  time.  Where  a  limiting  plastic  strain  of  .002  prevails,  this  limit 
can  be  attained  almost  instantaneously  upon  load  application  providing  stresses  are 
equivalent  to  the  yield  stress.  In  the  event  it  should  be  desired  that  permanent 
sets  arising  from  loading  alone  be  eliminated,  first  thought  would  lead  to  use  of 
loads  that  produce  stresses  at  or  below  the  proportional  limit.  By  way  of  the  ex¬ 
ample  of  Armco  PH  15-7  Mo,  Condition  RH  950  tested  at  800°F,  the  yield  strength 
given  by  the  producer  is  14-9  K.S.I.  and  the  proportional  limit  is  92.5  KSI.  The 
difference  in  magnitude  between  these  two  values  emphasizes  the  importance  of  good 
working  stress  level  selection  especially  as  a  three  way  compromise  between  struc¬ 
tural  weight,  the  initial  deformation  upon  loading  and  subsequent  creep  deformation 
is  involved  in  the  selection. 

In  relation  to  the  foregoing  example,  it  is  of  importance  that  the  deformation  pre¬ 
sent  under  actual  load  amounts  to  .040  ( 4.0$>  deformation)  strain  at  the  proportional 
limit  and  .079  (7 -9$  deformation)  at  yield  stress  as  indicated  by  the  Armco  Techni¬ 
cal  Data  Manual  information  referencing  to  tension  test  (stress-strain  diagram)  data. 

LIMITING  DEFORMATION:  CREEP  (TIME-DEFORMATION)  DIAGRAM  APPROACH 


In  relating  creep  curve  information  given  in  Figures  44  through  47  with  observa¬ 
tions  given  above  and  arrived  at  from  study  of  the  Armco  Technical  Data  Manual,  the 
disparity  of  deformation  information  from  the  two  sources  is  evident.  Whereas 
stress-strain  diagrams  from  tension  tests  indicate  .079  strain  at  149  KSI  yield 
stress,  the  creep  diagrams  showed  about  .028  and  .025  strain  at  170  KSI  stress  and 
about  .029  and  .026  strain  at  175  KSI  stress,  both  of  which  exceeded  the  yield 
stress  and  proportional  limit  stresses  reported  from  Armco  tension  tests  by  wide 
margins.  Although  rigorous  explanation  of  these  differences  is  not  immediately 
available,  mode  of  loading  suggests  itself  as  a  possible  contributor  to  this  situa¬ 
tion.  In  tension  testing  the  loads  are  gradually  applied  over  an  appreciable  length 
of  time,  and,  in  creep  tests  loading  is  accomplished  as  rapidly  as  possible  and  in 
one  increment.  In  passing  it  may  be  surmised  that  the  strain  differences  observed 
upon  loading  the  two  pairs  of  creep  specimens,  .028  and  .025  strain  at  170  KSI 
stress  and  .029  and  .026  at  175  KSI  stress,  arose  from  variations  in  speed  of  loading. 

STRAIN  MODES  APPEARING  IN  CREEP  TESTS 


Figure  56  relates  the  outcome  of  attempts  to  assess  creep  curves  in  terms  of  strain 
modes,  especially  up  to  the  time  where  secondary  creep  asserts  itself.  The  elastic 
strain  shown  in  Figure  56  is  not  elastic  deformation  taking  place  up  to  the  propor¬ 
tional  limit:  it  is  the  elastic  component  of  total  strain  taking  place  upon  loading. 
It  is  the  quotient  of  the  applied  stress  divided  by  that  modulus  of  elasticity  given 
for  the  appropriate  temperature  in  the  Armco  Data  Manual.  The  time  independent 
plastic  strain  shown  in  Figure  56  is  the  plastic  component  of  the  total  strain  taking 
place  upon  loading.  It  is  the  difference  between  the  total  strain  and  the  elastic 
strain  component .  The  primary  creep  strain  shown  in  Figure  5 6  is  strain  occurring 
between  time  of  completion  of  loading  to  beginning  of  secondary  creep. 
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Figure  57  shows  results  from  unloading  the  creep  specimens  at  the  end  of  test  run. 

The  elastic  recovery  shown  is  the  contraction  in  specimen  length  observed  upon  re¬ 
moval  of  the  load.  The  permanent  set  indicated  in  Figure  57  is  that  dilation  exper¬ 
ienced  by  the  specimen  during  test  which  remained  unchanged  after  unloading. 

Figures  58  an<i  59  compare  the  time  independent  deformations  occurring  upon  loading 
with  the  elastic  recovery  observed  after  unloading.  These  charts  show,  in  general, 
that  the  energy  input  into  a  specimen  in  the  form  of  elastic  and  plastic  action  is 
not  entirely  converted  into  potential  energy  in  the  form  of  stored  elastic  action. 

The  energy  losses  suggested  by  the  differences  in  height  of  the  two  bars  represent¬ 
ing  each  specimen  represent  a  part  of  the  energy  input  to  specimens  which  was  used 
up  in  plastic  deformation  action. 

The  point  is  that  loading  energy  input  is  used  up  in  both  storage  of  potential  energy 
in  a  test  bar  and  in  causing  it  to  deform  in  a  plastic  manner.  In  other  words  part 
of  the  energy  is  recoverable  and  that  which  is  not  can  contribute  to  structural 
damage  prior  to  the  onset  of  time  -  dependent  creep  damage  as  shown  in  Figures  60 
and  6l.  In  these  figures  the  time  independent  (elastic  and  inelastic)  deformation 
converted  to  permanent  set  upon  loading  represents  the  difference  between  the  total 
permanent  set  arising  from  the  combination  of  inelastic  and  creep  deformation,  and, 
the  creep  deformation.  The  charts  depicted  in  Figures  60  and  6l  show  that  at  all 
four  temperatures  of  test  the  portion  of  permanent  set  acquired  during  loading  was 
generally  of  a  magnitude  that  constituted  structural  damage  conditions.  Briefly, 
the  pronounced  tendency  for  time  independent  plastic  strain  to  cause  permanent  set 
in  excess  of  .002  (0.2 $  offset)  strain  was  evidenced.  In  view  of  this,  together 
^ith  the  fact  that  the  loads  used  were  greatly  in  excess  of  both  the  proportional 
and  yield  limits,  the  desirability  for  information  relative  to  those  loads  which 
would  result  in  negligible  structural  damage  asserted  itself. 

INCREMENTAL  LOADING  TESTS 


In  view  of  the  importance  of  stress  level  in  creep  tests,  especially  because  this 
level  influences  magnitude  of  time  independent  plastic  deformation  experienced  upon 
loading,  incremental  loading  tests  were  made  to  determine  approximate  stress  levels 
at  which  time  independent  plastic  deformation  asserts  itself,  its  increase  with 
respect  to  increasing  stress,  and,  its  repeatability  upon  complete  recycling. 

The  incremental  loading  tests  were  made  with  Armco  PH  15-7  Mo,  Condition  RH-950 
stainless  steel.  These  tests  were  run  in  Arcweld  Model  JE  creep  test  machines  using 
Arcweld  Model  400  extensometers  to  observe  deformations  taking  place.  Tests  were 
run  under  the  limiting  conditions  of  test  shown  in  Table  k-0.  Specific  test  per¬ 
formance  was  as  follows:  first,  a  specimen  having  a  .505  in.  wide  x  .025  in.  thick 
x  2  in.  gage  length  test  section  conforming  to  Federal  Standard  151A  was  loaded  to 
10$  of  the  limiting  stress;  then  this  load  was  released;  in  a  second  step  the  speci¬ 
men  was  loaded  to  20$  of  the  limiting  stress,  and  then  this  load  was  released. 
Following  this,  successive  loadings  at  30$,  40$  and  so  on  up  to  100$  of  limiting 
stress  were  applied  and  released  to  complete  a  single  "step"  loading  cycle.  In  each 
test  run  with  a  single  specimen,  five  complete  cycles  were  run  in  order  to  examine 
repeatability  of  performance  under  repetitive  conditions  of  test. 

INITIAL  INCREMENTAL  LOADING  CYCLES 

Figure  62  summarizes  the  findings  from  a  single  initial  loading  cycle.  In  this 
chart  the  line  AB  shows  the  total  deflections  obtained  throughout  the  loading  pro¬ 
gram.  The  line  ACD  shows  the  extensions  observed  upon  release  of  each  incremental 
load.  The  area  ABCD  shows  a  field  wherein  elastic  action  activity  is  manifested, 
the  points  along  the  line  AB  representing  extension  upon  application  of  load  and 
the  line  CD  representing  the  residual  extension  accruing  upon  release  of  load.  The 
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area  CDE  represents  a  field  of  permanent  set  accumulation  wherein  a  portion  of  the 
energy  input  arising  from  loading  is  converted  into  set  rather  than  into  stored 
potential  energy  which  is  reflected  in  elastic  recovery. 

If  the  permanent  set  is  considered  as  representing  a  plastic  component  of  strain, 
it  may  be  looked  upon  as  damaging  so  far  as  an  airframe  structure  is  concerned.  In 
the  instance  shown,  the  plastic  component  appears  at  a  stress  of  about  92.5  KSI, 
the  porportional  limit  stress  shown  in  the  Armco  Technical  Data  Manual.  In  addi¬ 
tion  the  stress  required  to  produce  0.2 'jo  deformation  is  l66  KSI,  a  number  in  excess 
of  the  149  KSI  yield  strength  figure  given  by  Armco  for  700 °F  yield  strength.  This 
latter  number,  166  KSI,  suggests,  in  comparison  with  the  former,  149  KSI,  an  effect 
arising  from  loading  rate.  In  the  case  of  incremental  tests,  loading  was  accom¬ 
plished  in  about  5  to  10  seconds.  In  tension  testing  about  75  to  90  seconds  are 
ordinarily  required  to  attain  the  yield  load. 

SECONDARY  INCREMENTAL  LOADING  CYCLES 


In  the  sense  used  herein,  the  term,  secondary  incremental  loading  cycle,  is  used  to 
denote  any  cycle  exclusive  of  an  initial  incremental  loading  cycle.  Figure  63  shows 
such  a  secondary  cycle  as  an  example.  Specifically  the  cycle  shown  in  Figure  63 
Immediately  followed  that  shown  in  Figure  62  during  test  and  both  represent  a  single 
test  run.  In  Figure  63  the  zero  deformation  shown  at  the  origin  of  coordinates 
actually  represents  final  permanent  set  attained  at  the  end  of  the  initial  incre¬ 
mental  loading  cycle.  Also,  the  line  AB  indicates  the  total  deflections  observed 
throughout  the  secondary  cycle  loading  program.  The  area  ABC  shows  a  field  wherein 
elastic  action  activity  is  manifested.  The  points  along  the  line  AB  represent 
extention  upon  application  of  load  and  the  line  AC  represents  the  residual  extention 
resulting  from  release  of  load.  The  area  ACD  represents  a  field  of  permanent  set 
accumulation  wherein  a  portion  of  the  energy  input  arising  from  loading  is  converted 
into  set  rather  than  stored  potential  energy  which  is  reflected  in  elastic  recovery. 
The  schematic  nature  of  the  line  AC  is  to  be  noted  and  this  circumstance  arises 
from  measurement  of  dimensions  in  the  order  of  0.0001  to  0.002  inch  with  mechanical 
instrumentation.  Accordingly  firm  evidence  is  not  believed  to  be  presented  in 
single  secondary  cycle  results  but  trends  may  be  roughly  deduced  from  multiple  cycle 
test  trends. 

In  the  example  shown  in  Figure  63  the  predominance  of  elastic  action  is  fairly 
evident.  The  plastic  action  indicated  appears  as  a  more  or  less  linear  growth  sug¬ 
gestive  of  secondary  creep.  Evidence  of  time  independent  plastic  action  or  primary 
creep  action  is  not  suggested  since  the  non-linear  relations  are  not  evidenced. 

The  suggestion  that  presents  itself  however  is  that  lattice  displacements  accom¬ 
plished  during  initial  cycle  loadings  were  adequate  to  offset  those  displacements 
ordinarily  expected  in  "normal"  primary  creep.  Consequently  the  material  may  be 
"triggered"  to  go  into  secondary  creep  during  the  second  cycle  upon  input  of  suffi¬ 
cient  incremental  energy  to  initiate  secondary  creep  activation. 

The  interesting  implication  from  the  foregoing  is  a  suggested  mechanism  for  "stab¬ 
ilizing"  a  metal  for  "creep"  service  by  hot  stretching,  for  example.  Second 
thought,  however,  suggests  that  such  an  operation  implies  conditioning  fully  hardened 
material.  The  prospects  for  fabricating  such  a  material  into  airframe  configurations 
appear  remote  and  presently  seem  beyond  ordinary  fabricating  capabilities. 

INCREMENTAL  LOADING  TESTS:  Repetitive  Loadings 

Figure  64  schematically  diagrams  the  program  approach  used  to  obtain  the  results 
shown  in  Figures  65,  66  and  67.  These  tests  were  run  at  a  constant  temperature  of 
700°F  with  a  full  load  of  185  KSI.  Figures  68  and  69  display  results  from  tests 
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run  at  800°F  with  a  full  load  of  175  KSI.  The  800°F  test  condition  is  represented 
by  a  single  test.  Two  tests  were  run  but  extensive  extensometer  slippage  negated 
the  second  test  run  which  was  cast  out  during  data  reduction.  The  results  from 
test  run  at  900°F  and  135  KSI  full  load  are  given  in  Figures  70>  71  and  72  and  those 
results  from  tests  run  at  1000°F  and  77  KSI  full  load  are  shown  in  Figures  73,  7^ 
and  75* 

Figures  65,  66,  68,  70  and  71*  and,  73  and  7^-  depict  the  deformation  occurring  during 
the  loadings  at  700°,  800°,  900°  and  1000 °F  respectively.  These  curves  resemble 
stress-strain  curves  and  show  the  work  hardening  taking  place  during  the  initial 
cycle  followed  by  elastic  working  of  this  cold  worked  material  in  subsequent  cycles 
The  chief  item  of  interest  shown  in  these  curves  is  the  "modulus  line"  shown  in  the 
straight  line  portions  of  component  portions  of  the  several  curves.  In  essence 
these  "modulus"  lines  serve  only  to  reiterate  the  well  understood  and  evaluated 
inverse  change  of  modulus  of  elasticity  as  temperature  increases.  Figures  67,  69, 

72  and  75  show  the  deformations  occurring  upon  unloading  subsequent  to  each  in¬ 
cremental  loading.  These  curves  reassert  the  damaging  effect  of  an  initial  loading 
cycle  and  at  least  in  the  instance  of  the  tests  run  at  700°F,  800°F  and  900°F,  all 
temperatures  below  the  950  °F  aging  temperature  used  to  produce  Condition  RH-950 
in  Armco  PH  15-7  Mo,  show  the  profound  influence  of  the  time  independent  plastic 
strain  in  the  deformation  outcome  as  it  relates  to  the  "secondary  creep"  taking 
place  after  the  initial  loading  cycle  has  done  its  work. 

Figure  76  compares  loading  strains  developed  by  incremental  loading  to  a  pre¬ 
determined  stress  level  at  the  end  of  an  initial  loading  cycle  with  those  developed 
at  the  end  of  a  fifth  incremental  loading  cycle.  In  actuality  those  strains  shown 
as  appearing  at  the  end  of  the  fifth  loading  cycle  represent  the  permanent  set  accrued 
during  the  second,  third,  fourth  and  fifth  loading  cycles.  The  purpose  of  this 
comparison  is  to  emphasize  the  strong  influence  exerted  by  original  loading  in 
relation  to  the  somewhat  minor  influence  exerted  by  succeeding  loadings.  For 
additional  comparative  information  elastic  action  components  calculated  from 
published  elastic  modulus  data  are  shown  in  Figure  76. 

Figure  77  compares  strains  observed  upon  unloading  the  incremental  loading  tests 
and  again  it  compares  strains  observed  at  the  end  of  the  first  with  those  observed 
at  the  end  of  the  fifth  loading  cycles.  In  this  chart,  the  total  deformation, 
which  includes  both  permanent  set  and  elastic  recovery  is  compared  with  the  elastic 
recovery  observed  in  each  test  case.  The  high  permanent  set  generation,  especially 
at  700°,  800 °F,  and  900°F,  certainly  indicates  the  influence  of  high  stress  level 
in  promoting  early  damage.  However,  as  previously  discussed,  incremental  step-load¬ 
ing  tests  indicate  the  ratio  of  permanent  set  to  elastic  recovery  is  possibly 
amenable  to  control  providing  low  stress  levels  in  the  neighborhood  of  elastic 
limit  stress  are  compatible  with  structure  weights.  In  the  comparisons  shown  in 
Figure  77  the  important  influence  of  the  first  loading  is  reasserted  in  that  the 
permanent  set  produced  in  this  step  approached  that  produced  in  the  succeeding  four 
loading  cycles. 

Figures  78  find  79  compare  the  total  strain  observed  as  a  result  of  incremental  load¬ 
ing  with  the  elastic  recovery  observed  upon  load  removal.  These  charts  also  compare 
the  calculated  elastic  action  strain  component  anticipated  upon  loading  with  the 
strain  recovery  observed  upon  unloading.  In  these  charts  effects  arising  in  the 
first  loading  cycle  are  shown  separately  from  those  accrued  in  the  second,  third, 
fourth  and  fifth  loading  cycles.  This  was  done  to  reassess  the  importance  of 
initial  loading  action.  This  point  is  brought  out  by  the  charts  which  essentially 
show  that  severe  elastic -plastic  deformation  experienced  in  initial  loading  serves 
to  condition  material  by  a  work  hardening  process  to  subsequently  serve  as  an 
elastic  material  when  subjected  to  additional  loadings.  Clean-cut  comparisons 
relating  elastic  action  input  and  output  are  not  shown  in  these  charts.  It  is  felt 
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that  to  obtain  these  comparisons,  which  appear  fundamental  to  a  sound  understanding 
of  the  efficiency  of  elastic  action  and  "frictional"  effects  encountered  in  its 
operation,  critical  experiment  with  highly  refined  apparatus  is  required.  This 
approach  appears  highly  desirable  because  in  repeated  loading  situations  it  seems 
quite  probable  that  in  addition  to  gliding  movement  usually  associated  with  creep 
observed  in  conventional  tests,  an  additional  "friction"  or  hysteresis  entering 
with  re-cycled  loadings  is  suggested.  This  implies  desirability  for  knowledge  of 
mode  and  rate  of  accrual  of  these  "run-down"  factors. 

REPEATED  LOADING  TESTS 


In  the  performance  of  incremental  loading  tests  the  gradual  build-up  of  stress 
level  was  considered  a  means  for  accumulating  "cold-working"  effects  at  a  slow 
rate.  On  the  other  hand  the  sudden  build-up  of  stress  level  in  a  single  loading 
was  considered  a  means  for  accumulating  "cold  work"  at  a  greater  rate  than  that 
experienced  in  incremental  load  build-ups.  Also,  because  of  singular  importance 
asserted  by  initial  incremental  loading  cycles  in  relation  to  subsequent  loading 
cycles,  time  suggested  itself  as  an  influence  and  it  appeared  desirable  to  learn 
from  repeated  loading  tests  whether  rapidity  of  load  application  was  conducive  to 
build-up  of  "cold  working"  effects  in  load  cycles  other  than  the  initial  one. 

The  repeated  loading  tests  discussed  here  were  made  with  Armeo  PH  15-7  Mo,  Condi¬ 
tion  RH  950  stainless  steel.  These  tests  were  run  in  Arcweld  Model  J.E.  creep 
machines  using  Arcweld  Model  400  extensometers  to  observe  deflections.  Tests  were 
run  under  the  conditions  shown  in  Table  40.  Specific  test  performance  consisted  of 
the  following:  first,  a  specimen  with  a  0.505"  dia.  x  0.025"  thick  x  2"  long  gage 
length  conforming  to  Federal  Standard  15IA  was  loaded  to  the  limiting  stress;  then 
the  load  was  released.  In  each  run  with  a  single  specimen,  five  complete  loadings 
were  successively  applied  and  released  to  provide  for  data  for  examining  repeata¬ 
bility  of  performance  in  the  various  specimens  under  repetitious  loading. 

Figure  80  schematically  diagrams  the  approach  used  to  obtain  the  results  shown  in 
Figures  8l  through  84  inclusive.  In  these  charts  the  curves  drawn  are  faired  and 
do  not  reflect  observations  except  at  the  data  points  indicated.  Figures  85 
through  88  reflect  the  conditions  obtaining  with  respect  to  permanent  set-elastic 
recovery  relations  at  the  end  of  each  loading  cycle  throughout  the  entire  series  of 
tests.  The  permanent  set  shown  in  these  charts  is  cumulative  set  acquired  as  tests 
proceeded  through  cycling  The  amounts  of  permanent  set  shown  are  large  and,  as 
the  charts  show,  disproportionate  with  the  elastic  recovery.  As  mentioned  in  dis¬ 
cussion  relative  to  incremental  loading  tests,  the  indications  are  that  to  arrive 
at  tolerable  amounts  of  permanent  set  low  stresses  are  necessary. 

Figures  89  through  92  chart  changes  in  permanent  set  and  elastic  recovery  observed 
throughout  repeated  loading  tests  These  charts  reiterate  the  observation  that 
initial  loading  produced  the  greatest  and  determining  set  observed  in  each  test 
case  studied.  The  data  charted  in  Figures  8y  through  92  do  not  clearly  show 
whether  loadings  subsequent  to  initial  loading  cause  definite  increases  on  de¬ 
creases  in  either  permanent  set  or  elastic  recovery.  The  intimation  of  the  data, 
however,  is  that  at  least  a  slight  tendency  for  both  permanent  set  and  elastic 
recovery  to  increase  as  load  cycling  progresses  is  present.  The  outcome  of  these 
tests  suggests  two  things:  the  need  for  critical  experiment  in  order  to  more 
accurately  define  the  trends  for  change  in  permanent  set  and  elastic  recovery 
magnitudes;  and,  the  need  for  extended  experiment  to  determine  whether  the  simul¬ 
taneous  increase  in  permanent  set  and  elastic  recovery  persist  beyond  a  primary 
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creep  phase  of  action.  Figures  93  and  94  show  comparisons  attempting  to  relate 
permanent  set  and  elastic  recovery  observed  subsequent  to  sequenced  load  appli¬ 
cations  with  the  temperatures  at  which  the  load  applications  were  made.  Because 
the  stresses  under  which  these  tests  were  carried  out  apparently  have  little  re¬ 
lation  to  each  other,  clean-cut  interpretation  of  these  chants  is  difficult.  The 
striking  impression  gleaned  from  these  charts  is  that  in  order  to  effect  correla¬ 
tion  of  test  variable  interrelations  a  base  line  of  departure  in  the  form  of  some 
physical  event  is  required  to  supplant  arbitrary  test  variable  selection. 

TRANSIENT  HEATING  TESTS 


The  introduction  of  repeated  heatings  into  creep  tests  is  recognized  as  the  intro¬ 
duction  of  a  complex  rather  than  a  simple  variable.  This  situation  distinguishes 
tests  of  this  nature  from  those  made  at  a  constant  temperature  for  several  reasons, 
among  which  are: 

1.  Decay  in  the  modulus  of  elasticity  in  direct  proportion  to  temperature 
increase 

2.  Decay  in  the  elastic  limit  in  direct  relation  to  temperature  increase 

3 .  Introduction  of  displacements  arising  from  thermal  expansion  effects 
which  vary  from  linearity  with  respect  to  temperature  ranges  present. 

4.  Non-linearity  of  the  heat  or  temperature  increase  which  apparently  proceeds 
in  accordance  with  an  exponential  law. 

The  transient  heating  tests  described  were  made  with  0.025"  thick  x  0.125"  wide  x 
3"  long  gage  length  strips  of  Armco  PH  15-7  Mo,  Condition  RH  950  stainless  steel. 

In  the  conduct  of  these  tests,  the  desired  stress  level  was  obtained  by  dead 
weight  loading  cold  specimens  prior  to  any  heating.  After  loading  heat  was  applied 
by  resistance  heating  a  0.500"  outside  diameter  x  0.010  wall  Inconel  tube  sur¬ 
rounding  the  specimen.  This  arrangement  provided  a  combination  of  radiant  and  con¬ 
vective  heat  transfer  for  specimen  heating  and  resulted  in  a  requirement  for  from 
about  2.75  to  4.0  minutes  to  heat  an  individual  specimen  from  room  to  test  tempera¬ 
ture.  Upon  attainment  of  desired  temperature  in  a  specimen,  the  heater  was  de¬ 
energized  and  allowed  to  cool  to  room  temperature.  This  cooling  consumed  about 
55  minutes  in  each  test  case.  In  each  test  case  the  heating  program  involved  five 
successive  heatings  and  coolings. 

Figure  95  schematically  describes  the  procedure  used  to  obtain  the  data  shown  in 
Figure  96.  This  chart  shows  plots  of  data  points  obtained  at  20-second  intervals 
during  initial  transient  heating  of  each  test  specimen.  Aside  from  display  of 
heating  curves  falling  into  a  scatter  band  delineated  in  the  greater  part  by 
curves  governed  by  the  temperature  level  and  corresponding  heating  rates,  this 
data  plot  yields  little  information.  Figures  97  and  98  are  replots  of  the  data 
shown  in  Figure  96.  These  plots  by  throwing  the  data  into  logarithmic  form  bring 
out  the  duplex  nature  of  the  curves  shown  in  Figure  96  and  thus  bring  out  indi¬ 
cation  of  a  transition  which  suggests  change  in  the  mode  of  response  of  the  material 
to  its  environment .  At  this  stage  of  investigation  it  does  not  appear  wise  to 
attempt  assessment  of  this  transition  in  terms  of  the  stress,  times  and  tempera¬ 
tures  involved.  However  it  does  appear  that  the  transition  points  may  be  viewed 
as  an  indication  of  shift  from  a  state  where  elastic  action  predominates  to  one 
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where  plastic  action  predominates.  Figure  99  attempts  to  relate  transition  points 
with  temperatures  at  which  they  occurred.  Because  of  the  undetermined  relation 
between  the  indicated  transition  and  an  "effective  elastic  limit,"  it  is  difficult 
to  discern  trends  from  the  meager  data  at  hand.  However^  judging  from  the  fact 
that  the  strain  at  which  the  elastic  limit  occurs  as  temperatures  increase  in 
conventional  short  time  tension  tests,  it  is  surmised  that  the  data  imply  a  down¬ 
ward  trend  as  indicated  in  Figure  99- 

Figure  100  displays  data  similar  to  that  shown  in  Figure  96  and  represents  tran¬ 
sient  heating  tests  run  at  700°F  at  various  stress  levels.  These  data  also  present 
the  confused  situation  reflected  in  Figure  96  which  is  discussed  above. 

Figures  101  and  102  display  replots  of  the  data  shown  in  Figure  100.  These 
logarithmic  plots  also  show  the  existence  of  transition  points  in  the  time-deforma¬ 
tion  curves  which  suggest  the  shift  from  elastic  to  plastic  action  fields  already 
discussed.  Figure  103  relates  the  occurrence  of  the  transition  points  with  res¬ 
pect  to  the  strains  at  which  they  occur  at  differing  stress  levels.  The  upward 
trend  indicated  in  Figure  103  appears  to  be  consistent  with  observations  connected 
with  short  time  tension  test  data. 

The  foregoing  experimentation  indicates  in  its  way  a  requirement  for  critical 
experiment  which  provides  continuous  records  of  concurrently  occurring  strains 
and  temperatures  which  appear  in  transient  heating  tests.  Such  experiment  will 
lead  to  the  determination  of  "effective  elastic  limits"  operative  under  these 
environmental  conditions. 

Figures  104  through  110  display  results  obtained  from  five  cycle  transient  heating 
runs.  In  all  of  these  charts,  the  dominant  position  of  the  initial  heating  cycle 
stands  out  both  as  it  affected  the  deformations  appearing  upon  loading  and  upon 
release  of  load.  A  general  tendency  for  a  steady  permanent  set  condition  to  es¬ 
tablish  itself  in  the  initial  cycle  is  generally  evidenced  in  these  curves .  The 
permanent  set  observed  beyond  the  initial  cycle  accumulated  in  a  linear  manner  in 
general.  However  the  outcome  shown  in  Figure  107  was  such  that  it  suggested  the 
possible  influence  of  a  "primary"  creep  action.  This  primary  creep  action  first 
occurred  at  the  highest  stress  level  used  at  700°F  and  is  of  extreme  interest 
because  it  hints  at  possibility  of  a  "strain  aging"  mechanism  coming  into  play. 

In  general,  however,  the  deformations  occurring  upon  "loading"  appeared  to  coin¬ 
cide  in  curve  form,  shape  and  trend  with  the  permanent  set  curves.  The  gradual 
strain  increase  indicated  both  upon  "loading"  and  in  permanent  set  observed  upon 
"load"  release  was  indicative  of  an  action  akin  to  secondary  creep  action  and 
probably  reflected  a  time  effect  introduced  in  the  normal  course  of  heating  and 
cooling . 

PULSED  CYCLE  HEATING  TESTS 


The  complex  conditions  arising  from  heatings  of  specimens  maintained  at  a  constant 
stress  level  in  pulsed  cycle  heating  tests  are  the  same  as  those  that  apply  in 
transient  heating  cycles.  The  pulsed  cycle  heating  tests  were  accomplished  under 
the  same  conditions  as  those  used  for  transient  heating  tests.  The  distinguishing 
difference  between  the  two  series  of  tests  appeared  in  the  maintenance  of  the 
pulsed  cycle  tests  at  temperature  for  specific  time  intervals  prior  to  initiation 
of  cooling.  Figure  111  schematically  diagrams  the  procedure  followed  in  pulsed 
cycle  heating  tests  and  Figures  112  through  115  plot  the  permanent  set  trends  for 
each  test. 
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Figures  ll6,  117,  and  118  compare  the  effects  of  differing  stress  level  at  the 
various  temperatures  used  in  the  5-minute  pulse  pulsed  heating  tests.  These  curves 
reiterate  the  generalities  found  in  "conventional"  creep  tests  in  that  they  show 
the  existence  of  primary  and  secondary  stages  of  creep  together  with  the  tendency 
for  higher  stresses  to  accelerate  both  primary  and  secondary  creep.  Figure  119  dis¬ 
plays  the  effects  of  differing  stress  level  at  1000°F  upon  tests  using  a  25-minute 
heating  pulse.  Figures  120,  121  and  122  show  the  same  kind  of  comparative  data 
for  those  tests  run  with  a  125-minute  heating  pulse.  The  curves  shown  in  Figures 
119  through  122  reiterate  the  generalizations  given  in  relation  with  Figures  116 
through  118. 

Figures  123  through  126  compare  the  effects  of  varying  pulse  time  periods  at  equi¬ 
valent  stress  levels  and  temperatures.  The  characteristic  shown  in  these  compari¬ 
sons  is  the  tendency  for  creep  strains  to  increase  markedly  together  with  an  accel¬ 
erating  creep  rate  as  the  time  pulse  duration  increases.  Ready  explanation  of  this 
feature  does  not  appear  to  arise  from  total  energy  input  since  the  number  of  cycles 
of  heating  at  the  125  minute  pulse  interval  are  less  than  those  at  a  5  minute  pulse 
interval.  This  arises  because  the  charts  show  time  only  in  terms  of  pulse  times 
exclusive  of  the  times  during  the  heatings-up  and  cooling-down  where  variable  heat 
input  is  present.  About  the  only  surmise  that  can  be  advanced  without  detailed 
study  is  that  a  finite  time  subsequent  to  the  attainment  of  a  given  temperature  is 
required  to  activate  the  creep  processes  and  set  them  in  motion  at  a  maximum  steady 
rate  consistent  with  the  environmental  conditions  at  hand. 

Figures  127  and  128  compare  the  effects  of  the  various  temperatures  at  equivalent 
pulse  times  upon  the  observed  strain.  Largely  because  of  the  lack  of  comparability 
between  the  various  stress  levels  used  at  the  differing  temperature  levels,  these 
charts  do  not  serve  to  reveal  contrasting  features. 

COMPARISON  OF  EFFECTS  OF  "CONVENTIONAL"  CREEP  WITH  THOSE  OF  "PULSED  STEP- LOAD" 

CREEP 

These  test  were  made  to  compare  the  effects  of  conventional  creep,  that  is,  creep 
which  takes  place  under  conditions  of  steady  stress  application  at  constant  temper¬ 
ature,  with  those  effects  observed  as  a  result  of  creep  under  pulsed  step  loading 
at  constant  temperature.  The  pulsed  step  loading  used  consisted  of  applying  a  50KSI 
load,  a  stress  level  slightly  above  that  at  which  creep  activation  was  expected  at 
the  temperature  used,  maintaining  this  stress  for  fifteen  minutes  and  then,  without 
load  release,  increasing  the  stress  to  70KSI.  Upon  reaching  the  70KSI  stress,  this 
stress  level  was  maintained  for  fifteen  minutes  whereupon  the  loads  were  released 
to  complete  a  single  step  load  pulse  cycle.  In  the  test  programming  the  effects 
resulting  from  sixteen  step  load  pulses  which  were  equivalent  to  eight  hours  of 
creep  exposure  were  sought.  These  effects  were  for  comparison  with  those  arising 
from  four  hours'  exposure,  separately  applied,  at  each  of  the  two  stress  levels 
mentioned.  In  the  interpretation  of  test  outcomes,  comparisons  of  effects  arising 
from  four  hours  to  70KSI  added  to  those  arising  from  four  hours  at  50KSI  were  to  be 
made  with  those  arising  from  eight  hours  of  pulsed  loadings. 

The  tests  mentioned  were  run  in  a  Tinius  Olsen  Electromatic  Universal  Testing 
Machine  Olsen  S-l  recording  extensometers  acting  in  conjunction  with  an  Olsen 
Model  1050  recorder  were  used  to  observe  the  strains  taking  place.  The  universal 
testing  machine  itself  employs  a  Thymotrol  speed  control  unit  which  provides 
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infinitely  variable  speeds  and  this  control  provided  means  for  slowing  the  testing 
machine  motor  to  a  point  where  only  enough  "take  up"  to  accomplish  constant  load 
maintenance  was  obtained.  This  control  also  provided  means  for  rapidly  changing 
from  one  load  level  to  another  and  finally  for  releasing  loads  rapidly  without 
introducing  undue  shock  and  inertial  effects.  Figure  40  shows  a  schematic  repre¬ 
sentation  for  the  constant  temperature-constant  load  test  program  and  Figure  129 
shows  a  similar  representation  for  a  pulsed  step-load  test  program. 

Figure  130  displays  and  compares  the  results  obtained  from  the  above  tests.  Those 
tests  run  at  70KSI  load  and  1000°F  temperature  resulted  in  the  steeply  rising 
creep  curves-identif ied  Mark  A  and  Mark  B  in  Figure  130.  These  particular  tests 
were  not  run  for  the  full  four  hours  planned,  but  they  were  run  out  to  the  limit 
of  extension  measurable  by  the  equipment  used.  The  tests  run  at  50KSI  load  and 
1000 °F  resulted  in  the  flat  creep  curves  identified  Mark  C  and  D  in  Figure  130. 
Those  creep  curves  identified  Mark  E  and  Mark  F  in  Figure  13 0  show  creep  curves 
resulting  from  the  pulsed  step-load  tests  also  run  at  1000  F  under  the  loading 
conditions  described  above. 


In  comparing  results  from  "continuous"  creep  with  those  obtained  from  "pulsed" 
creep,  the  coordinates  representing  the  terminal  observations  from  the  two  tests 
run  at  70KSI  and  1000°F  (Mark  A  and  Mark  B)  were  added  to  the  terminal  observation 
obtained  in  the  first  50KSI,  1000°F  creep  test  (Mark  C).  This  resulted  in  locating 
the  summation  points  shown  by  Mark  G  and  Mark  H.  Likewise  the  coordinates ^repre¬ 
senting  the  terminal  observations  from  the  two  tests  run  at  fOKSI  and  1000  F 
(Mark  A  and  Mark  B)  were  added  to  the  terminal  observations  obtained  in  the  second 
50KSI,  1000°F  creep  test  (Mark  D).  This  resulted  in  locating  the  summation  points 

shown  by  Mark  J  and  K . 


In  consideration  of  the  data  shown  in  Figure  130  the  following  observations  can  be 
made: 

A.  The  first  (Mark  C)  50KSI,  1000°F  constant  stress  test  acquired  0.15$ 
strain  in  245  minutes  and  the  first  (Mark  A)  7CKSI,  1000  F  constant 
stress  test  acquired  3-15$  strain  in  100  minutes.  These  outcomes  sum 
up  to  3.3O $  strain  acquired  in  345  minutes  by  constant  stress  loadings. 
This  compares  with  the  lesser  2.35$  strain  acquired  in  355  minutes 
through  pulsed  step-loading. 


B. 


C. 


D. 


The  first  (Mark  c)  50KSI,  1000°F  constant  stress  test  acquired  0.15$ 
strain  in  245  minutes  and  the  second  (Mark  B)  70KSI,  1000°F  constant 
stress  test  acquired  3.56$  strain  in  13 0  minutes.  These  outcomes  sum 
up  to  3.71$  strain  acquired  in  385  minutes  by  constant  stress  loadings. 
This  compares  with  the  lesser  2.65$  strain  acquired  in  385  minutes 
through  pulsed  st ep- loadings . 


The  second  (Mark  D)  50SKI,  1000°F  constant  stress  test  acquired  0.54$ 
strain  in  255  minutes  and  the  first  (Mark  A)  70KSI,  1000°F  constant 
stress  test  acquired  3-15$  strain  in  100  minutes.  These  outcomes  sum 
up  to  3.69$  strain  in  355  minutes  by  constant  stress  loadings.  This 
compares  with  the  lesser  2.53$  strain  acquired  in  355  minutes  through 
pulsed  step-loading. 


The  second  (Mark  D)  50KSI,  1000°F  constant  stress  test  acquired  0.54$ 
strain  in  255  minutes  and  the  second  Mark  B  70KSI,  1000°F  constant 
stress  test  acquired  3-56$  strain  in  13 0  minutes.  These  outcomes  sum 
up  to  4.10$  strain  acquired  in  385  minutes  by  constant  stress  loadings. 
This  compares  with  the  lesser  2.76$  stress  obtained  in  385  minutes 
through  pulsed  step-loading. 
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Figure  131  compares  the  several  numbers  obtained  by  slimming  outcomes  of  constant 
stress  tests  and  comparing  them  with  outcomes  of  pulsed  step-load  tests.  The 
suggested  trend  lines  intimate  that  constant  stress  loading  induces  greater  amounts 
of  creep  than  pulsed  step-loading  because  the  strain  it  produces  accumulates  more 
rapidly.  These  trend  lines  bring  focus  on  the  relaxation  effects  observed  in 
preceding  test  studies  wherein  both  a  "work-hardening"  effect  which  tended  to 
reflect  in  increasing  amounts  of  elastic  action  and  in  a  partition  of  loading  input 
energy  into  elastic  recovery  and  permanent  set  on  strain  was  observed.  In  view 
of  the  unrefined  test  techniques  used  together  with  the  high  stress  levels  and 
temperatures  also  used  to  bring  focus  on  phenomenon  trends,  a  need  for  critical 
and  refined  experiment  is  expressed  in  the  foregoing  outcome.  This  approach  bears 
consideration  because,  on  one  hand,  evaluations  of  "material  constants"  necessary 
to  the  practical  operation  of  creep  prediction  procedures  will  of  necessity  need 
be  gathered  rapidly  under  simplified  experimental  approaches,  and,  on  the  other 
hand,  the  cyclic  nature  of  aircraft  loadings  require  realistic  assessments  to 
permit  correlations  of  "materials  evaluation  test  data"  with  actualities  of  material 
performance  under  flight  service  conditions  in  order  to  render  utility  into  the 
prediction  schemes. 

SUMMARY  OF  TEST  FINDINGS 


Those  tests  which  isolated  the  stress  variable  for  study  repeatedly  asserted  the 
importance  of  initial  load  application  in  establishing  deformation  patterns.  The 
deformation  attained  upon  initial  load  application  is  comprised  of  elastic  and 
inelastic  or  plastic  components.  The  elastic  component,  within  limits,  operates 
under  conditions  of  reversibility  whereas  the  inelastic  component  does  not.  The 
reversible  action  appears  related  in  extent  and  perhaps  is  limited  to  the  elastic 
limit  established  by  the  temperature  at  which  the  material  is  operated.  The 
plastic  or  irreversible  action  taking  place  in  initial  loading  appeared  to  be 
somewhat  time  dependent  since  loading  rates  appeared  to  influence  its  extent. 

Beyond  the  initial  loading,  the  extent  of  the  deformation,  which  is  inelastic  in 
nature,  is  dependent  upon  the  temperature  and  stress  at  which  loads  are  applied 
for  the  establishment  of  the  rate  of  accumulation  of  irreversible  deformation,  and, 
the  times  for  which  these  rates  are  allowed  to  accumulate. 

Those  tests  which  isolated  the  temperature  variable  for  study  asserted  the  com¬ 
plexity  of  events  taking  place  upon  heating  to  temperatures  and  a  need  for  under¬ 
standing  of  these  so  that  the  extent  of  reversible  action  taking  place  could  be 
assessed.  The  plastic  action  taking  place  upon  attainment  of  temperature  appeared 
to  be  dependent  also  on  the  temperature  and  stress  at  which  loads  are  applied  for 
establishment  of  rate  of  the  accumulation  of  the  irreversible  deformation,  and,  the 
times  for  which  these  rates  are  allowed  to  accumulate. 

The  test  findings  assert  need  for  critical  or  precise  study  of  (1)  the  effect  of 
load  rate  in  initial  load  applications  upon  the  partition  of  deformation  action 
into  elastic  and  inelastic  parts  especially  with  respect  to  the  temperature  present 
and  the  elastic  limit  of  given  materials;  (2)  the  interrelations  existing  between 
times  and  strains  appearing  with  varying  stress  levels  at  given  temperatures  as  a 
result  of  accumulation  of  plastic  movements;  (3)  the  interrelations  existing 
between  times  and  temperatures  required  to  produce  equivalent  strains  at  specific 
stress  levels;  and  (4)  the  combining  of  the  cumulative  distribution  of  strain 
function  (step  2)  with  the  exponential  growth  pattern  (step  3)  to  produce  a  tool 
for  prediction  of  inelastic  or  plastic  action  increase  .  This  in  combination  with 
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considered  relations  derived  from  loading  rate  relations  can  result  in  a  compre¬ 
hensive  basis  for  predicting  basic  material  deformation  patterns  under  conditions 
of  elastic  and  inelastic  behavior  anticipated  with  varying  loading  rates,  stress 
levels  of  operation,  temperature  and  times  of  use. 
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TABLE  40  CONDITIONS  OF  TEST  AND  PERTINENT  DATA  POINT  VALUES 
FOUND  IN  "STANDARD"  CREEP  TESTS 


TEMP 

°F 
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ft 

tt 
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ft 

It 
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tf 
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DEFORMATION-PERCENT 


MATERIAL:  ARM  CO  PH  15-7  MO-RH  950 

TEMPERATURE:  700°  F 
LOAD:  180  KSI. 
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Figure  40.  Time-Deformation  Data,  Armco  PH  15-7  Mo  700°F. , 
180  KSI,  "Standard"  Creep  Test 
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DEFORMATION-PERCENT 


MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  700°  F 
LOAD:  180  KSI. 
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Figure  41.  Time- Deformation  Data,  Armco  PH  15-7  Mo  700°F. , 
180  KSI,  "Standard"  Creep  Test 
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MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  700°  F 

LOAD:  182  KSI. 
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Figure  42.  Time- Deformation  Data,  Armco  PH  15-7  Mo  700°F. ,  182  KSI,  "Standard"  Creep  Test 
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MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  700°  F 
LOAD:  182  KSI. 


Figure  43.  Time-Deformation  Data,  Armco  PH  15-7  Mo  700°F. , 
182  KSI,  "Standard"  Creep  Test 
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DEFORMATION-PERCENT 


MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  800°  F 


Figure  44.  Time-Deiormation  Data,  Armco  PH  15-7  Mo  800°F. , 
170  KSI,  "Standard"  Creep  Test 
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DEFORMATION-PERCENT 


3. 


MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  800°  F 
LOAD:  170  KSI. 
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Figure  45.  Time- Deformation  Data,  Armco  Ph  15-7  Mp  800°F. , 
170  KSI,  "Standard"  Creep  Test 
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DEFORMATION-PERCENT 
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MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  800°  F 

LOAD:  175  KSI. 
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Figure  46.  Time- Deformation  Data,  Armco  PH  15-7  Mo  800°F. , 
175  KSI,  "Standard"  Creep  Test 
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DEFORMATION-PERCENT 
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MATERIAL: 

TEMPERATURE: 

LOAD: 


ARMCO  PH  15-7  MO-RH  950 
800°  F 
175  KSI. 


Figure  47.  Time- Deformation  Data,  Armco  PH  15-7  Mo  800°F 
175  KSI,  "Standard"  Creep  Test  *  ’ 


WADD  TR  60-411  Part  I 


151 


DEFORMATION-PERCENT 
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MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  900°  F 

LOAD:  130  KSI. 
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Figure  48.  Time-Deformation  Data,  Armco  PH  15-7  Mo  900°F. , 
130  KSI,  "Standard"  Creep  Test 
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Figure  49.  Time-Deformation  Data,  Armco  PH  15-7  Mo  900°F. ,  130  KSI,  "Standard"  Creep  Test 
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Figure  50.  Time-Deformation  Data,  Armco  PH  15-7  Mo  900°F. ,  135  KSI,  "Standard"  Creep  Test 
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Figure  51.  Time-Deformation  Data,  Armco  PH  15-7  Mo  900°F. ,  135  KSI,  "Standard"  Creep  Test 
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MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  1000°  F 

LOAD:  55  KSI. 


Figure  52.  Time- Deformation  Data,  Armco  PH  15-7  Mo  1000° F. , 
55  KSI  "Standard”  Creep  Test 
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Figure  53.  Time-Deformation  Data,  Armco  PH  15-7  Mo  1000°F. ,  55  KSI,  "Standard"  Creep  Test 
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MATERIAL:  ARM  CO  PH  15-7  MO-RH  950 

TEMPERATURE:  1000°  F 
LOAD:  77  KSI. 
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Figure  54.  Time- Deformation  Data,  Armco  pH  15-7  Mo  1000° F. , 
77  KSI  "Standard-’  Creep  Test 
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DEFORMATION-PERCENT 
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MATERIAL:  ARMCO  PH  15-7  MO-RH  950 

TEMPERATURE:  1000°  F 

LOAD:  77  KSI. 


Figure  55.  Time- Deformation  Data,  Armco  PH  15-7  Mo  1000°F.  , 
77  KSI  "Standard"  Creep  Test 
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Figure  56.  "Make-up"  of  Loading  Strains  in  "Standard”  Creep  Tests  with  Armco  PH  15-7  Mo. 
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LOAD  180  KSI.  I  I  182  KSI.  170  KSI.  I  1 175  KSI.  130  KSI.  I  I  135  KSI.  55  KSI.  I  I  77 

TEMPI* - 700°  F - H  - - 800°  F - i  - - 900°  F' - —I  I— —  1000°F- 

Figure  57.  "Make-up"  of  Unloading  Strains  in  "Standard"  Creep  Tests  with  Armco  PH  15-7  Mo. 
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TOTAL  TIME  INDEPENDENT  DEFORMATION 
(ELASTIC  &  INELASTIC)  UPON  LOADING 
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Figure  59.  Comparison  of  Time  Independent  Strains  Observed  upon  Loading  and  Unloading  "Standard 
Creep  Tests  with  Armco  PH  15-7  Mo  at  900°F  and  1000°F. 
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Figure  60.  Comparison  of  TotaW^lastic  Recovery  Observed  upon  Unloading  "Standard"  Creep  Tests  with 
Armco  PH  15-7  Mo  at  700°F  and  800°F  with  the  Input  Strain  Converted  into  Permanent  Set  up- 
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Figure  61.  Comparison  of  Total  Elastic  Recovery  Observed  upon  Unloading  "Standard”  Creep  Tests  with 
Armco  PH  15-7  Mo  at  900°F  and  1000°F  with  the  Input  Strain  Converted  into  Permanent  Set  up¬ 
on  Loading 
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Figure  62.  Typical  Field  of  Deformation  Activity  Observed  upon  Increment¬ 
al  Step  Loading  Armco  PH  15-7  Mo  at  Elevated  Temperature 
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Figure  63.  Typical  Field  of  Deformation  Activity  Observed  upon  Unloading 
Incremental  Step  Loaded  Armco  PH  15-7  Mo  at  Elevated  temp¬ 
eratures 
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Figure  64.  Program  for  Incremental  Step  Loaded  Test  (Schematic) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950 

TEMPERATURE:  700°  F 

FULL  LOAD:  185  KSI. 

o  -LOAD  POINT  DLOAD  RELEASE  POINT 
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Figure  65.  Time-Deformation  Data,  Armco  PH  15-7  Mo  700°F. ,  185  KSI  Incremental  Step  Loading. 
(Loading  Onuifon) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950 

TEMPERATURE:  700°  F 

FULL  LOAD:  185  KSI. 

o  LOAD  POINT  □  LOAD  RELEASE  POINT 
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Figure  66.  Time  Deformation  Data,  Armco  PH  15-7  Mo  700°F.  ,  185  KSI,  Incremental  Step  Loading 
(Loading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950  FULL  LOAD:  185  KSI 
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Figure  69.  Time-Deformation  Data,  Armco  PH  15-7  Mo  800°F. ,  175  KSI„  Incremental  Step  Loading 
(Unloading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950 

TEMPERATURE:  900°  F 

FULL  LOAD:  135  KSI. 

o  LOAD  POINT  □  LOAD  RELEASE  POINT 
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Figure  7(X  Time-Deformation  Data,  Armco  PH  15-7  Mo  900°F. ,  135  KSL  Incremental  Step 
(Loading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950 

TEMPERATURE:  900°  F 

FULL  LOAD:  135  KSI. 

o  LOAD  POINT  □  LOAD  RELEASE  POINT 
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Figure  7L  Time- Deformation  Data,  Armco  PH  15-7  Mo  900°F.,  135  KSL  Incremental  Step  Loading 
(Loading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950  FULL  LOAD:  135  KSI. 

TEMPERATURE:  900°  F 
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L0AD-(SEE  NOTE) 

Figure  72.  Time-Deformation  Data,  Armco  PH  15-7  Mo  900°F. ,  135  KSI.  Incremental  Step  Loading 
(Unloading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950 

TEMPERATURE:  1000°  F 

FULL  LOAD:  77  KSI. 

o  LOAD  POINT  □  LOAD  RELEASE  POINT 
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Figure  73.  Time-Deformation  Data,  Armco  PH  15-7  Mo  1000°F. ,  77  KSL  Incremental  Step  Loading 
(Loading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950 

TEMPERATURE:  1000°  F 

FULL  LOAD:  77  KSI. 

o  -LOAD  POINT  DLOAD  RELEASE  POINT 
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Figure  74.  Time-Deformation  Data,  Armco  PH  15-7  Mo  1000°F. ,  77  KSL  Incremental  Step  Loading 
(Loading  Condition) 


MATERIAL:  ARMCO  PH  15-7  MO.  COND.  RH  950  FULL  LOAD:  77  KSI. 

TEMPERATURE:  1000°  F 
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L0AD-(SEE  NOTE) 

Figure  75.  Time-Deformation  Data,  Armco  PH  15-7  Mo  1000°F. ,  77  KSL  Incremental  Step  Loading 
(Unloading  Condition) 
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Cycle  5  rotation  indicates  sum  of  action 
observed  In  cycles  2,  3,  4  and  5  only. 

See  loading  program,  Figure  64. 

Material!  Armco  PH  15*7  Mo 
Condition;  R  H950  Steel 


700”  F  II  800°  F  _ I  I  900°  F  I  |  1000°  F _ I 
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182  KSI 
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"Make-up"  of  Loading  Strains  in  Increment  Loaded  Creep 
Tests  with  Armco  PH15-7  Mo. 
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'NOTE:  1st  fast  and  2nd  fast  indicates  result  obtained 
from  a  single  test  specimen. 

Cycle  5  rotation  Indicates  sum  of  action 
observed  in  cycles  2,  3,  4  and  5  only. 

Sea  loading  program,  Figure  64. 

Material:  Armco  PH  15-7  Mo 
Condition:  RH  950  Steel 
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"Make-up"  of  Unloading  Strains  in  Increment  Loaded 
Creep  Tests  with  Armco  PH  15 -7  Mo. 
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OF  LOADING  STRAIN 
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Figure  78.  Comparison  of  Strains  Incidental  to  Loading  With  Elastic 

Strain  Recovery  Occurring  Upon  Load  Release  (700°F  and  800°F). 
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Figure  79.  Comparison  of  Strain  Incidental  to  Loading  with  Elastic 

Strain  Recovery  Occurring  Upon  Load  Release  (900°F  &  1000°F) 
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Figure  80.  Program  for  Repeated  Load  Tests. 


STRAIN-PERCENT 


TIME-MINUTES 

Figure  81.  Strains  Resulting  From  Successive  Applications  of  182 
KSI  Load  to  Armco  PH15-7  Mo  at  700°F. 
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NOTE:  SEE  LOADING  PROGRAM,  FIGURE  80 
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Figure  82.  Strains  Resulting  from  Successive  Applications  of  175  KSI 
Load  to  Armco  PH  15-7  Mo  at  800°F. 


WADD  TR  60-411  Part  I 


186 


STRAIN-PERCENT 


Figure  83.  Strains  Resulting  from  Successive  Applications  of  135  KSI 
Load  to  Armco  PH  15-7  Mo  at  900°  F. 


WADD  TR  60-411  Part  I 


187 


STRAIN-PERCENT 
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NOTE:  SEE  LOADING  PROGRAM,  FIGURE  80 

MATERIAL:  ARMCO  PH  15-7  MO,  CONDITION  RH  950  STEEL 
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Figure  84.  Strains  Resulting  From  Successive  Applications  of  77  KSI 
Load  to  Armco  PH  15-7  Mo  at  1000°F, 
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Figure  85.  “Make-up”  of  Strain  Patterns  Occurring  Upon  Load  Re¬ 
lease  in  Repeated  Loading  Tests  (700°F  182  KSI  Stress). 
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Figure  86.  "Make-up”  of  Strain  Patterns  Occurring  Upon  Load  Re¬ 
lease  in  Repeated  Loading  Tests  (800°F  175  KSI  Stress). 
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Figure  87.  Make-up"  of  Strain  Patterns  Occurring  Upon  Load  Re¬ 
lease  in  Repeated  Loading  Tests  (900°F  135  KSI). 


WADD  TR  60-411  Part  I 


191 


STRAIN-PERCENT 


Figure  88.  “Make-up”  of  Strain  Patterns  Occurring  Upon  Load  Re¬ 
lease  in  Repeated  Loading  Tests  (1000  F  77  KSI). 
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Figure  90.  Strain  Changes  Incidental  to  Repeated  Load  Tests 
(800°F,  175  KSI). 
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Figure  91.  Strain  Changes  Incidental  to  Repeated  Load  Tests 
(900 °F,  135  KSI). 
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Figure  92.  Strain  Changes  Incidental  to  Repeated  Load  Tests 
(1000°  F,  77  KSI). 
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Figure  93.  Permanent  Set  Accumulation  Observed  in  Repeated 
Load  Tests  Run  at  Various  Temperatures. 
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LOAD  APPLICATION  SEQUENCE 

Figure  94.  Elastic  Recovery  Accumulation  Observed  in  Repeated 
Load  Tests  Run  at  Various  Temperatures. 
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C0NDI1I0N  OF  TEST:  TRANSIENT  TEMPERATURE  APPLICATIONS 

AND  CONSTANT  LOADING 


2.2 


HEATING  TIME- MINUTES 

Figure  96.  Initial  Cycle,  Transient  Heating  Creep  Curves. 
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O  FIRST  SPECIMEN 
□  SECOND  SPECIMEN 

SEE  LOADING  PROGRAM,  FIGURE  95 


Figure  97.  Initial  Cycle,  Transient  Heating  Creep  Curves 
(Logarithmic  700°  and  800°F). 
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o  FIRST  SPECIMEN 
□  SECOND  SPECIMEN 

SEE  LOADING  PROGRAM  FIGURE  95 
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Figure  98.  Initial  Cycle  Transient  Heating  Creep  Curves 
(Logarithmic  900°  and  1000°  F) 
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Figure  99.  Strain  Transition  -  Temperature  Relations  in  Initial 
Transient  Heating  Cycles. 
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Figure  100.  Initial  Cycle  Transient  Heating  Creep  Curves  for 
Constant  Temperature, Varying  Stress  Conditions. 
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o  FIRST  SPECIMEN 
□  SECOND  SPECIMEN 

SEE  LOADING  PROGRAM  FIGURE  95 


Figure  101„  Initial  Cycle  Transient  Heating  Creep  Curves  for  Constant 

Stress,  Varying  Temperature  Conditions  (700°F,  182  &  185  KSI)„ 
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Figure  102 . 


O  FIRST  SPECIMEN 
□  SECOND  SPECIMEN 
SEE  LOADING  PROGRAM  FIGURE  95 


Initial  Cycle  Transient  Heating  Creep  Curves  for  Constant 
Stress,  Varying  Temperature  Conditions  (700°F,  188  &  191  KSI). 
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Figure  103.  Strain  Transition  -  Temperature  Relations  in  Initial 
Transient  Heating  Cycles. 
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Figure  105.  Transient  Heating  Creep  Curves  (700°F,  185  KSI) 
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Figure  106,  Transient  Heating  Creep  Curves  (700°F,  188  KSI) 
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Figure  107.  Transient  Heating  Creep  Curves  (700°F,  191  KSI) 
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Figure  108,  Transient  Heating  Creep  Curves  (800°F,  175  KSI) 
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Figure  109„  Transient  Heating  Creep  Curves  (900°F,  135  KSI) 
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Figure  110o  Transient  Heating  Creep  Curves  (1000°F,  77  KSI) 
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Figure  111.  Program  for  Pulsed  Heating  Tests. 
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Figure  112.  Pulsed  Heating  Creep  Curves  (700°F). 
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Figure  113.  Pulsed  Heating  Creep  Curves  (800°F). 
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Figure  114.  Pulsed  Heating  Creep  Curves  (900°F). 
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TIME  AT  TEMPERATURE-MINUTES 
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Figure  116.  Effect  of  Stress  Variation  at  700°F  Upon  5  Minute  Pulsed  Heating  Creep. 
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Figure  117.  Effect  of  Stress  Variation  at  900°F  Upon  5  Minute  Pulsed  Heating  Creep. 
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Figure  118.  Effect  of  Stress  Variation  at  1000°F  Upon  5  Minute  Pulsed  Heating  Creep. 
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Figure  119.  Effect  of  Stress  Variation  at  1000°F  Upon  Pulsed  Heating  Creep. 
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Figure  120.  Effect  of  Stress  Variation  at  700°F  Upon  125  Minute  Pulsed  Heating  Creep. 
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Figure  121.  Effect  of  Stress  Variation  at  900°F  Upon  125  Minute  Pulsed  Heating  Creep. 
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Figure  122.  Effect  of  Stress  Variation  at  1000°F  Upon  Pulsed  Heating  Creep  at  Various 
Stress  Levels. 
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Figure  123.  Effect  of  Pulse  Time  Variation  at  Constant  Stress  at  700°F  on  Pulsed  Creep  (180  KSI). 
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Figure  124.  Effect  of  Pulse  Time  Variation  at  Constant  Stress  at  700°F  on  Pulsed  Creep  (185  KSI). 
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Figure  125.  Effect  of  Pulse  Time  Variation  at  Constant  Stress  at  1000°F  on  Pulsed 
Heating  Creep  (77  KSI). 
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Figure  126.  Effect  of  Pulse  Time  Variation  at  Constant  Stress  at  1000°F  on  Pulsed 
Heating  Creep  (60  KSI). 
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Figure  127.  Effect  of  Stress  and  Temperature  Variation  at  5  Minute  Pulse  Intervals 
Upon  Creep. 
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Figure  12«.  Effect  of  Stress  and  Temperature  Variation  at  125  Minute  Pulse  Intervals 
Upon  .Creep. 
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Figure  129.  Program  for  Pulsed  Step  Load  Tests. 
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Figure  130.  Comparisons  of  Creep  Resulting  From  Constant  Stress  -  Constant  Temperature 
and  Pulsed  Step  Loading  -  Constant  Temperature  Applications. 
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Figure  131.  Comparison  of  Creep  Trends  in  Constant  Stress 
and  Pulsed  Step  Load  Tests. 
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PREDICTION  OF  CREEP  EFFECTS 
IN 

AIRCRAFT  STRUCTURES 


APPENDIX  II  -  MATERIAL  DATA 


TABLE  4l  EFFECT  OF  VARIOUS  TEMPERATURES  ON  THE  MODULUS 
OF  ELASTICITY  OF  ALUMINUM  ALLOYS 


TEMPERATURE 

°F 

APPROXIMATE  VALUE  OF  MODULUS  OF  ELASTICITY  IN  TERMS 

OF  75 °F  VALUES,  (l) 

75S 

Alloy 

i4s,  17s,  24s,  25s,  a 51s,  52s 

and  6lS  Alloys 

-  320 

112  $ 

112$ 

-  112 

10 5$ 

105  $ 

-  18 

102$ 

102$ 

75 

100  $ 

100$ 

212 

95$ 

98$ 

300 

00 

Oo 

-cR 

95$ 

400 

80$ 

90$ 

500 

■eR- 

00 

VO 

00 

0 

-cR. 

6oo 

50$ 

70$ 

NOTE:  (1) 

Ref.  ANC-5, 

March  1955* 
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TABLE  42  AVERAGE  COEFFICIENT  OF  THERMAL  EXPANSION  FOR  SOME 

ALUMINUM  ALLOYS 


ALLOY  AVERAGE  COEFFICIENT 

OF  THERMAL  EXPANSION 
FAHR.  (1) 

IN.  PER  DEGREE 

68°  -  212  °F 

TEMPERATURE  RANGE 
68°  -  392  °F 

68°  -  572 °F 

l4s 

0 . 0000122 

0.0000130 

O.OOOOI38 

17S 

0 . 0000122 

0.0000130 

0 . 0000138 

25S 

0.0000122 

0.0000130 

0.0000138 

A51S 

0 . 0000111 

0 . 0000113 

0.0000119 

52S 

0.0000111 

0.0000113 

0.0000119 

6lS 

0.0000133 

0 . 0000138 

0 . ooooi44 

75S 

0.0000129 

0.0000135 

o.ooooi44 

NOTE:  (1) 

Ref.  "Alcoa  Alloys",  1944. 

TABLE  43 

AVERAGE  COEFFICIENT  OF  THERMAL  EXPANSION  FOR  VARIOUS 

METALS 

ALLOY 

AVERAGE  COEFFICIENT  OF  THERMAL  EXPANSION  IN.  PER  DEGREE 

FAHR.  (1) 

TEMPERATURE  RANGE 

68°  -  212 °F . 

Brass 

0.0000097  to  0.0000107 

Cast  Iron 

0.0000059 

Copper 

O.OOOOO93 

Lead 

0.0000170 

Monel 

O.OOOOO78 

Nickel 

0.0000072 

Steel 

0.0000070 

Zinc 

0.0000180 

NOTE:  (1)  Ref.  "Alcoa  Alloys"  1944. 
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